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Unclassified 

StCUSITT  CLASSIFICATION  OF  THIS  FAC t 


FOREWORD 


by 

Dwain  K.  Butler 

Flow  of  fluids  through  a  porous  soil  on  rock  generates  electrical  volt¬ 
ages  (potentials)  through  a  process  known  as  electrokinesis.  These  potentials 
are  called  streaming  or  self  potentials  (SP) .  The  magnitude  of  the  SP  depends 
on  the  electrical  resistivity,  dielectric  constant  and  viscosity  of  the  fluid, 
on  a  coupling  constant  between  the  fluid  and  the  soil/rock,  and  on  the  pres¬ 
sure  drop  along  the  flow  path.  The  SP  anomaly  caused  by  the  flow  can  be  mea¬ 
sured  on  the  surface  above  the  flow  path;  this  is  the  basis  of  the  SP  method 
for  seepage  detection  and  mapping. 

The  US  Army  Engineer  Waterways  Experiment  Station  (WES)  has  successfully 
applied  the  self-potential  (SP)  and  other  geophysical  methods  to  detect,  map, 
and  monitor  anomalous  seepage  conditions  at  water  retention  and  hazardous 
waste  disposal  sites  throughout  the  United  States.  The  keystone  of  this  suc¬ 
cessful  methodology  has  been  the  self-potential  method,  which  has  been  applied 
using  permanent  arrays  of  inexpensive  copper-clad  steel  electrodes  (cut  from 
common  grounding  rod  stock).  Use  of  the  metallic  electrodes  for  SP  measure¬ 
ments  is  contrary  to  commonly  accepted  geophysical  practice;  however,  cost, 
ease  of  installation  and  maintenance,  and  general  success  considerations 
seemed  to  outweigh  other  factors.  Data  processing  techniques  were  developed 
which  attempted  to  compensate  for  electrode  polarization  effects  and  sensi¬ 
tivity  to  environmental  variables  (see  References  below).  In  spite  of  the 
general  success  of  the  methodology,  cases  were  encountered  where  the  data  were 
extremely  noisy,  and  straightforward  interpretation  was  not  possible.  In 
recent  years,  partly  as  a  result  of  WES  successes,  geotechnical  applications 
of  the  SP  method  have  increased  in  the  United  States,  and  nonpolarizing  elec¬ 
trodes  which  are  more  rugged  and  maintenance-free  have  been  developed  and 
applied.  Although  there  has  always  been  the  promise  for  quantitative  inter¬ 
pretation  of  the  SP  data  to  give  flow  rates  and  depths,  data  quality  has  gen¬ 
erally  not  permitted  it,  and  the  primary  need  was  for  mapping  flow  paths  in 
plan. 
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This  report  is  pivotal  in  the  development  of  geotechnical  applications 
of  the  SP  method  for  anomalous  seepage  detection,  mapping,  and  monitoring. 

Many  of  the  environmental  factors  which  can  affect  SP  measurements  with  metal¬ 
lic  electrodes  are  compared  in  detail  with  the  general  lack  of  effect  of  the 
same  factors  on  nonpolarizing  electrodes.  The  application  of  nonpolarizing 
electrodes  at  Beaver  Dam,  Arkansas,  "side-by-side"  with  a  metallic  electrode 
monitoring  network  illustrates  that  the  metallic  electrode  data  are  not  just 
shifted  up  or  down  in  magnitude,  while  showing  qualitative  agreement,  as  has 
been  postulated  earlier.  In  many  locations,  the  metallic  electrodes  are 
clearly  responding  to  phenomena  not  detected  at  all  by  the  nonpolarizing 
electrodes.  Also,  nonpolarizing  electrodes  now  exist  which  can  be  installed 
in  a  monitoring  network  that  will  require  no  maintenance  over  a  several  year 
period.  Thus,  the  advantages  of  nonpolarizing  electrodes  seem  to  outweigh  the 
lower  cost  and  ease  of  installation  of  metallic  electrodes  for  seepage 
monitoring. 

Appendix  B  presents  a  summary  of  recommended  field  procedures  and  data 
reduction  methods  for  SP  surveys  with  nonpolarizing  electrodes.  This  proce¬ 
dure  is  for  a  single  pass  survey,  and  monitoring  applications  would  require 
repetition  of  the  survey  at  later  times.  Details  of  field  procedure  for 
long-term  monitoring  permanent  arrays  of  nonpolarizing  electrodes  have  been 
developed  and  proposed  for  application  at  the  Beaver  Dam  site. 

Also,  this  report  presents  initial  attempts  to  model  quantitatively  SP 
survey  results  acquired  at  a  Corps  damsite.  The  computer  program  documented 
in  this  report  utilizes  geometric  SP  source  models,  and  data  can  be  inter¬ 
preted  by  iterative  adjustment  of  source  model  parameters.  These  geometric 
source  models  are  relatively  simple  compared  with  quantitative  modeling  which 
treats  SP  from  first  principle  considerations  of  streaming  potential 
generation  from  flow  through  a  porous  medium.  Computer  programs  which  model 
cross-coupled  phenomena  in  a  two-  or  three-dimensional  medium  are  complex  to 
use  and  require  considerable  computer  memory  and  execution  time;  efforts  are 
underway  to  streamline  and  simplify  these  sophisticated  computer  programs. 

This  report  presents  the  results  of  a  preliminary  run  of  the  quantitative 
modeling  program  on  a  "Beaver  Dam-like"  problem.  More  detailed  results  of 
modeling  the  Beaver  Dam  SP  data  were  presented  at  an  International  Symposium 
on  Detection  of  Subsurface  Flow  Phenomena  at  Karlsruhe,  Germany. 
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PREFACE 


The  work  described  herein  was  performed  during  the  period  August  1986  to 
February  1988  by  Dr.  Robert  W.  Corwin,  Consulting  Geophysicist  to  Harding- 
Lawson  Associates,  Navato,  California.  The  work  was  performed  under  Contract 
No.  DACW39-86-C-0059 ,  and  funded  by  the  Repair,  Evaluation,  Maintenance,  and 
Rehabilitation  (REMR)  Research  Project  "Geophysical  Techniques  for  Assessment 
of  Existing  Structures  and  Structural  Foundations,"  Work  Unit  No.  WU  1-32315. 

This  report  was  prepared  by  Dr.  Corwin.  Principal  Investigator  for  the 
REMR  Research  Project  was  Dr.  Dwain  K.  Butler,  Earthquake  Engineering  and 
Geophysics  Division  (EEGD) ,  Geotechnical  Laboratory  (GL) ,  WES.  Dr.  Butler 
also  prepared  the  Foreword  and  conducted  an  extensive  technical  review  of  the 
report.  The  field  work  was  performed  at  Beaver  Dam,  Arkansas,  in  close 
coordination  with  other  work  at  the  site  performed  by  EEGD,  Little  Rock 
District,  Corps  of  Engineers,  and  other  contractors.  The  report  was  prepared 
under  the  general  supervision  of  Dr.  Arley  G.  Franklin,  Chief,  EEGD,  and 
Dr.  William  F.  Marcuson,  Chief,  GL.  Senior  technical  reviewer  of  this  report 
was  Dr.  Paul  F.  Hadala,  Assistant  Chief,  GL.  Mrs.  Joyce  H.  Walker,  Informa¬ 
tion  Products  Division,  Information  Technology  Laboratory,  edited  the  report. 

COL  Dwayne  G.  Lee,  EN,  is  Commander  and  Director  of  WES.  Dr.  Robert  W. 
Whalin  is  Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 


lx 


To  Obtain 


atmospheres  (standard) 
degrees  (angle) 
Fahrenheit  degrees 


101.325 

0.01745329 

5/9 


feet 

inches 

pounds  (mass) 


0.3048 

2.54 

0.4536 


kilopascals 

radians 

Celsius  degrees 
or  Kelvins* 

metres 

centimetres 

kilograms 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C  =  (5/9) (F  -  32).  To  obtain  Kelvin  (K)  read¬ 
ings,  use:  K  =  (5/9) (F  -  32)  +  273.15. 
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DEVELOPMENT  OF  SELF-POTENTIAL  INTERPRETATION 


TECHNIQUES  FOR  SEEPAGE  DETECTION 


PART  I:  INTRODUCTION 


1.  This  report  describes  the  results  of  a  research  program  conducted  to 
develop  and  document  field  techniques,  analytical  methods,  and  computer  pro¬ 
grams  to  improve  data  acquisition  and  interpretation  procedures  for  the  use  of 
the  self-potential  (SP)  method  for  dam  seepage  investigations.  As  the  third 
in  a  series  of  reports  on  geotechnical  applications  of  the  SP  method  (Erchul 
1988;  Erchul  and  Slifer  1989),  this  report  addresses  key  fundamental  questions 
related  to  past,  present,  and  future  practice.  The  key  questions  relate  to 
(a)  electrode  comparisons  and  long-term  stability  concerns;  (b)  data  acquisi¬ 
tion  procedures;  and  (c)  modeling  and  quantitative  interpretation  of  SP  data. 

2.  The  research  program  included  four  major  components: 

a.  Studies  of  electrodes  to  determine  their  suitability  and  charac¬ 
teristics  for  long-term  SP  monitoring  (Part  II). 

b.  Field  investigations  conducted  as  part  of  a  large-scale  geophy¬ 
sical  study  at  Beaver  Lake  Dam,  Arkansas  (Part  III). 

c.  Development  of  computer  programs  for  interpretation  of  SP  data 
(Part  IV  and  Appendix  A) . 

d.  Compilation  of  a  bibliographic  data  base  of  publications  relat¬ 
ing  to  SP  methods  for  dam  seepage  investigations  (Part  VI) . 

The  results  developed  for  each  of  these  components  are  described  in  the  indi¬ 
cated  parts  of  this  report.  A  summary  of  the  findings  of  this  study  and 
recommended  procedures  is  presented  in  Part  V.  Appendix  B  is  a  stand-alone 
guide  to  SP  field  procedures  and  data  reduction.  Appendix  C  is  a  reprint  of  a 
paper  which  gives  details  of  the  Beaver  Dam,  Arkansas,  site,  as  a  reference 
for  Part  III. 


7 


PART  II:  ELECTRODE  STUDIES 


Introduction 


3.  The  objective  of  this  component  of  the  research  program  was  to 
determine  the  suitability  of  various  electrode  types  for  long-term  monitoring 
of  SP  variations  related  to  seepage  flow,  and  to  establish  the  responses  of 
these  electrodes  to  environmental  noise  sources  such  as  rainfall,  temperature, 
and  soil  property  variations.  Investigations  within  this  component  included 
studies  of  previous  publications  on  this  topic  and  laboratory  and  field  mea¬ 
surements.  Results  of  field  electrode  studies  at  Beaver  Lake  Dam,  Arkansas, 
are  included  in  Part  III  of  this  report. 

4.  As  discussed  in  the  early  publications  on  telluric  current  monitor¬ 
ing  referenced  below,  the  electrodes  are  the  most  critical  component  of  a  sys¬ 
tem  designed  to  monitor  SP  variations  over  long  periods  of  time.  Contact 
potentials  caused  by  interactions  between  the  electrodes  and  local  soil  con¬ 
ditions  may  be  larger  than  the  SP  voltages  generated  by  streaming  potentials 
related  to  seepage  flow,  and  time  variations  of  these  contact  potentials  may 
be  greater  than  those  related  to  seepage.  Therefore,  selection  of  the  proper 
electrode  type  and  installation  procedures  is  crucial  to  the  successful  per¬ 
formance  of  a  seepage  monitoring  array. 

5.  Initial  plans  had  been  to  study  a  wide  variety  of  polarizing  and 
nonpolarizing  electrode  types  (e.g.,  copper-copper  sulfate,  silver-silver 
chloride,  cadmium-cadmium  chloride,  lead-lead  chloride,  etc.)  as  well  as 
several  types  of  metallic  electrodes.  Also,  it  was  planned  to  study  nonpolar¬ 
izing  electrodes  in  which  the  liquid  electrolyte  and  porous  junction  are 
replaced  with  a  block  of  plaster  of  pairs  mixed  with  electrolyte  solution. 
However,  results  of  initial  studies,  time  constraints,  and  the  desirability  of 
focusing  on  previous  usage  led  to  a  change  in  emphasis  in  which  laboratory  and 
field  investigations  were  confined  to  three  electrode  types:  copper-copper 
sulfate  (for  brevity,  these  are  subsequently  referred  to  as  copper  sulfate) , 
copper-clad  steel  (CCS)  stakes,  and  metallic  lead.  Copper  sulfate  was 
selected  in  preference  to  other  nonpolarizing  electrode  types  because  of  com¬ 
mercial  availability  and  extensive  previous  usage  for  field  SP  measurements. 
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Previous  Investigations 


6.  Installed  pairs  or  arrays  of  electrodes  have  been  used  for  long¬ 
term  monitoring  of  natural  earth  currents  (tellurics)  related  to  magneto- 
spheric  activity,  SP  variations  that  could  be  earthquake  precursors,  and  of  SP 
signals  caused  by  possible  seepage-related  ground-water  flow  in  the  vicinity 
of  dams,  dikes,  and  other  containment  structures. 

7.  Van  Nostrand  and  Cook  (1966)  give  an  excellent  summary  of  the  work 
of  early  investigators  such  as  Matteucci,  Mauchley,  Gish,  and  Rooney  in  mea¬ 
suring  telluric  currents  using  permanently  installed  arrays  of  lead  elec¬ 
trodes.  These  references  are  included  in  the  bibliography  of  Part  V. 

8.  Lead  electrode  pairs,  carefully  prepared  and  installed,  were  found 
to  be  superior  to  nonpolarizing  electrodes  for  this  type  of  measurement 
(Mauchley  1918),  which  is  similar  in  concept  to  the  type  of  long-term  monitor¬ 
ing  desired  for  seepage  detection.  However,  as  discussed  later,  factors  other 
than  those  important  for  telluric  monitoring  must  be  considered  in  designing  a 
practical  seepage  monitoring  array. 

9.  Electrode  dipole  pairs  and  arrays  intended  to  measure  possible  SP 
precursors  to  earthquakes  have  been  installed  in  the  United  States,  the  Soviet 
Union,  and  China.  An  installation  of  this  type  is  described  in  Corwin  and 
Morrison  (1977).  Pairs  of  steel  plates  were  used  as  current  transmitters,  and 
lead  and  copper  sulfate  electrodes  were  used  as  potential  receivers  in  an 
experiment  to  determine  whether  electrical  resistivity  variations  preceded 
earthquakes  in  the  study  area  (near  Hollister,  California).  Before  and  after 
the  resistivity  measurements,  SP  was  continuously  monitored  across  most  of  the 
transmitter  and  receiver  electrode  pairs.  This  provided  a  nearly  continuous 
record  of  SP  variations  across  steel,  copper  sulfate,  and  lead  electrodes  for 
a  period  of  several  years.  Similar  records  probably  are  available  from 
investigators  involved  in  such  studies  in  other  locations,  but  little  of  the 
data  has  been  published. 

10.  Electrode  arrays  designed  to  detect  and  monitor  SP  variations 
related  to  seepage  flow  have  been  installed  at  a  number  of  locations  by  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES).  Descriptions  of  some  of 
these  installations  are  given  by  Llopis  (1987),  Llopis  and  Butler  (1988) 
Koester,  Butler,  Cooper  and  Llopis  (1984),  Cooper,  Koester,  and  Franklin 
(1982),  and  US  Army  Corps  of  Engineers,  Little  Rock  District  (1987).  All  of 
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these  arrays  utilized  copper-clad  steel  (CCS)  stakes  as  measuring  electrodes. 
CCS  stakes  also  were  used  in  a  monitoring  array  to  detect  sinkhole  drainage 
patterns  in  Virginia  (Erchul  1986) .  An  array  designed  to  monitor  seepage 
through  a  dike,  installed  by  the  Tennessee  Valley  Authority  (TVA)  at  Colbert, 
Alabama,  used  both  metal  stakes  and  copper  sulfate  electrodes  (R.  Hopkins*). 

11.  An  experiment  of  particular  interest  to  this  study  is  presently 
being  conducted  by  the  University  of  California  at  a  field  site  in  Richmond, 
California  (Wilt  et  al.  1987).  The  purpose  of  the  experiment  is  to  determine 
variations  of  resistivity  and  SP  in  response  to  downhole  saltwater  injection. 
Instrumentation  for  this  experiment  consists  of  arrays  of  surface  and  downhole 
transmitter  and  receiver  electrodes  along  with  a  computerized  data  acquisition 
and  control  system.  This  installation  is  very  similar  in  concept  and  instru¬ 
mentation  to  one  which  could  be  used  for  long-term  dam  seepage  monitoring. 

12.  The  receiver  electrode  array  at  the  Richmond  field  site  consists  of 
about  80  gelled  copper  sulfate  electrodes,  6  of  which  are  installed  in  bore¬ 
holes  and  the  remainder  in  a  surface  array  buried  at  depths  of  12-  to  18-in.** 
These  copper  sulfate  electrodes,  which  are  of  the  same  type  as  those  described 
below,  have  been  in  place  for  about  1  year  and  have  not  shown  any  signs  of 
deterioration. 

13.  Results  of  these  previous  studies  indicate  the  following  important 
points  regarding  selection  of  electrodes  for  long-term  seepage  monitoring 
arrays : 

a.  Lead  electrodes  for  telluric  monitoring  arrays  consist  of  large 
(several  feet)  grids  of  chemically  pure  lead  wire  buried  at 
depths  of  4-8  ft  (Rooney  2932;  1937).  These  preparation  and 
burial  requirements  increase  the  cost  of  lead  electrode  instal¬ 
lations  to  the  point  where  they  would  not  be  competitive  with 
copper  sulfate  or  CCS  stake  installations.  Even  with  such 
elaborate  installations,  yearly  drift  of  lead  electrodes  may  be 
several  to  several  tens  of  millivolts  (mV)  per  year.  As  indi¬ 
cated  by  the  results  of  the  field  experiment  described  later, 
electrodes  fabricated  of  industrial  (rather  than  laboratory) 
purity  lead  and  installed  at  relatively  shallow  depth  may  not 
give  stable  readings  for  long-term  measurements. 

b.  Copper  sulfate  (and  other  similar  types  of  nonpolarizing  elec¬ 
trodes  consisting  of  a  metal  rod  inserted  in  an  electrolyte 
solution  consisting  of  a  salt  of  the  metal  and  connected  to 


*  Personal  Communication,  1987,  R.  Hopkins,  TVA. 

**  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  6. 
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the  soil  through  an  inert  porous  membrane)  gives  very  stable 
potentials  even  in  the  presence  of  environmental  variations 
such  as  rainfall  and  temperature  changes.  Drawbacks  of  such 
electrodes  include  relatively  high  initial  cost,  care  needed  in 
installation,  and  the  possibility  of  electrolyte  leakage  or 
freezing  and  subsequent  failure  of  the  electrode.  Rooney 
(1937)  states 

It  is  sometimes  possible,  using  certain 
reversible  [nonpolarizing]  electrodes  - 
metallic  copper  in  copper-sulphate  solution, 
for  instance  -  to  reduce  the  absolute  values 
of  potential  recorded  between  two  given 
points  in  the  ground,  and  it  should  be  possi¬ 
ble,  theoretically,  to  keep  such  electrodes 
more  constant  than  the  simple  metallic  ones. 

However,  it  is  found  in  practice  that  the 
advantages  of  reversible  electrodes  for 
earth-current  work  are  not  readily  realized. 

In  order  to  make  contact  with  the  ground 
there  must  inevitably  be  a  slow  seepage  of 
the  electrode-solution  into  the  ground  about 
the  electrode.  This  results  in  a  constantly 
changing  electrode-environment  and  the  con¬ 
stancy  of  such  an  electrode  turns  out  to  be 
pretty  much  a  myth." 


c.  CCS  electrodes  offer  relatively  low  initial  cost,  ease  of 

installation,  and  low  maintenance  (assuming  that  the  electrodes 
are  not  damaged  by  on-site  activity) .  The  data  published  by 
WES  indicate  that  CCS  electrodes  may  in  some  way  amplify  the  SP 
signals  generated  by  seepage-related  streaming  potentials,  and 
that  geologic  noise  levels  (point-to-point  contact  potential 
variations  caused  by  local  soil  conditions)  are  considerably 
greater  than  those  for  copper  sulfate  and  other  nonpolarizing 
electrodes . 

14.  The  most  important  questions  raised  by  the  results  of  these  pre¬ 
vious  studies  are: 

a.  Are  the  careful  installation  procedures  necessary  for  lead 
electrodes  cost-effective  for  seepage  monitoring  arrays? 

b.  Is  the  larger  amplitude  of  the  variations  observed  on  CCS 
electrode  arrays  a  true  amplification  of  seepage-related  sig¬ 
nals  or  an  apparent  effect  due  to  other  causes? 

c.  If  the  CCS  electrodes  provide  true  signal  amplifications,  what 
is  its  cause? 

d.  How  do  signal-to-noise  ratios  compare  for  CCS,  copper  sulfate, 
and  lead  electrodes? 
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e.  Considering  the  above  factors,  which  electrode  type  provides 
the  best  combination  of  accuracy  and  cost-effectiveness  for 
seepage  monitoring  arrays. 

The  studies  described  below  were  designed  to  help  answer  these  questions. 


Long-Term  Monitoring  Study 


15.  This  section  describes  the  results  of  a  study  in  which  three  elec¬ 
trode  pairs  (copper  sulfate,  metallic  lead,  and  CCS  stakes)  were  installed  in 
a  location  protected  from  disturbance  in  a  suburban  residential  area  and 
monitored  for  a  period  of  about  1  year.  The  purposes  of  the  study  were  to 
determine  the  time  needed  for  the  electrode  pairs  to  stabilize,  the  response 
of  the  different  electrode  types  to  environmental  variations  such  as  rainfall 
and  temperature  changes,  and  the  long-term  stability  of  the  different  elec¬ 
trode  types. 

Description  of  installation 

16.  Plate  1  shows  the  arrangement  of  the  electrodes.  Electrode  types 
were  as  follows: 

a.  Copper  sulfate:  Tinker  and  Rasor  Model  6B,  1-in.  diam  by  6-in. 
length  plastic  body  with  1-in.  diam  porous  ceramic  junction, 
initially  filled  with  saturated  copper  sulfate  solution  (as 
discussed  later,  the  liquid  electrolyte  was  replaced  with 
gelled  electrolyte  after  about  6  months  of  operation) .  The 
electrodes  were  fitted  with  submersible  waterproof  adaptors 
(Model  W-7)  to  allow  full  burial.  The  relatively  small  porous 
junction  area  and  O-ring  seals  of  these  electrodes  should  have 
rendered  them  as  leakproof  as  possible  for  commercially  avail¬ 
able  off-the-shelf  copper  sulfate  electrodes.  The  electrodes 
were  installed  vertically,  with  the  tips  at  a  depth  of  about 

1  ft,  and  completely  covered  with  soil. 

b.  Copper-clad  steel  (CCS)  stakes:  These  electrodes  were  1/2-in. 
diam  by  2-ft-long  copper-clad  steel  grounding  rods.  They  were 
driven  to  a  depth  of  1  ft  into  the  soil,  and  stranded  cooper 
lead-in  wires  were  clamped  to  the  tops  of  the  rods  with  stain¬ 
less  steel  hose  clamps.  The  CCS  stakes  were  left  partially 
exposed  (and  not  completely  buried)  to  duplicate  previous  Corps 
of  Engineers  field  procedure.  The  clamps  were  coated  with 
insulating  compound  (Scotchkote)  to  minimize  galvanic  corrosion 
between  the  stainless  steel  clamps  and  the  copper  rod  cladding. 

c.  Metallic  lead  electrodes:  These  electrodes  were  cut  to  a  size 
of  2  in.  by  4.5  in.  from  99  percent  pure  lead  sheet  of 
0.042-in.  thickness  (2.5-lb/sq  ft  roof  flashing).  Stranded 
copper  lead-in  wire  was  soldered  to  the  top  of  the  long  side  of 
each  electrode  and  the  solder  joint  coated  with  Scotchkote. 
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The  electrodes  were  buffed  with  steel  wool,  washed  with  soap 
and  water,  and  rinsed  before  installation.  They  were  buried 
with  the  4.5-in.  side  vertical,  with  their  center  points  at  a 
depth  of  about  8  in. 

The  chronology  of  the  installation  was  as  follows: 

a.  31  January  1986:  Installation  of  copper  sulfate  and  CCS 
electrodes . 

b.  4  March  1987:  Installation  of  lead  electrodes. 

c.  22  August  1987:  As  the  data  were  indicating  problems  with  the 
copper  sulfate  electrodes,  they  were  removed  for  inspection  and 
found  to  be  dry  (i.e.,  all  the  electrolyte  solution  had  leaked 
out) .  The  electrolyte  was  replaced  with  a  gel  consisting  of 
two  Knox  gelatins  (by  weight)  added  to  heated  saturated  copper 
sulfate  solution  and  allowed  to  set.  This  conductive  gel  had  a 
consistency  similar  to  that  of  Jell-0,  which  it  was  hoped  would 
eliminate  or  greatly  retard  electrolyte  leakage. 

d.  15  January  1988:  The  copper  sulfate  electrodes  were  removed 
for  inspection  and  found  to  be  in  good  condition,  with  no 
evidence  of  electrolyte  leakage  or  junction  deterioration. 

17.  The  potentials  across  these  three  electrode  pairs  were  monitored 
using  a  combination  of  strip  chart  recorders  and  digital  multimeters.  A  rain 
gage  installed  as  shown  in  Plate  1  was  used  to  monitor  precipitation  beginning 
in  the  fall  of  1987.  Due  to  occasional  equipment  problems  and  absence  of 
personnel  from  the  site  the  records  are  not  fully  continuous.  However, 
sufficient  data  were  obtained  to  meet  the  objectives  of  the  experiment. 

Results 

18.  A  strip  chart  record  of  the  data  from  the  copper  sulfate  and  CCS 
electrodes  for  a  period  of  about  6.5  min  is  shown  in  Plate  2.  The  record 
shows  continuous  variations  with  amplitudes  ranging  from  a  few  tenths  of  a 
millivolt  (mV)  to  a  few  mV,  with  periods  of  a  few  seconds  to  several  tens  of 
seconds.  These  variations  are  caused  by  stray  currents  from  the  electrically 
operated  Bay  Area  Rapid  Transit  System  (BART),  located  about  1  mile  east  of 
the  monitoring  electrode  (the  daily  records  discussed  below  show  the  elimina¬ 
tion  of  these  variations  when  BART  is  not  operating). 

19.  The  important  conclusion  reached  from  this  record  is  that  copper 
sulfate  and  CCS  electrodes  show  essentially  identical  responses  to  electric 
fields  generated  by  outside  sources  (the  response  of  the  lead  electrodes  to 
the  BART-generated  signals  also  was  identical  to  that  of  the  copper  sulfate 
and  CCS  electrodes) .  Thus  if  SP  signals  generated  by  seepage  flow  are  of 


similar  nature  to  those  generated  by  stray  currents,  the  SP  voltage  measured 
by  copper  sulfate,  CCS,  or  lead  electrodes  should  be  identical. 

20.  A  sample  record  of  copper  sulfate  and  lead  electrode  potentials  for 
a  period  of  about  3  days  is  shown  in  Plate  3.  Note  the  strong  reduction  of 
the  short-period  noise  level  between  about  11  p.m.  and  3  a.m.  each  day,  when 
the  BART  system  is  not  operating.  The  intrinsic  short-period  noise  level  for 
both  electrode  pairs  is  of  the  order  of  about  1  mV  during  these  intervals. 
Similar  results  were  seen  for  the  CCS  electrodes. 

21.  Intermittent  rain  fell  during  the  period  shown  in  Plate  3,  which 
illustrates  the  response  of  lead  and  copper  sulfate  electrodes.  The 
responses  of  the  copper  sulfate  and  lead  electrodes  during  this  rainfall  were 
representative  of  that  seen  during  most  of  the  experiment:  little  or  no 
variation  for  the  copper  sulfate  electrodes  and  the  intermittent  appearance  of 
DC  offsets  on  the  lead  electrodes.  The  copper  sulfate  electrodes  responded 
only  to  very  heavy  rainfall,  with  DC  levels  returning  to  previous  values  once 
the  rain  ended.  The  response  of  the  CCS  electrodes  is  similar  to  that  of  the 
lead  electrodes.  Note  that  CCS  and  copper  sulfalte  electrode  response  to 
rainfall  is  also  discussed  in  Part  III  of  this  report. 

22.  Plate  4  shows  the  effects  of  temperature  variations  on  the  copper 
sulfate  and  CCS  electrode  pairs.  Potential  for  both  the  copper  sulfate  and 
CCS  electrodes  was  relatively  stable  on  6  November,  which  was  cool  and  over¬ 
cast.  On  7  November,  a  warm  and  sunny  day,  the  CCS  electrodes  showed  a  posi¬ 
tive  excursion  of  about  15  mV  corresponding  to  the  heating  of  the  electrodes 
by  direct  sunlight.  The  buried  copper  sulfate  (and  lead)  electrodes  showed  no 
such  diurnal  response.  Although  of  variable  magnitude,  this  temperature 
response  of  the  CCS  electrodes  was  seen  on  most  sunny  days.  Diurnal  varia¬ 
tions  of  this  type  would  not  be  critical  for  long-term  monitoring  over  periods 
of  several  months  or  more,  but  seasonal  DC  offsets  caused  by  electrode  temper¬ 
ature  variations  could  be  of  concern.  Complete  burial  of  the  CCS  electrodes 
would  probably  reduce  the  SP  "noise”  due  to  temperature  variation. 

23.  Summarized  data  for  the  entire  term  of  the  monitoring  experiment 
are  shown  in  Plate  5.  The  plotted  values  are  estimates  of  the  daily  DC  poten¬ 
tial  levels  that  a  long-term  seepage  monitoring  network  would  be  designed  to 
record.  Gaps  in  the  records  are  periods  when  the  recording  instruments  mal¬ 
functioned  and/or  no  personnel  were  available  to  record  the  data. 

24.  Several  conclusions  may  be  reached  from  the  data  shown  in  Plate  5: 
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The  copper  sulfate  electrodes  reached  a  stable  potential  value 
almost  immediately  (within  a  few  minutes)  after  installation. 
The  CCS  and  lead  electrodes  showed  large  potential  excursions 
(several  hundred  mV)  over  a  period  of  1-  to  2-months  after 
installation  before  reaching  more  stable  values. 

During  the  dry  period  between  about  mid-May  and  mid-October 
1987,  the  lead  electrodes  showed  little  daily  variation  and  a 
long-term  negative  drift  of  about  25  mV.  Their  DC  value 
appeared  to  be  relatively  constant  at  about  +75  mV  for  the 
final  3  months  of  this  period.  The  CCS  electrodes  showed 
somewhat  greater  daily  variation  than  the  lead  electrodes 
during  this  period,  with  occasional  large  (50-70  mV)  positive 
excursions  lasting  for  a  few  days.  The  DC  level  of  the  copper 
sulfate  electrodes  was  stable  within  a  few  mV  of  its  initial 
value  from  the  time  of  its  installation  through  mid-August, 
when  the  electrolyte  solution  had  leaked  out.  After  being 
reinstalled  with  gelled  electrolyte,  the  DC  level  of  these 
electrodes  returned  to  its  previous  value  and  stayed  constant 
through  mid-October. 

Seasonal  rainfall  began  in  mid-October  1987  and  continued 
through  the  end  of  the  monitoring  period  in  mid-January  1988. 
This  rainfall  appeared  to  produce  the  following  effects: 

(1)  The  DC  potential  of  the  lead  electrodes  dropped  from  about 
+75  mV  to  about  +10  mV  between  mid-October  and  the  end  of 
October,  then  rose  to  about  +30  MV  between  the  end  of 
November  1987  and  mid-January  1988.  Day-to-day  variations 
were  somewhat  greater  than  those  seen  during  the  preceding 
dry  period. 

(2)  The  DC  potential  of  the  CCS  electrode  pair  dropped  from  an 
average  value  of  about  +20  mV  during  the  dry  period  to 
about  -5  mV  following  the  onset  of  the  rainy  season. 
Day-to-day  variations  appeared  to  be  somewhat  lower  during 
the  rainy  period,  with  no  consistent  response  to  individ¬ 
ual  periods  of  rain.  The  occasional  large  positive  excur¬ 
sions  seen  during  the  dry  period  were  less  frequent  during 
the  rainy  season. 

(3)  The  DC  value  of  the  potential  between  the  copper  sulfate 
electrode  pair  remained  within  about  +/- 2  mV  of  its  pre¬ 
rainfall  mean  of  about  -3  mV.  However,  individual  periods 
of  heavy  rain  often  were  followed  by  negative  excursions 
of  about  -10  mV,  lasting  1-2  days.  As  the  electrode  pair 
was  installed  on  a  hillside  with  the  negative  electrode  at 
about  4  ft  lower  elevation  than  the  positive  electrode, 
these  excursions  have  the  correct  polarity  for  a  streaming 
potential  generated  by  downslope  movement  of  near-surface 
ground  water.  If  this  is  the  case,  it  is  not  apparent  why 
similar  excursions  were  not  seen  on  the  CCS  and  lead 
electrode  pairs. 
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Discussion 


25.  From  the  results  described  above,  the  most  significant  conclusions 
of  the  long-term  monitoring  study  are  as  follows: 

a.  Stabilization  period:  The  copper  sulfate  electrode  pair 
reached  a  stable  DC  potential  value  within  a  few  minutes  of 
installation  and  remained  at  that  value  for  the  duration  of  the 
study.  These  results  are  similar  to  those  seen  for  previous 
copper  sulfate  monitoring  arrays  (Corwin  and  Morrison  1977) . 

The  CCS  electrodes  showed  large  potential  excursions  for  a 
period  of  about  3  months  following  installation,  after  which 
they  remained  in  a  band  of  about  30  mV  with  the  exception  of 
occasional  +50-  to  +70-mV  excursions  of  a  few  days  duration 
during  the  dry  season.  Additional  monitoring  would  be  neces¬ 
sary  to  determine  whether  the  CCS  electrodes  have  actually 
attained  a  stable  DC  value.  The  lead  electrodes  also  required 
a  period  of  about  3  months  after  installation  to  settle  to  a 
relatively  constant  value,  but  drift  continued  even  through  the 
4-month  dry  season.  Continuing  changes  of  the  DC  level  of 
these  electrodes  during  the  rainy  sea^  n  indicate  that  a  stable 
DC  level  was  not  reached  within  the  11  months  following  their 
installation. 

b.  Response  to  external  electric  fields:  The  experimental  results 
indicate  that  all  three  electrode  types  have  identical  response 
to  external  electric  fields  in  the  earth  that  have  periods  of  a 
few  seconds  to  a  few  minutes. 

c.  Day-to-day  variations:  Disregarding  the  effects  of  rainfall, 
the  copper  sulfate  electrodes  showed  the  lowest  level  of  day- 
to-day  variability,  with  most  of  the  variation  due  to  the  dif¬ 
ficulty  of  estimating  DC  levels  in  the  presence  of  several 
millivolts  of  background  noise  caused  by  the  BART  system. 
Day-to-day  variation  of  the  lead  electrodes  during  the  dry  sea¬ 
son  was  almost  as  low  as  that  of  the  copper  sulfate  electrodes, 
probably  because  their  burial  below  the  surface  isolated  them 
from  direct  heating  by  sunlight. 

d.  The  CCS  electrodes  showed  the  highest  day-to-day  variability, 
in  part  due  to  direct  heating  effects.  It  also  appears  that, 
especially  during  dry  weather,  the  CCS  electrodes  were  subject 
to  significant  potential  excursions,  probably  caused  by  the 
buildup  and  subsequent  decay  of  periodic  corrosion  reactions, 

e.  Reliability:  The  CCS  and  lead  electrodes  needed  no  maintenance 
after  installation.  The  copper  sulfate  electrodes  with  liquid 
electrolyte,  even  though  well  sealed,  leaked  out  all  their 
electrolyte  within  6-1/2  months  of  installation.  Inspection  of 
the  gelled  copper  sulfate  electrodes  after  about  16  months  of 
operation  showed  no  electrolyte  loss.  Although  it  should  be 
kept  in  mind  that  the  soil  was  saturated  during  about  half  of 
this  period,  which  would  tend  to  retard  leakage,  the  gelled 
electrodes  appear  to  be  capable  of  long-term,  unmaintained 
burial. 
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26.  Summarizing  the  findings  of  previous  investigations  and  the  field 
experiment  described  above,  copper  sulfate  electrodes  appear  to  stabilize  much 
more  rapidly  than  do  CCS  or  lead  electrodes  and  to  be  capable  of  maintaining  a 
stable  DC  value  indefinitely.  Lead  electrodes  made  of  commercial  grade  lead 
and  buried  at  depths  practical  for  multiple  electrode  arrays  appear  to  perform 
poorly,  showing  DC  offsets  in  response  to  rainfall  even  after  almost  1  year  of 
burial.  Although  lead  is  the  electrode  material  of  choice  for  long-term 
telluric  monitoring  arrays  consisting  of  a  few  carefully  prepared  electrodes 
constructed  of  expensive  chemically  pure  material  and  buried  at  depths  of  4- 
to  8-ft,  it  does  not  appear  to  be  practical  for  arrays  of  tens  or  hundreds  of 
electrodes  constructed  of  affordable  commercial  grade  material  to  be  buried  at 
depths  of  2  ft  or  less. 

27.  CCS  electrodes  are  considerably  noisier  and  less  stable  than  copper 
sulfate.  If  the  response  of  CCS  electrodes  to  SP  signals  generated  by  seepage 
flow  is  the  same  as  that  of  copper  sulfate,  the  choice  between  them  would  be 
based  on  the  superior  signal-to-noise  (S/N)  ratio  of  copper  sulfate  versus  the 
lower  initial  cost  and  possibly  lower  maintenance  costs  of  CCS.  However, 
there  is  evidence  from  WES  dam  seepage  studies  using  CCS  electrodes  that  SP 
signal  levels  as  well  as  noise  levels  are  greater  for  CCS  electrodes  as  com¬ 
pared  with  copper  sulfate,  giving  comparable  S/N  ratios  for  the  two  electrode 
types.  In  this  context,  "noise"  is  any  electrical  potential  variation  that  is 
not  related  to  the  phenomena  of  interest,  i.e.,  flow  of  water  (seepage)  in  the 
subsurface.  A  number  of  the  laboratory  tests  described  in  the  following  sec¬ 
tion  were  conducted  to  attempt  to  determine  whether  CCS  electrodes  give 
greater  SP  signal  levels  for  a  given  seepage  flow  than  do  copper  sulfate 
electrodes . 


Laboratory  Measurements 

28.  The  laboratory  measurements  described  in  this  section  were  con¬ 
ducted  to  study  the  response  of  CCS  and  copper  sulfate  electrodes  to  environ¬ 
mental  variations  such  as  temperature  and  soil  moisture  content  and  to  examine 
the  SP  response  of  these  electrode  types  to  the  flow  of  water  through  a  porous 
medium. 
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Freezing  of  copper  sulfate  electrodes 

29.  Freezing  of  the  liquid  or  gelled  electrolyte  of  a  nonpolarizing 
electrode  could  result  in  damage  to  the  electrode  and/or  degradation  of 
electrode  performance.  To  examine  the  results  of  electrolyte  freezing,  one  of 
the  gelled  copper  sulfate  electrodes  described  previously  was  held  overnight 
in  a  freezer  at  8°  F. 

30.  Inspection  of  the  electrode  showed  that  the  electrolyte  was  frozen 
solid,  but  that  there  was  no  apparent  physical  damage  to  the  electrode. 

Values  of  potential  and  contact  resistance  between  the  frozen  electrode  and  an 
identical  nonfrozen  electrode  in  a  bath  of  copper  sulfate  solution  at  various 
stages  in  the  experiment  are  tabulated  as  follows.  The  negative  lead  of  the 
digital  multimeter  used  for  the  measurements  was  connected  to  the  unfrozen 


electrode. 

Contact 

Potential  Resistance 

_ Condition _  (mV)  (kohm) 

Before  freezing  +1.7  1.71 

Positive  electrode  partly  frozen  -9.8  86 

Positive  electrode  fully  frozen  -16.1  770 

Positive  electrode  fully  thawed  +1.1  1.50 


After  conclusion  of  the  freezing  experiment,  the  electrodes  were  reinstalled 
in  the  earth.  Their  performance  after  reinstallation  appeared  to  be  identical 
to  that  before  freezing. 

31.  These  results  indicate  that  freezing  and  subsequent  thawing  of  this 
type  of  electrode  do  not  appear  to  damage  the  electrode  or  affect  its  perfor¬ 
mance.  However,  the  field  measuring  system  must  be  designed  with  high  enough 
input  impedance  to  handle  the  increased  contact  resistance  for  frozen  elec¬ 
trodes.  Alternatively,  the  gel  could  include  a  percentage  of  antifreeze  solu¬ 
tion  (such  as  glycol)  to  retard  freezing  at  the  cost  of  somewhat  greater 
unfrozen  resistance. 

32.  The  18-mV-negative  potential  change  for  the  60°  F  temperature  drop 
is  typical  for  copper  sulfate  electrodes.  In  a  field  installation,  maximum 
temperature  differences  would  be  considerably  less,  and  maximum  temperature 
effects  probably  would  be  limited  to  a  few  millivolts.  Temperature  effects 
are  discussed  in  more  detail  in  later  sections. 
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Streaming  potential  measurements 

33.  The  streaming  potential  cell  shown  in  Plate  6  was  designed  to 
examine  the  SP  response  of  CCS  and  copper  sulfate  electrodes  to  water  flow. 

The  cell  was  made  from  thin-walled  plastic  tubing,  and  all  fittings  and 
connections  also  were  of  plastic  to  eliminate  possible  corrosion  potentials 
from  metal  components.  The  cell  was  filled  with  medium-grained  sand  of  mixed 
mineralogy  and  saturated  with  tap  water  of  about  300-ohm-ft  resistivity. 

34.  The  CCS  electrodes  were  6-in.  lengths  of  the  same  1/2-in.  diameter 
copper-clad  steel  grounding  rods  described  previously,  and  the  copper  sulfate 
electrodes  were  Tinker  &  Rasor  Model  2A  "pocket  cell"  electrodes  (3/8-in.  out¬ 
side  diameter)  filled  with  gelled  copper  sulfate  solution  as  described  above. 
The  CCS  electrodes  were  allowed  to  equilibrate  in  the  saturated  sand  for  about 
1  week  prior  to  beginning  the  experiment. 

35.  Flow  rates  for  the  test  runs  were  kept  low  (0.008  to  0.396  gpm)  in 
an  effort  to  maintain  laminar  flow.  The  relationship  between  flow  rate  and 
pressure  drop  across  the  electrodes  was  determined  by  replacing  the  electrodes 
with  a  pair  of  manometer  tubes  and  measuring  the  water  height  difference  in 
the  tubes  at  a  known  flow  rate. 

36.  The  first  test  run  was  made  using  the  CCS  electrodes  in  the  config¬ 
uration  shown  in  Plate  6,  with  the  electrodes  extending  to  within  1/4  in.  of 
the  opposite  wall  of  the  tube.  For  the  second  set  of  runs,  both  CCS  and 
copper  sulfate  electrodes  were  installed  with  their  lower  (inner)  ends  approx¬ 
imately  flush  with  the  bottoms  of  the  plastic  fittings.  For  the  final  set  of 
test  runs,  both  electrode  types  were  completely  isolated  from  the  water  flow 
by  means  of  a  small  plug  of  clay  inside  the  plastic  fitting.  The  purpose  of 
these  successive  removals  of  the  electrodes  from  the  flow  was  to  attempt  to 
separate  effects  caused  by  interaction  between  the  moving  water  and  the  elec¬ 
trode  surface  (so-called  "flow  potentials")  from  those  due  to  streaming  poten¬ 
tials  generated  by  flow  through  the  porous  medium. 

37.  Results  are  tabulated  below.  Each  listed  value  of  measured  poten¬ 
tial  is  the  average  of  several  runs  at  different  flow  rates.  Because  the  zero 
values  of  the  measured  potentials  for  the  CCS  electrodes  were  not  stable, 
reproducibility  of  the  listed  values  probably  is  about  ±25  percent.  Copper 
sulfate  values  probably  are  reproducible  within  about  ±10  percent. 
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Measured  potential 

_ (mV/atm) _ 

Copper 

Condition  CCS  Sulfate 


Electrodes  in  flow  4,670 
Electrodes  at  bottom  of  fitting  1,260 
Electrodes  isolated  by  clay  plug  510 


420 

380 


38.  Potentials  measured  by  the  two  electrode  types  isolated  from  the 
flow  by  the  clay  plug  were  not  significantly  different.  The  magnitude  of 
about  400  mV/atm  is  similar  to  that  measured  by  other  investigators  for  simi¬ 
lar  test  configurations  (e.g.,  Ogilvy,  Ayed,  and  Bogoslovsky  1969;  Tuman 
1963).  When  the  electrodes  were  partially  isolated  from  the  flow,  the  poten¬ 
tial  change  for  the  copper  sulfate  electrodes  was  not  significant,  but  the  CCS 
electrode  reading  increased  by  a  factor  of  about  2.5.  When  the  CCS  electrodes 
were  arranged  to  span  nearly  the  full  width  of  the  tube,  exposing  them  to  the 
maximum  flow,  their  potential  increased  to  a  value  about  nine  times  greater 
than  that  for  the  isolated  condition. 

39.  The  value  of  about  4,670  mV/atm  for  the  CCS  electrodes  is  much 
greater  than  that  expected  for  a  true  streaming  potential.  This  discrepancy, 
along  with  the  large  potential  differences  between  isolated  and  nonisolated 
flow,  imply  that  the  CCS  electrodes  are  subject  to  significant  "flow  poten¬ 
tials"  caused  by  the  movement  of  water  along  the  electrode  surface.  These 
flow  potentials  for  metal  electrodes  have  been  observed  in  previous  studies 
(Ogilvy,  Ayed,  and  Bogoslovsky  1969).  Conversely,  nonpolarizing  electrodes 
appear  to  be  free  from  flow  potential  effects  (Ogilvy,  Ayed,  and  Bogoslovsky 
1969;  Corwin  and  Conti  1973). 

40.  With  the  electrodes  isolated  from  the  flow,  as  would  be  the  case 
for  a  typical  field  SP  measurement,  streaming  potential  response  appears  to  be 
about  the  same  for  CCS  and  copper  sulfate  electrodes.  This  raises  the  ques¬ 
tion  of  the  origin  of  the  large-amplitude  anomalies  (greater  than  100  mV) 
observed  on  CCS  electrode  arrays  that  appear  to  be  related  to  subsurface  seep¬ 
age  flow  but  are  not  seen  on  measurements  taken  at  the  same  stations  with 
copper  sulfate  electrodes  (see  Part  III). 

41.  One  possibility  is  that  some  of  these  anomalies  may  be  generated  by 
flow  potentials  caused  by  movement  of  water  through  surface  soils,  with  this 
movement  related  in  turn  to  the  presence  of  faults  or  fractures  that  extend 
from  the  surface  to  the  zone  of  seepage.  As  discussed  below,  soil  property 


20 


variations  associated  with  surface  expression  of  subsurface  features  also 
could  generate  anomalies  on  CCS  electrodes  that  are  indicative  of  seepage 
zones . 

Temperature  effects 

42.  Diurnal  SP  variations  on  a  CCS  electrode  pair  that  appeared  to  be 
related  to  electrode  temperature  differences  were  described  above.  While  such 
variations  probably  would  tend  to  average  out  during  the  course  of  long-term 
monitoring,  they  could  generate  spurious  anomalies  in  day-to-day  measurements. 
Therefore,  it  is  useful  to  have  some  idea  of  the  magnitude  of  SP  changes  gen¬ 
erated  by  electrode  temperature  differences. 

43.  The  temperature  response  of  nonpolarizing  electrodes  can  be  derived 
analytically  (Ives  and  Janz  1961;  Maclnnes  1961),  and  has  been  measured  for  a 
number  of  different  electrode  types.  Such  measurements  have  been  performed 
for  copper  sulfate  electrodes  by  Semenov  (1974),  Ewing  (1939),  and  Morrison 

et  al.  (1979a).  For  a  variety  of  experimental  conditions,  the  temperature 
coefficient  of  copper  sulfate  electrodes  was  found  to  average  about  +0.5  mV 
per  degree  F  (i.e.,  heating  the  positive  electrode  caused  a  positive  SP 
change).  It  should  be  noted  that  temperatures  do  not  refer  to  the  soil  around 
the  electrode,  but  to  the  metallic  element  of  the  nonpolarizing  electrode. 
Thus,  there  is  a  time  lag  between  heating  of  the  soil  and  appearance  of  an  SP 
change  due  to  heating  of  the  electrode  element. 

44.  Because  no  such  measurements  appear  to  have  been  made  for  CCS  elec¬ 
trodes,  a  laboratory  study  was  performed  to  determine  their  temperature  coef¬ 
ficient.  The  electrode  pair  used  for  the  streaming  potential  measurement 
described  above  was  installed  in  a  plastic  container  filled  with  the  same 
saturated  sand  used  in  the  streaming  potential  cell.  The  electrodes  were 
immersed  to  a  depth  of  3  in.  in  the  sand  and  separated  by  7  in.  As  for  the 
streaming  potential  cell,  a  digital  multimeter  and  strip  chart  recorder  were 
connected  across  the  electrode  pair,  with  their  positive  terminals  to  the 
heated  electrode. 

45.  Thermometers  were  strapped  to  the  two  electrodes  to  measure  their 
temperature  while  one  of  the  electrodes  was  heated  with  a  hot  air  gun. 

Response  to  temperature  changes  was  almost  immediate,  and  averaged  about  +1  mV 
per  degree  F  for  a  series  of  temperature  differences  ranging  up  to  77°  F. 

46.  This  coefficient  is  about  twice  that  measured  for  copper  sulfate 
electrodes.  Because  the  maximum  diurnal  difference  in  surface  soil 
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temperature  between  any  two  electrodes  in  a  monitoring  array  may  be  of  the 
order  of  several  degrees,  soil  temperature  differences  would  be  expected  to 
produce  diurnal  SP  variations  of  only  a  few  millivolts  for  electrodes  not 
exposed  to  direct  sunlight,  as  would  be  the  case  for  buried  electrodes. 

47.  For  CCS  stake  electrodes  exposed  to  direct  sunlight,  however, 
diurnal  temperature  differences  between  exposed  and  shaded  stakes,  and  con¬ 
sequent  diurnal  SP  variations,  could  be  much  greater.  Examples  of  such  varia¬ 
tions  are  seen  in  Plate  4.  Their  magnitude  is  large  enough  that  they  could 
produce  significant  spurious  long-period  anomalies  by  aliasing  of  daily  read¬ 
ings  under  conditions  of  changing  solar  electrode  heating.  Therefore,  for 
arrays  of  metal  stakes  exposed  to  solar  heating,  care  must  be  taken  that  tem¬ 
perature  effects  are  accounted  for.  Results  of  the  long-term  monitoring 
experiment  described  above  indicate  that  temperature  effects  for  either  metal 
or  nonpolarizing  electrodes  buried  more  than  a  few  inches  beneath  the  surface 
are  small  and  generally  can  be  ignored. 

Soil  moisture  variations 

48.  A  series  of  laboratory  and  field  studies  performed  by  Morrison 

et  al.  (1978,  1979a,  1979b)  indicated  that  the  most  significant  soil  property 
affecting  field  SP  readings  taken  with  nonpolarizing  electrodes  was  the  mois¬ 
ture  content  of  the  soil  at  the  electrode  measurement  station.  Changes  in 
moisture  content  were  found  to  have  a  greater  effect  than  variations  of  soil 
type  or  of  soil  moisture  salinity,  pH,  or  Eh.  For  copper  sulfate  electrodes 
in  a  variety  of  soil  types,  the  average  response  to  soil  moisture  changes  was 
about  +1  mV  per  percent  increase  in  soil  moisture  content  (i.e.,  the  electrode 
in  the  wetter  soil  became  more  positive  with  respect  to  that  in  the  dryer 
soil)  . 

49.  Because  similar  measurements  do  not  appear  to  have  been  made  for 
CCS  electrodes,  a  laboratory  study  of  the  response  of  CCS  electrodes  to  soil 
moisture  variations  was  performed.  The  electrodes  used  were  those  described 
above  for  the  streaming  potential  and  temperature  measurements,  and  the  soil 
was  a  silty  clay  obtained  from  the  Beaver  Lake  Dam  site  described  in  Part  Ill 
(the  sample  is  representative  of  the  soil  in  which  most  of  the  CCf-  electrodes 
were  installed) . 

50.  Samples  of  nominal  (in  situ  moisture  content)  and  wet  soil  were 
placed  in  contact  in  a  small  plastic  container.  (Water  from  exit  6  at  Beaver 
Dam  was  added  to  the  nominal  soil  to  obtain  the  wet  sample.)  Moisture  content 
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for  the  two  samples  was  determined  by  weighing  the  wet  samples,  drying  them 
over  low  heat,  and  reweighing.  Moisture  contents  by  weight  for  the  wet  and 
dry  samples  were  about  34  and  30  percent,  respectively.  The  value  for  the  wet 
sample  may  be  a  few  percent  low  due  to  the  presence  of  some  gravel  in  the 
sample . 

51.  A  digital  multimeter  and  strip  chart  recorder  were  connected  to  the 
electrodes,  which  were  allowed  to  equilibrate  in  the  nominal  soil.  The  posi¬ 
tive  electrode  then  was  moved  to  the  wet  soil,  and  SP  values  allowed  to  equil¬ 
ibrate.  The  positive  electrode  then  was  moved  back  to  the  nominal  soil  and 
the  process  repeated  several  times.  For  comparison,  this  same  procedure  also 
was  carried  out  with  copper  sulfate  electrodes. 

52.  It  proved  impossible  to  obtain  a  reproducible  response  from  the  CCS 
electrodes  for  the  first  few  trials.  Zero  levels  changed  up  to  100  mV  each 
time  the  positive  electrode  was  reinstalled  in  the  nominal  soil,  and  the 
response  to  placement  in  the  wet  soil  drifted  for  hours  without  reaching  a 
steady  value.  Quasi-steady  values  ranged  from  within  a  few  mV  of  the  value 
with  both  electrodes  in  the  nominal  soil  to  +/-100  mV  of  this  value.  After 
several  trials,  however,  somewhat  more  reproducible  readings  could  be 
obtained,  for  which  the  quasi-steady  potential  appeared  to  be  about  50  mV  more 
negative  with  Che  positive  electrodes  in  the  nominal  soil. 

53.  In  contrast,  the  copper  sulfate  electrodes  exhibited  a  reproducible 
response  of  +6-mV  difference  between  the  wet  and  nominal  soils  (i.e.,  the 
electrode  in  the  wet  soil  was  positive  with  respect  to  that  in  the  dry  soil). 
This  value  is  reasonably  consistent  with  the  measured  change  in  soil  moisture 
content  (given  that,  as  discussed  above,  this  change  may  actually  be  somewhat 
greater  than  4  percent). 

54.  The  results  of  this  experiment  indicate  that  CCS  electrodes  are 
extremely  sensitive  to  surface  effects  caused  by  movement  of  the  electrode 
within  the  soil  (it  was  observed  that  even  slight  disturbance  of  one  of  the 
electrodes  often  produced  SP  changes  of  several  tens  of  mV  or  more).  These 
surface  effects  may  be  related  to  formation  and  removal  of  corrosion  reaction 
products.  The  results  also  indicate  that  steady-state  response  of  CCS  elec¬ 
trodes  to  soil  moisture  changes  may  be  of  the  order  of  about  -10  mV  per  per¬ 
cent  moisture  change. 

55.  In  a  clay-rich  soil  such  as  that  from  the  Beaver  Dam  site  moisture 
content  may  vary  by  as  much  as  about  10  to  15  percent  between  very  wet  and 


very  dry  stations.  Thus,  SP  variations  measured  by  CCS  electrodes  due  to 
changes  in  moisture  content  in  such  soils  may  be  as  much  as  100  or  150  mV. 
However,  it  does  not  appear  that  moisture  content  effects  alone  are  the  source 
of  station-to-station  SP  variations  of  a  few  hundred  mV  or  more. 

Conclusions  of  Electrode  Study 

56.  Of  the  three  electrode  types  studied,  electrodes  of  commercial- 
grade  metallic  lead  appear  to  be  the  least  suitable  for  long-term  monitoring 
of  SP  variations  related  to  seepage  flow.  Buried  at  a  depth  of  less  than 

1  ft,  lead  electrodes  do  not  appear  to  achieve  a  stable  potential  during  dry 
weather,  and  exhibit  unstable  DC  shifts  in  response  to  rainfall  even  after  a 
year  of  burial. 

57.  Copper  sulfate  electrodes  with  gelled  electrolyte  appear  to  be 
capable  of  surviving  at  least  a  few  years  of  burial  without  maintenance  or 
significant  deterioration  of  physical  properties  or  performance.  Response  of 
these  electrodes  to  changes  of  environmental  parameters  such  as  temperature  or 
soil  moisture  content  is  relatively  small,  consistent,  and  reversible.  These 
responses  are  of  the  order  of  a  few  mV,  compared  with  typical  streaming  poten¬ 
tial  anomalies  at  dam  sites  of  a  few  tens  of  mV. 

58.  Thus,  copper  sulfate  electrodes  appear  to  be  technically  acceptable 
for  SP  monitoring  applications.  Their  major  drawback  in  comparison  with  CCS 
electrodes  is  their  considerably  higher  initial  cost  (about  $20  to  $30  each 
compared  with  a  few  dollars  for  a  CCS  stake)  and  additional  effort  involved  in 
installation  (burial  and  careful  insulation  of  connections) . 

59.  For  conventional  streaming  potential  measurements,  CCS  electrodes 
appear  to  have  a  S/N  ratio  considerably  lower  than  that  of  copper  sulfate 
electrodes.  Field  and  laboratory  measurements  indicate  that  responses  of  CCS 
electrodes  to  environmental  disturbances  are  one  or  two  orders  of  magnitude 
greater  than  those  for  copper  sulfate  electrodes.  This  higher  noise  level  is 
due  both  to  the  direct  exposure  of  the  metal  electrode  to  the  soil  (as  opposed 
to  the  isolation  provided  by  the  construction  of  a  nonpolarizing  electrode)  as 
well  as  the  exposure  of  the  unburied  portion  of  the  stake  to  solar  heating  and 
rainfall . 

60.  As  streaming  potential  response  for  CCS  electrodes  isolated  from 
direct  contact  with  flowing  fluid  is  about  the  same  as  that  for  copper 
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sulfate,  there  would  appear  to  be  no  technical  advantage  to  using  CCS  elec¬ 
trodes  for  SP  monitoring.  However,  the  apparent  correlations  between  subsur¬ 
face  seepage  flow  and  SP  anomalies  of  more  than  100  mV  measured  by  field 
arrays  of  CCS  electrodes  imply  that  some  secondary  mechanism  may  generate 
these  large-amplitude  SP  anomalies.  Such  mechanisms  might  include  soil  prop¬ 
erty  variations  and/or  surface  water  movement  associated  with  faults  or  frac¬ 
tures  that  extend  from  the  surface  to  a  seepage  path  at  depth.  Additional 
research,  possibly  including  field  measurements  of  soil  properties  at  elec¬ 
trode  stations  for  a  CCS  installation,  would  be  necessary  to  determine  whether 
such  mechanisms  exist. 
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PART  III:  SELF-POTENTIAL  INVESTIGATIONS  AT 
BEAVER  LAKE  DAM,  ARKANSAS 


Introduction 


61.  This  section  describes  results  of  geophysical  investigations  per¬ 
formed  at  Beaver  Lake  Dam,  Arkansas,  at  pool  elevations  of  1,116  ft  (August 
1986)  and  1,120  ft  (February  1987).  The  geotechnical  problem  at  the  Beaver 
Dam  site  is  anomalous  underseepage  in  the  foundation  of  Dike  1.  The  dike  is 
constructed  across  a  graben  with  vertical  displacements  of  200  ft  (61  m) .  A 
detailed  description  of  the  site  and  the  geotechnical  problem  is  presented  in 
Appendix  C.  The  purposes  oi  the  present  investigations  were  (a)  to  determine 
the  relationship  between  SP  readings  and  seepage  flow  rates  and  patterns;  and 
(b)  to  study  the  performance  of  copper-copper  sulfate  and  copper-clad  steel 
electrodes. 

62.  Investigations  performed  at  the  field  site  included: 

a.  Measurement  of  SP  on  permanently  installed  copper-clad 
steel  (CCS)  stake  electrodes  (referenced  to  the  CCS  base 
electrode  at  the  north  end  of  line  B) .  Note  that  no  measure¬ 
ments  were  made  for  line  A,  on  which  the  electrodes  were 
installed  under  water,  on  the  floor  of  the  reservoir. 

b.  Measurement  of  SP  profiles  using  copper  sulfate  electrodes 
along  stake  lines  B,  C,  D,  and  E,  as  well  as  additional  sta¬ 
tions  in  areas  of  interest. 

c.  Acquisition  of  electrical  resistivity  data  at  three  stations 
along  line  C  (February  1987)  for  comparison  with  previous  data. 

d.  Monitoring  of  SP  noise  levels  between  CCS  stake  and  copper 
sulfate  electrode  pairs  over  periods  of  several  days  to  obtain 
information  on  time  variation  of  SP  readings  unrelated  to  seep¬ 
age  flow. 

e.  Field  and  laboratory  measurements  of  water  and  soil  resistivi¬ 
ties  to  determine  their  effect  on  the  SP  readings. 

Results  of  these  investigations  are  described  in  the  following  sections. 


Summary  of  Results 

63.  SP  contours  from  data  measured  at  low  pool  (1,116  ft)  in  August 
1986  and  at  high  pool  (1,120  ft)  in  February  1987  showed  negative  anomalies 
associated  with  subsurface  seepage  flow  paths  and  positive  anomalies  associ¬ 
ated  with  areas  of  seepage  outflow.  SP  anomaly  patterns  were  influenced  by 
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topographic  variations,  but  topographic  effects  do  not  appear  to  make  major 
contribution  to  the  observed  SP  anomalies. 

64.  Changes  between  the  low  and  high  pool  measurements  included  a 
reduction  of  most  anomaly  amplitudes  at  high  pool  and  significant  shifts  of 
anomaly  patterns  at  several  locations.  The  reduced  anomaly  amplitudes  proba¬ 
bly  were  associated  with  reduced  subsurface  resistivity  at  high  pool  caused  by 
increased  poor  water  saturation  levels.  The  shifted  SP  anomaly  patterns 
appear  to  be  associated  with  changes  in  seepage  flow  patterns. 

65.  Studies  of  electrode  performance  and  background  noise  levels  at  the 
Beaver  Dam  site  indicate  that  buried  nonpolarizing  electrodes  give  satisfac¬ 
tory  S/N  levels  for  long-term  monitoring  of  SP  variations  associated  with 
seepage  flow.  As  discussed  in  Part  II,  less  expensive  CCS  stake  electrodes 
also  may  prove  satisfactory  for  this  purpose  if  their  possible  greater 
response  to  SP  seepage  anomalies  offsets  their  observed  higher  noise  level 
caused  by  variations  in  environmental  conditions  such  as  temperature  and  soil 
moisture  content. 

66.  Measured  S/N  ratios  for  CCS  electrodes  appear  to  be  lower  than 
those  for  copper  sulfate  electrodes,  and  the  two  electrode  types  measured  dif¬ 
fering  SP  profiles  along  the  same  survey  lines.  The  CCS  data  profiles  were 
not  simply  "amplified"  versions  of  the  copper  sulfate  data  profiles,  but, 
instead,  the  two  electrode  types  appear  to  be  responding  to  different  input 
parameters . 


Results  of  Field  Measurements 


Self-potential  measure- 

ments  on  copper  stake  electrodes 

67.  A  network  of  CCS  was  installed  at  the  site  by  WES  and  Little  Rock 
District  (LRD)  personnel  in  March  1985.  Installation  details  and  the  subse¬ 
quent  measurement  program  for  these  stake  electrodes  are  described  in  the  Dam 
Safety  Assurance  Program  Reconnaissance  Report,  Supplement  No.  1,  of  April 
1986;  the  Feature  Design  Memorandum  (FDM)  of  September  1987,  and  a  thesis  by 
Llopis  (1987). 

68.  As  part  of  the  present  program,  SP  measurements  were  made  on  the 
CCS  stake  electrodes  on  15  February  1987.  Because  previous  construction  at 
the  site  had  resulted  in  considerable  damage  to  the  permanently  installed 
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connecting  cables,  measurements  were  made  using  a  reel  of  insulated  wire  con¬ 
nected  between  the  stake  base  electrode  at  the  north  end  of  line  B  and  the 
stake  to  be  measured.  Readings  were  taken  using  a  Fluke  Model  8020A  digital 
multimeter.  Results  are  summarized  in  Table  1  and  are  discussed  in  the  fol¬ 
lowing  sections. 

Self-potential  measurements 

using  copper-copper  sulfate  electrodes 

69.  Results  of  measurements  performed  at  low  pool  (1,116  ft)  in  August 
1986  are  described  in  the  final  report  for  Phase  I  of  this  study.*  Field 
procedures  and  equipment  for  these  measurements  and  for  those  taken  at  high 
pool  are  described  in  Appendix  B. 

70.  The  measuring  stations  used  for  the  August  1986  survey  were  reoc¬ 
cupied  in  February  1987  (with  some  minor  changes  of  station  locations  in  the 
southeastern  portion  of  the  • urvey  area) .  Station  locations  and  contoured 
data  for  the  August  1986  and  February  1987  surveys  are  shown  in  Plates  7 
and  8,  respectively. 

71.  In  later  discussions,  line  S-l-S-3  refers  to  a  north-south  SP  pro¬ 
file  that  runs  through  the  locations  of  Piezometers  S-l  and  S-3,  extending 
from  about  150-ft  north  of  S-l  and  a  similar  distance  south  of  S-3.  Line  E4.5 
refers  to  a  profile  that  begins  at  sta  E5  and  extends  about  450  ft  to  the 
south  (Plate  8).  Portions  of  both  of  these  lines  were  measured  in  August 
1986,  and  complete  data  profiles  were  taken  in  February  1987. 

72.  SP  profiles  for  the  two  copper  sulfate  electrode  surveys,  together 
with  topography  along  each  north-south  survey  line,  are  shown  in  Plate  9;  and 
Plate  10  shows  the  differences  between  the  readings  for  the  two  surveys 
(February  1987  values  minus  August  1986  values)  along  each  north-south  survey 
line. 

Measurement  of  time  varia- 
tlons  of  self-potential  readings 

73.  The  report  for  Phase  I  of  this  study  described  measurements  of  time 
variations  across  selected  pairs  of  CCS  stake  and  copper  sulfate  electrodes 

*  Corwin,  Robert  F. ,  1986,  "Development  of  Self  Potential  Interpretation 
Techniques  for  Seepage  Detection,"  Phase  I  Final  Letter  Report,  Contract 
DACW39-86-C-0059,  US  Army  Engineer  Waterways  Experiment  Station, 

Vicksburg,  MS. 
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Table  1 

Self-Potential  Measurements  on  Copper-Clad  Steel  Stake  Electrodes 

15  February  1987 
Pool  Level  1,120  ft 


Line  B  Line  C  Line  C  Line  D  Line  E 


Station 

(mV) 

Station 

(mV) 

Station 

(mV) 

Station 

(mV) 

Station 

(mV) 

1 

-361 

1 

+13 

36 

-256 

1 

+9 

1 

* 

2 

-555 

2 

+32 

37 

-18 

2 

i k 

2 

* 

3 

-478 

3 

-12 

38 

-148 

3 

-17 

3 

* 

4 

-522 

4 

+11 

39 

-33 

4 

-90 

4 

* 

5 

-544 

5 

+34 

40 

-66 

5 

* 

5 

-8 

6 

* 

6 

-116 

41 

-136 

6 

-90 

6 

* 

7 

-544 

7 

+  18 

42 

-137 

7 

-58 

7 

-86 

8 

-561 

8 

-170 

43 

-37 

8 

-69 

8 

+42 

9 

-532 

9 

-93 

44 

-441 

9 

* 

10 

-599 

10 

-26 

45 

-497 

10 

+14 

11 

-567 

11 

-18 

46 

-324 

11 

-84 

12 

-6 

12 

-12 

47 

-7 

12 

+6 

13 

-513 

13 

-56 

48 

-471 

13 

-10 

14 

-217 

14 

-381 

49 

-207 

14 

+24 

15 

* 

15 

-17 

50 

+  16 

15 

* 

16 

-555 

16 

-392 

51 

-31 

16 

* 

17 

-509 

17 

* 

52 

+21 

17 

* 

18 

-409 

18 

-279 

53 

-27 

19 

-318 

19 

-24 

54 

-578 

20 

-477 

20 

-423 

55 

-314 

21 

-414 

21 

-232 

56 

-93 

22 

-443 

22 

-317 

57 

-8 

23 

-520 

23 

* 

58 

-34 

24 

-442 

24 

-367 

59 

* 

25 

-244 

25 

-436 

60 

* 

26 

-85 

26 

-489 

61 

-26 

27 

-8 

27 

-85 

62 

+  10 

28 

-28 

28 

-306 

29 

-32 

29 

-39 

30 

-34 

30 

-306 

31 

-23 

31 

-226 

32  -239 

33  -476 

34  -354 

35  -596 


Note:  Reference  at  stake  base  electrode  at  north  end  of  line  B. 
*  Indicates  missing  stake. 
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performed  in  August  1986.  The  readings  across  the  CCS  stake  electrode  were 
made  using  the  permanent  cable  installed  for  those  electrodes,  and  there  was 
some  question  as  to  whether  observed  variations  associated  with  rainfall  might 
have  been  due  to  grounding  of  the  cable.  Therefore,  similar  measurements  were 
made  in  February  1987  using  an  independent  insulated  cable  between  stakes  C42 
and  C62. 

74.  The  potential  between  these  stakes  was  monitored  continuously  on  a 
strip  chart  recorder  from  1,000  on  15  February  1987  to  2,000  on  18  February 
1987.  The  strip  chart  record  is  shown  in  Plate  11.  The  significance  of  the 
recorded  data  is  discussed  below. 

Electrical  resistivity  soundings 

75.  During  the  course  of  the  field  investigation,  it  became  apparent 
that  information  about  possible  changes  of  subsurface  electrical  resistivity 
vales  would  be  needed  to  properly  interpret  the  SP  data.  Resistivity  data  had 
been  obtained  by  WES  personnel  in  August  of  1986  at  sta  C21,  sta  C27,  and 

sta  C33.  For  comparison  with  these  data,  vertical  electric  soundings  were 
taken  at  sta  C21,  sta  C28,  and  sta  C33  in  February  1987  (C27  could  not  be 
reoccupied  due  to  piezometer  installation). 

76.  The  data  were  taken  using  a  Soiltest  "Strata  Scout"  Model  R-40C 
resistivity  meter  (20  mA  maximum  output).  The  Schlumberger  array  was  used  at 
sta  C28  and  the  Wenner  array  at  sta  C21  and  sta  C33. 

77.  Observed  data  and  interpreted  layering,  determined  with  an  auto¬ 
mated  computer  inversion  program,  are  summarized  in  Table  2.  Because  the 
instrument  signal  level  was  marginal  for  the  desired  depth  of  investigation, 
interpreted  layering  should  be  considered  approximate.  Interpreted  layering 
for  the  August  1986  and  February  1987  measurements  are  compared  in  Plate  12. 
The  significance  of  the  observed  resistivity  variations  is  discussed  below. 

For  a  discussion  of  the  electrical  method,  see  Department  of  the  Army  (1979) 
and  Butler  et  al.  (1982). 

Measurements  of 
ground-water  resistivity 

78.  Because  soil  resistivity  and  SP  anomaly  source  strength  both  depend 
strongly  on  the  resistivity  of  ground  water,  both  in  situ  and  laboratory 
measurements  of  ground-water  resistivity  were  performed.  Measurements  of 
temperature  and  resistivity  for  reservoir  surface  water  and  seepage  flow  water 
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Table  2 


Measured 

Data  and  Interpreted  Layering  for 

Vertical 

Electric  Soundings, 

,  February  1987 

Measured 

Data 

Station  C21 

Station 

C28 

Station  C33 

(Wenner  Array) 

(Schlumberger  Array) 

(Wenner  Array) 

a  pa 

AB/2 

pa 

a 

pa 

(ft)  (ohm-ft) 

(ft) 

(ohm-ft) 

(ft) 

(ohm-ft) 

5  213 

4.64 

189 

5 

116 

10  181 

6.81 

234 

10 

147 

20  152 

10.0 

284 

16 

190 

30  171 

14.7 

329 

20 

214 

50  195 

21.5 

391 

30 

294 

100  185 

31,6 

454 

50 

390 

46 .4 

484 

64 

373 

68.1 

426 

100 

520 

100.0 

247 

147.0 

122 

Interpreted  Layering 

C21 

C28 

C33 

p^  (ohm  ft) 

230.0 

173.0 

110.0 

t1  (ft) 

8.1 

4 . 6 

10.2 

d1  (ft) 

8.1 

4.6 

10.2 

P2  (ohm-ft) 

41.6 

728.0 

916.0 

t2  (ft) 

5.0 

31.8 

20.7 

d2  (ft) 

13.1 

46.4 

30.9 

P3  (ohm-ft) 

773.0 

22.8 

144.0 

t3  (ft) 

11.0 

26.1 

d3  (ft) 

24.1 

57.0 

p^  (ohm-ft) 

136.0 

5,000.0 
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from  exit  6  were  taken  on  19  February  1987  for  comparison  with  previous  data 
taken  by  WES  personnel.  Samples  of  each  of  these  waters  were  taken  and 
laboratory  measurements  of  conductivity  (the  inverse  of  resistivity)  as  a 
function  of  temperature  were  performed  using  a  Yellow  Springs  Model  33  S-CT 
temperature-conductivity  meter.  Plots  of  measured  conductivity  versus 
temperature  are  shown  in  Plate  13.  These  data  are  discussed  below. 


Interpretation  of  SP  Data  Taken  with 
Copper  Sulfate  Electrodes 


79.  For  purposes  of  this  discussion,  we  will  refer  to  the  data  taken  in 
August  1986  (pool  level  1,116  ft)  as  "low  pool"  data,  and  the  February  1987 
data  (pool  level  1,120  ft)  as  "high  pool"  data. 

80.  Comparison  between  SP  profiles  (Plates  9  and  10)  and  contours 
(Plates  7  and  8)  taken  at  low  and  high  pool  indicates  that,  although  the  gen¬ 
eral  SP  pattern  is  similar  at  low  and  high  pool  levels,  significant  detail 
changes  have  taken  place. 

a.  The  general  level  of  SP  activity  appears  to  be  somewhat  less  at 
high  pool  than  at  low  pool.  This  effect  is  most  clearly  seen 
in  line  D  (Plate  9) ,  where  the  SP  anomaly  profile  shapes  are 
similar  at  high  and  low  pool,  but  peak-to-peak  amplitude  is 
about  100  mV  at  low  pool  and  about  50  mV  at  high  pool.  This 
effect  also  is  seen  at  the  south  end  of  line  B,  as  well  as  in 
the  seepage  areas  around  and  to  the  east  of  exit  6,  where  maxi¬ 
mum  anomaly  amplitudes  are  about  +60  mV  at  low  pool  and  about 
+40  mV  at  high  pool  (exit  6  is  located  about  100-ft  south  of 
sta  D14-D17) . 

b.  From  theoretical  considerations,  we  would  expect  the  amplitude 
of  seepage-related  SP  anomalies  to  increase  with  increasing 
seepage  flow  rates  at  high  pool  level.  Such  increases  in  SP 
anomaly  levels  have  been  observed  at  other  dams  where  measure¬ 
ments  were  made  at  low  and  high  pool  level.  Therefore,  the 
observed  apparent  general  decrease  in  SP  activity  level  at  high 
pool  must  be  explained  to  understand  the  true  relationship 
between  seepage  rates  and  SP  anomaly  amplitude  at  Beaver  Dam. 
Possible  reasons  for  this  observed  decrease  are  discussed  in 
later  sections. 

c.  Significant  shifts  in  locations  of  major  anomaly  features  are 
seen.  For  example,  the  negative  closure  centered  near  sta  C45 
at  low  pool  appears  to  have  migrated  about  100  ft  to  the  north 
at  high  pool,  and  the  negative  closure  centered  near  sta  D5  at 
low  pool  shows  a  considerable  change  of  orientation  at  high 
pool.  A  further  change  is  seen  in  the  northern  portions  of 
lines  B  and  C,  where  high  pool  anomaly  levels  are  more  positive 
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on  this  portion  of  line  B,  but  more  negative  in  the  same  area 
of  line  C.  This  results  in  a  significant  shift 
of  the  contour  pattern  in  this  area. 

81.  In  an  effort  to  interpret  the  observed  profiles  and  changes  in 
terms  of  seepage  flow,  several  aspects  of  the  observed  data  were  investigated. 
These  aspects  included  the  relation  between  SP  anomalies  and  observed  seepage 
patterns  and  geology,  the  effects  of  topography  on  the  SP  readings,  and  the 
relation  between  SP  anomaly  amplitudes  and  soil  and  ground-water  resistivity. 

Relation  between  SP  anom- 

alies,  geology,  and  observed  seepage 

82.  In  the  absence  of  preferred  seepage  paths,  the  SP  pattern  around  an 
impoundment  structure  such  as  a  dam  or  dike  generally  will  show  maximum  nega¬ 
tive  values  on  the  crest  of  the  structure,  with  values  becoming  uniformly  more 
positive  in  the  downstream  direction.  This  positive  downstream  gradient  is 
due  to  the  uniform  seepage  flow  beneath  the  structure.  This  same  "negative 
summit"  effect  also  is  seen  in  areas  of  topographic  relief,  where  downhill 
flow  of  ground  water  often  generates  SP  signals  that  mirror  image  the 
topography. 

83.  In  areas  where  uniform  seepage  is  interrupted  by  channeling  along 
preferential  flow  paths,  negative  anomalies  are  seen  to  follow  along  the  sub¬ 
surface  seepage  flow  paths,  and  positive  anomalies  form  in  areas  where  the 
seepage  flow  is  upward  toward  the  surface.  The  presence  of  subsurface  fea¬ 
tures  such  as  faults,  dikes,  contacts,  or  artificial  drainage  structures  also 
will  have  a  characteristic  effect  on  the  SP  pattern. 

8A.  The  Beaver  Dam  data  at  both  high  and  low  pool  exhibit  the  behavior 
discussed  above.  Negative  anomalies  along  the  southern  fault  zone  correlate 
with  known  subsurface  seepage  flow  paths  and  positive  anomalies  are  seen  in 
the  vicinity  of  the  seepage  exits  in  the  southeastern  portion  of  the  survey 
area.  A  negative  anomaly  centered  near  sta  D5  may  indicate  seepage  along  a 
fault  zone  inferred  from  seismic,  radar,  and  resistivity  data,  with  positive 
anomalies  corresponding  to  the  emergence  of  seepage  water  in  the  south  ravine. 
Although  data  coverage  in  the  vicinity  of  sta  64+00  on  line  C  is  not  suffi¬ 
cient  to  allow  reliable  contouring,  a  strong  SP  gradient  apparent  both  in  the 
contours  and  in  the  profiles  of  Plates  9  and  10  in  this  area  appears  to  be 
related  to  seepage  flow  along  the  southern  boundary  of  the  northern  fault 
zone. 
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85.  Soil  conditions  for  the  August  1986  survey  ranged  from  very  dry 
over  most  of  the  survey  area  to  saturated  in  the  seep  areas.  Because  previous 
studies  have  shown  that  SP  readings  made  with  copper  sulfate  electrodes  often 
become  more  positive  as  soil  moisture  content  increases  (see  Part  II) ,  there 
was  some  concern  that  the  positive  SP  anomalies  seen  at  the  seep  areas  could 
be  caused  primarily  by  the  higher  soil  moisture  content  in  those  areas. 

86.  For  the  February  1987  survey,  previous  rainfall  had  uniformly  satu¬ 
rated  the  surface  soil  throughout  the  survey  area.  Therefore,  the  positive  SP 
anomalies  at  the  seepage  exit  areas  appear  to  be  related  mainly  to  seepage 
flow  rather  than  to  variations  in  soil  moisture  content. 

87.  Results  of  a  brief  study  of  the  source  depths  for  the  SP  anomalies 
along  line  C  are  shown  in  Plates  14  and  15.  The  study  was  conducted  using 
computer  program  SP1.BAS,  discussed  in  Part  IV  and  included  in  Appendix  A. 

88.  Plate  14  shows  the  original  SP  field  data  for  line  C,  taken  in 
August  1986.  In  Plate  15,  the  data  have  been  "reduced"  for  interpretation 
purposes  by  removing  a  least-squaros  trend  line  and  subtracting  a  constant 
shift  of  10  mV.  This  allowed  a  better  visual  fit  between  the  field  data  and 
the  calculated  model  curve. 

89.  The  model  parameters  are  listed  in  Plate  15.  The  model  consisted 
of  four  negative  line  sources  and  one  positive  line  source,  all  trending 
east-west  and  100  ft  in  length.  The  source  at  x  =  6,398  ft  is  at  a  depth  of 
2  ft,  representing  a  near-surface  effect  that  may  be  one  of  the  boundaries  of 
the  northern  fault  zone.  The  negative  sources  at  x  =  6,900,  7,250,  and 
7,650  ft  appear  to  represent  zones  of  horizontal  seepage  ilow,  while  the 
positive  source  at  x  =  7,458  ft  could  represent  a  zone  of  upward  water  flow. 
The  maximum  depth  of  about  50  ft  for  these  four  sources  places  them  close  to 
or  somewhat  above  the  high-velocity  bedrock  shown  in  Figure  31  of  Llopis 
(1987) . 

90.  The  locations  of  these  SP  sources  are  in  good  general  agreement 
with  high-resistivity  peaks  shown  in  Figure  18  of  Llopis  (1987)  and  with  the 
locations  of  fracture  zones  inferred  from  seismic  and  radar  data  discussed  in 
Appendix  A  of  the  Feature  Design  Memorandum  (FDM)  of  September  1987.  As  dis¬ 
cussed  in  Part  IV  of  this  report,  a  more  detailed  SP  modeling  effort  using 
program  SPXCPL  presently  is  being  conducted. 

91.  Thus,  regardless  of  the  relations  between  SP  anomaly  amplitude  and 
seepage  flow  rate,  and  between  SP  patterns  and  topography  as  discussed  below, 
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the  contoured  SP  data  at  Beaver  Dam  appear  to  be  strongly  correlated  with  know 
seepage  flow  paths.  These  correlations  are  discussed  in  more  detail  in  Appen¬ 
dix  A  of  the  September  1987  FDM. 

Topographic  effects 

92.  The  effects  of  topography  on  SP  readings  are  discussed  briefly 
above.  Inspection  of  Plate  9  shows  the  mirror  image  effect  between  SP  and 
topographic  profiles  along  lines  D  and  S1-S3,  and  between  sta  7300  and 

sta  7500  on  line  C.  However,  it  is  of  considerable  interest  that  the  opposite 
effect  is  seen  at  the  southern  ends  of  lines  C  and  E4.5  and  at  the  northern 
ends  of  lines  B  and  C. 

93.  Close  inspection  of  the  profiles  for  lines  D  and  S1-S3  indicates 
that  the  peaks  of  the  topographic  and  SP  profiles  are  offset  from  each  other. 
Cross-plots  of  SP  versus  elevation  on  these  lines  (Plates  16  and  17)  indicate 
the  expected  negative  correlation  between  SP  and  elevation,  but  the  correla¬ 
tion  is  not  as  strong  as  would  be  expected  if  the  observed  anomalies  were  due 
only  to  topographic  effects. 

94.  This  weak  correlation,  along  with  the  positive  correlation  seen  on 
other  lines,  suggests  that  although  topographic  effects  contribute  to  observed 
SP  anomalies  at  Beaver  Dam,  they  are  not  necessarily  the  major  source  of 
anomalies  even  where  the  topography  is  steep. 

Effect  of  ground-water  resistivity 

95.  The  magnitude  of  the  streaming  potential  generated  by  seepage  flow 
is  directly  proportional  to  the  resistivity  of  the  flowing  ground  water. 
Therefore,  with  other  conditions  equal,  SP  anomaly  amplitudes  will  increase 
linearly  as  ground-water  resistivity  increases.  Also,  for  a  given  source 
amplitude,  SP  anomalies  measured  at  the  surface  will  increase  linearly  with 
the  resistivity  of  the  soil  between  the  seepage  zone  and  the  surface. 

96.  The  soil  and  water  resistivity  measurements  described  previously 
were  made  to  determine  whether  these  resistivity  values  had  changed  between 
low  and  high  pool  levels.  As  shown  in  Plate  12,  earth  resistivity  values 
changed  significantly  between  the  low  pool  measurements  in  August  1986  and  the 
high  pool  measurements  in  February  1987.  At  depths  of  less  than  5  to  10  ft, 
resistivities  generally  were  somewhat  less  at  low  pool  than  at  high  pool.  At 
greater  depths,  there  is  considerable  variability  in  the  observed  resistivity 
changes.  Particularly  significant  are  the  decreases  at  high  pool  from  12,700 
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to  136  ohm-ft  for  the  basement  resistivity  value  at  sta  C21  and  from  3,700-  to 
22.8-ohm-ft  at  sta  C27-C28. 

97.  Previous  repeated  earth  resistivity  measurements  in  other  areas 
show  little  seasonal  variation  (Morrison,  Corwin,  and  Chang  1977).  Therefore, 
it  appears  that  the  observed  resistivity  variations  at  Beaver  Dam  are  related 
to  the  seepage  flow  at  depths  of  about  30-  to  50-ft.  Note  that  sta  C33  and 
C27-C28,  where  the  greatest  resistivity  decreases  at  depth  were  observed,  are 
just  upstream  of  the  large  negative  SP  anomaly  that  was  previously  interpreted 
as  related  to  seepage  flow  along  a  fault  inferred  by  seismic  reflection, 
radar,  and  resistivity  data.  Also,  this  is  an  area  where  a  significant  change 
in  the  SP  anomaly  contour  pattern  was  observed  between  low  and  high  pool. 

98.  There  are  two  possible  explanations  for  the  generally  lower 
observed  resistivities  at  high  pool:  increased  pore  saturation  due  to  the 
greater  seepage  flow,  and  decreased  pore  water  resistivity  due  to  changes  in 
the  temperature  of  the  seepage  water.  To  investigate  the  effect  of  tempera¬ 
ture  changes  on  pore  water  resistivity,  laboratory  measurements  as  described 
previously  were  performed. 

99.  The  data  shown  in  Plate  13  indicate  that,  as  expected  from  theore¬ 
tical  considerations,  water  conductivity  increases  with  increasing  tempera¬ 
ture.  (Note  chat  resistivity  is  the  inverse  of  conductivity,  and  resistivity 
in  ohm-ft  is  equal  to  32,800  divided  by  conductivity  in  micromhos/cm.)  Thus, 
from  the  data  in  Plate  13,  seepage  water  resistivity  at  a  typical  winter  tem¬ 
perature  of  10°  C  is  about  300  ohm-ft,  and  decreases  to  about  230  ohm-ft  at  a 
typical  summer  water  temperature  of  20°  C. 

100.  As  this  effect  is  opposite  from  the  observed  general  decrease  in 
soil  resistivity  during  the  winter,  and  as  the  observed  soil  resistivity 
changes  are  considerably  greater  than  the  100-ohm-ft  range  expected  from  tem¬ 
perature  changes,  our  conclusion  is  that  the  observed  soil  resistivity  changes 
are  due  primarily  to  changes  in  levels  of  pore  water  saturation  caused  by 
infiltration  of  surface  precipitation  and  by  changes  in  seepage  flow  rates. 

101.  Considerable  work  remains  to  be  done  in  relating  the  observed 
resistivity  variations  between  low  and  high  pool  to  piezometric  profiles  and 
seepage  patterns.  However,  the  observed  resistivity  variations  have  two 
important  implications.  First,  it  appears  that  periodic  measurements  of 
resistivity  variations  could  provide  a  useful  technique  for  monitoring  pore 
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saturation  levels  and  seepage  patterns.  Second,  the  observed  resistivity 
decreases  are  large  enough  to  readily  account  for  the  observed  general 
decrease  in  SP  anomaly  amplitude  at  high  pool  levels.  Even  though  SP  anomaly 
source  amplitudes  would  increase  due  to  the  higher  seepage  flow  rates  at  high 
pool  and  the  higher  water  resistivity  during  the  winter,  the  effect  of  greatly 
reduced  soil  resistivity  caused  by  increased  pore  water  saturation  could  cause 
a  reduction  in  the  SP  anomaly  amplitude  measured  at  the  earth's  surface. 

Data  Acquisition 

102.  In  the  previous  discussions  it  was  assumed  that  the  measured  SP 
data  represented  the  actual  electric  field  at  the  surface  of  the  earth. 

This  is  strictly  true  only  if  no  errors  are  contributed  by  the  measuring 
electrodes,  and  if  short-period  SP  variations  unrelated  to  seepage  flow 
either  are  insignificant  relative  to  seepage  anomalies  or  can  be  corrected  as 
part  of  the  data  reduction  process.  Investigations  of  both  of  these  possible 
error  sources  are  discussed  below. 

Short-period  variations 
of  the  earth's  electric  field 

103.  Because  seepage-related  SP  anomalies  would  be  expected  to  vary 
over  periods  of  several  days  or  more,  any  SF  variation  occurring  over  a 
period  of  a  few  days  or  less  would  be  considered  as  short-period  for  the  pur¬ 
poses  of  this  study.  Also  defined  as  short-period  variations  are  those  not 
related  to  the  electrode  effects  discussed  below,  but  instead  to  actual  vari¬ 
ations  of  the  earth's  electric  field.  Understanding  of  short-period  varia¬ 
tions  is  important  for  SP  seepage^monitoring  studies  because  such  variations 
are  superimposed  on  the  anomalous '’'P  field  generated  by  seepage  flow,  and 
thus  are  a  source  of  noise  during  the  measuring  process. 

104.  Short-period  variations  have  two  sources:  natural  and  artificial. 
Natural  variations  are  caused  by  time  fluctuations  of  the  earth's  magnetic 
field  that  generate  electric  currents  in  the  earth  through  electromagnetic 
coupling.  These  currents  in  turn  generate  electric  potentials  in  the  earth 
that  are  called  telluric  voltages.  Periods  of  these  telluric  voltages  range 
from  less  than  1  sec  up  to  hours  or  days.  The  magnitude  of  the  telluric 
voltage  field  depends  on  the  resistivity  of  the  earth  beneath  the  measuring 
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point  as  well  as  the  magnitude  of  the  magnetic  field  fluctuations.  Telluric 
voltages  in  a  uniform  earth  will  increase  linearly  with  increasing  separation 
between  the  measuring  points. 

105.  The  measurements  at  Beaver  Dam  indicate  that  telluric  variations 
with  periods  of  about  10  sec  to  several  minutes  can  range  up  to  about  ±10  mV 
per  1,000  ft  of  electrode  separation.  Examples  of  such  variations  are  seen 
in  Plate  11  between  0830  and  1200  on  16  February,  and  at  about  1800  on 

18  February.  A  longer-period  telluric  variation  of  about  1-hr  duration  and 
10-mV  amplitude  is  seen  between  0300  and  0400  on  18  February.  Information 
such  as  that  shown  in  Plate  11  allows  selection  of  appropriate  measurement 
procedures  and  data  reduction  methods  to  remove  the  effects  of  short-term 
variations  and  allow  measurement  of  the  true  steady-state  SP  field. 

106.  It  should  be  noted  that  lightning  strikes,  even  when  located 
several  hundred  miles  from  the  measuring  site,  can  cause  very  large  short- 
period  fluctuations.  However,  the  period  of  such  fluctuations  is  so  short 
(less  than  1  sec)  that  they  rarely  will  affect  the  SP  measurements.  Many 
examples  of  voltage  "spikes"  generated  by  lightning  strikes  are  seen  in  the 
data  of  Plates  5  and  6  of  the  Phase  I  report. 

107.  Short-period  SP  variations  also  can  be  generated  by  artificial 
sources,  particularly  overhead  power  lines  and  grounds  of  electrical  machin¬ 
ery.  Overhead  powerlines  induce  60-Hz  noise  into  SP  measuring  systems,  often 
causing  drift  and  irregular  fluctuations  of  the  measured  values.  This  prob¬ 
lem  was  present  at  Beaver  Dam,  where  it  was  found  that  60-Hz  noise  severely 
affected  the  strip  chart  recorders  used  for  long-term  monitoring.  Installa¬ 
tion  of  capacitors  across  the  chart  recorder  inputs  greatly  reduced  the  60-Hz 
noise.  This  implies  that  appropriate  input  filtering  should  be  part  of  any 
future  SP  monitoring  system  installed  at  electrically  noisy  sites  such  as 
Beaver  Dam. 

108.  An  example  of  short-period  noise  generated  by  electrical  machinery 
is  shown  in  the  upper  illustration  of  Plate  18.  As  described  in  more  detail 
in  Part  II,  the  measurements  were  made  across  CCS  stake  and  copper  sulfate 
electrode  pairs  separated  by  28  ft  and  located  about  1  mile  from  an  electri¬ 
cal  substation  for  the  electrically  operated  Bay  Area  Rapid  Transit  (BART) 
system  in  the  San  Francisco  Bay  Area  of  California. 

109.  The  noise  level  ranges  from  about  1-  to  5-mV  over  the  28-ft  elec¬ 
trode  separation,  which  extrapolates  to  about  40-  to  200-mV  per  1,000  ft. 
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The  association  of  the  noise  with  operation  of  the  BART  system  is  shown  in  the 
lower  illustration  of  Plate  18,  Noise  levels  during  the  nonoperational  period 
from  about  1  a.m.  to  4  a.m.  are  much  lower  than  those  during  the  day.  Even 
with  this  very  high  noise  level,  average  daily  potentials  for  both  electrode 
pairs  are  stable  (during  good  weather  conditions)  and  with  appropriate 
filtering  of  short-period  variations  it  would  be  possible  to  monitor  long-term 
seepage  potentials  even  under  such  noisy  conditions. 

110.  Note  that  in  Plate  18,  the  response  of  both  the  copper  sulfate  and 
the  CCS  stake  electrodes  to  short-period  variations  is  virtually  identical. 
This  contrasts  with  their  differing  response  to  changes  in  environmental 
conditions  below  and  in  Part  II. 

Electrode  effects 

111.  General  discussion.  The  steady-state  potential  measured  between 
an  electrode  pair  in  the  earth  is  the  sum  of  four  major  contributions: 

a.  Polarization  potentials  caused  by  electrochemical  differences 
between  the  two  electrodes  (this  potential  may  be  constant  or 
may  change  with  time) . 

b.  Environmental  effects  caused  by  different  temperature  or  soil 
conditions  (moisture  content,  chemistry,  etc.)  at  the  two 
electrode  locations. 

c.  Artificial  sources. 

d.  The  anomaly,  if  present,  that  is  to  be  measured. 

Additionally,  installation  defects  such  as  a  grounded  connecting  cable  or  poor 
contact  between  an  electrode  and  the  soil  can  contribute  large  offsets  to  the 
measured  potentials. 

112.  Electrode  effects  for  "nonpolarizing"  electrodes  (in  which  the 
metal  sensing  element  is  isolated  from  the  soil  in  a  bath  of  solution  contain¬ 
ing  a  salt  of  the  metal  element;  e.g.,  CCS)  have  been  studied  extensively  (see 
Part  II).  Results  of  these  studies  indicate  that  long-term  variations  between 
pairs  of  nonpolarizing  electrodes  buried  in  the  earth  are  only  a  few  mV,  and 
generally  do  not  exceed  10  mV.  A  sample  record  between  a  pair  of  copper  sul¬ 
fate  electrodes  installed  at  Beaver  Dam  is  shown  in  Plate  5  of  the  Phase  I 
report.  Maximum  long-period  drift  between  these  electrodes  was  about  5  mV. 

113.  Electrode  effects  for  CCS  stake  electrodes  have  not  previously 
been  studied  in  detail.  It  is  known  that  the  magnitude  of  point-to-point 
potential  variations  measured  with  CCS  stake  electrodes  is  much  larger  than 
that  measured  by  nonpolarizing  electrodes  (for  example,  see  Plates  7  and  8  of 
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the  Phase  I  report).  The  contribution  of  electrode  effects  to  these  varia¬ 
tions  is  discussed  in  Part  II. 

114.  Long-term  measurements  (over  periods  of  several  hours  to  a  few 
days)  between  pairs  of  CCS  stake  electrodes  at  Beaver  Dam  are  described  in  the 
Phase  I  report.  Results  of  the  studies  of  Part  II  indicate  that  much  of  the 
observed  variation  was  caused  by  diurnal  temperature  effects.  As  discussed  in 
Part  II,  rainfall  appears  to  have  a  very  significant  effect  on  the  potential 
between  a  pair  of  CCS  stake  electrodes.  Plate  6  of  the  Phase  I  report  shows 
an  offset  of  more  than  50  mV  between  stakes  at  sta  El  and  sta  E8  that  appears 
to  be  related  to  rainfall.  Because  this  offset  could  possibly  have  been 
caused  by  grounding  of  the  connecting  cable  between  the  electrodes,  a  similar 
measurement  was  made  in  February  1987,  using  an  independent  insulated  connect¬ 
ing  cable. 

115.  Results  of  this  measurement  are  shown  in  Plate  11  of  this  report. 

A  strip  chart  recorder  connected  between  stakes  at  sta  C42  and  sta  C62  showed 
a  change  from  about  -120  mV  to  about  -200  mV  occurring  at  2200,  15  February. 
The  change  followed  a  moderate  1-day  rainstorm  that  ended  a  few  hours  before 
the  variation.  Following  the  offset,  the  potential  slowly  drifted  back  to 
about  -140  mV  over  a  period  of  about  3  days,  when  the  measurement  was 
terminated. 

116.  Similar  rainfall-related  variations  have  been  observed  on  a  CCS 
stake  pair  installed  as  part  of  a  long-term  monitoring  experiment  (the 
experimental  installation  also  includes  pairs  of  copper  sulfate  and  buried 
metallic  lead  (Pb)  electrodes) .  The  results  of  this  monitoring  experiment  are 
described  in  Part  II. 

117.  These  results  indicate  that  electrode  effects  due  to  temperature 
and  soil  moisture  variations  are  considerably  larger  for  CCS  stake  electrodes 
than  for  copper  sulfate  electrodes.  However,  there  is  some  evidence  that  sig¬ 
nal  levels  (generated  by  seepage  flow)  may  also  be  greater  for  CCS  electrodes 
(see  Plates  7  and  8  of  the  Phase  I  report) .  A  field  study  of  this  possible 
"amplification"  effect  is  discussed  below. 

118.  Comparison  of  copper  sulfate  and  CCS  data.  A  major  consideration 
in  the  selection  of  electrodes  for  long-term  seepage  monitoring  installations 
is  the  question  of  whether  CCS  electrodes  measure  a  greater  SP  signal  level 
for  a  given  seepage  flow  than  do  nonpolarizing  electrodes  such  as  copper 
sulfate.  Laboratory  studies  of  this  question  are  described  in  Part  II.  Below 
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are  compared  in  detail  some  of  the  data  for  copper  sulfate  and  CCS  electrodes 
obtained  from  the  Beaver  Dam  survey  of  February  1987.  As  described  previously 
in  this  section,  both  the  copper  sulfate  and  the  CCS  data  for  this  survey  were 
taken  using  the  same  field  procedure.  Therefore,  differences  between  the  data 
sets  should  be  due  only  to  electrode  effects. 

119.  Plate  19  shows  unsmoothed  copper  sulfate  and  CCS  data  for  line  C. 
It  is  apparent  that  the  total  range  of  SP  activity  is  about  10  times  greater 
for  the  CCS  electrodes  than  for  the  copper  sulfate  electrodes  (about  700  mV 
versus  about  50  mV) ;  and  that  the  magnitude  of  point-to-point  variations  as  a 
percentage  of  total  range  is  greater  for  the  CCS  electrodes  (i.e.,  the  stan¬ 
dard  deviation  for  the  CCS  data  profile  is  greater) .  There  appear  to  be  some 
areas  of  correlation  between  the  two  profiles,  but  the  variability  of  the  data 
makes  it  difficult  to  be  certain  of  this. 

120.  In  an  effort  to  make  visual  comparison  between  the  two  profiles 
easier,  the  field  data  were  smoothed  using  a  three-point  running  mean.  The 
smoothed  profiles  are  shown  in  Plate  20.  In  general,  the  profiles  appear  to 
be  correlated  poorly,  if  at  all.  The  CCS  profile  shows  a  generally  negative 
central  area  that  is  not  apparent  on  the  copper  sulfate  profile,  and  the  loca¬ 
tions  of  individual  positive  and  negative  peaks  on  the  two  profiles  are  not 
coincident.  The  major  apparent  similarity  between  the  two  smoothed  profiles 
is  the  predominant  spatial  wavelength  of  about  200  ft  in  the  central  portion 
of  the  profiles.  However,  the  negative  anomalies  of  this  wavelength  for  the 
CCS  profile  are  more  sharply  peaked,  even  after  smoothing,  than  those  for  the 
copper  sulfate  profile.  This  suggests  a  shallower  source  depth  for  the  CCS 
variations,  since  the  source  depth  is  directly  related  to  the  width  of  the 
anomaly  at  half  its  maximum  value. 

121.  Comparison  of  unsmoothed  data  profiles  for  line  B  is  shown  in 
Plate  21.  Because  the  data  for  this  line  generally  were  less  noisy  than  those 
for  line  C,  smoothing  did  not  affect  the  visual  correlation  of  the  profiles. 
Although  some  apparent  correlations  are  seen  between  the  two  profiles  (e.g., 
the  positive  trend  at  the  right  side  and  the  relative  low  centered  near 

sta  7200),  in  general  the  data  do  not  appear  to  be  significantly  correlated. 

122.  Based  on  these  comparisons,  the  CCS  profiles  do  not  appear  to 
represent  an  "amplified"  version  of  the  copper  sulfate  profiles.  This  conclu¬ 
sion  is  in  agreement  with  the  results  discussed  in  Part  II,  which  indicate 
that  streaming  potential  response  for  CCS  electrodes  is  no  greater  than  that 


for  copper  sulfate  electrodes.  Because  these  results  also  indicate  that  the 
geologic  noise  level  is  considerably  greater  for  CCS  electrodes,  it  is  con¬ 
cluded  that  the  S/N  ratio  for  copper  sulfate  electrodes  is  greater  than  that 
of  CCS  electrodes  for  seepage  monitoring  applications. 

123.  However,  the  successful  results  of  previous  studies  by  WES  and  by 
Erchul  (1986)  using  CCS  electrodes  for  seepage  monitoring  indicate  that  useful 
data  regarding  seepage  flow  can  be  obtained  from  CCS  electrode  arrays.  There¬ 
fore,  it  is  possible  that  some  mechanism  not  apparent  from  the  present  studies 
actually  does  amplify  the  field  response  of  CCS  electrodes  to  streaming  poten¬ 
tials.  As  discussed  in  Part  II,  more  detailed  laboratory  studies,  along  with 
careful  long-term  monitoring  of  coincident  arrays  of  CCS  and  nonpolarizing 
electrodes  under  known  seepage  conditions,  would  be  necessary  to  determine  any 
such  mechanism. 


Conclusions  of  Beaver  Dam  Study 

124.  The  SP  contours  taken  with  copper  sulfate  electrodes  at  the  site 
are  strongly  related  to  seepage  flow  patterns.  Even  though  the  SP  data  are 
affected  to  some  extent  by  topography  and  by  seasonal  changes  of  soil  and  pore 
water  resistivity,  it  appears  that  negative  SP  anomalies  are  associated  with 
subsurface  seepage  flow  paths  and  that  positive  SP  anomalies  are  associated 
with  areas  of  seepage  outflow. 

125.  Although  the  effects  of  seasonal  subsurface  resistivity  variations 
complicate  the  interpretation  of  SP  changes  due  to  pool  level  variations, 
there  appear  to  be  some  significant  SP  variations  between  pool  levels  of  1,116 
and  1,120  ft  that  could  be  related  to  changes  in  seepage  flow  patterns.  The 
most  important  of  these  are  the  variations  measured  in  the  northern  portions 
of  lines  B  and  C,  seen  in  the  contours  of  Plates  7  and  8  and  the  profiles  of 
Plates  9  and  10.  The  proximity  of  these  variations  to  the  northern  fault  zone 
indicate  that  significant  seepage  flow  variations  between  high  and  low  pool 
levels  may  have  occurred  along  this  fault  zone  or  its  boundaries.  This  inter¬ 
pretation  is  supported  by  the  possible  appearance  of  fines  in  the  seepage  flow 
from  this  area  in  February  1987. 

126.  Other  significant  changes  in  SP  patterns  appear  to  be  shifts  in 
the  locations  of  the  negative  SP  anomalies  centered  near  sta  C45  and  sta  D5. 
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Comparison  of  these  shifts  in  SP  anomaly  locations  with  piezometer  data  should 
establish  whether  they  are  related  to  changes  in  seepage  flow  patterns. 

127.  Survey  results  indicate  that  monitoring  of  subsurface  resistivity 
data  is  necessary  for  complete  interpretation  of  SP  data.  The  resistivity 
data  also  should  prove  very  useful  for  monitoring  of  soil  saturation  changes 
due  to  surface  infiltration  and  seepage  flow. 

128.  Measurements  made  during  the  two  field  studies  have  established 
the  nature  and  magnitude  of  SP  noise  due  to  telluric,  artificial,  and  environ¬ 
mental  sources  at  the  Beaver  Dam  site.  The  S/N  level  for  copper  sulfate  elec¬ 
trodes  clearly  is  large  enough  to  permit  reliable  monitoring  of  long-term  SP 
variations  having  magnitudes  of  greater  than  a  few  millivolts.  As  both  anom¬ 
aly  amplitudes  and  variations  related  to  pool  level  appear  to  exceed  10  mV, 
copper  sulfate  or  similar  nonpolarizing  electrodes  should  be  suitable  for 
long-term  monitoring  of  seepage-related  SP  anomalies. 

129.  CCS  stake  electrodes  are  less  expensive  to  purchase,  install,  and 
maintain  than  are  nonpolarizing  electrodes.  Although  their  noise  response  to 
environmental  changes  such  as  temperature  or  rainfall  is  greater  than  that  of 
nonpolarizing  electrodes,  they  have  been  used  successfully  for  seepage-related 
SP  measurements.  However,  results  from  this  study  indicate  that  the  two  elec¬ 
trode  types  may  be  responding  to  different  input  parameters.  Factors  that 
determine  which  electrode  type  is  more  suitable  from  both  economic  and  techni¬ 
cal  considerations  for  long-term  SP  monitoring  are  discussed  in  more  detail  in 
Part  II. 
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PART  IV:  METHODS  FOR  INTERPRETATION  OF  SEEPAGE-RELATED  SP  DATA 


Introduction 


130.  In  this  section,  a  catalog  and  description  techniques  are  pre¬ 
sented  for  interpretation  of  SP  data  for  seepage  investigations.  It  is 
assumed  that  the  field  data  are  of  good  quality,  and  that  the  effects  of  geo¬ 
logic,  artificial,  topographic,  time-varying,  and  other  noise  sources  have 
been  accounted  for. 

131.  SP  data  may  be  interpreted  qualitatively,  geometrically,  or  quan¬ 
titatively.  The  interpretation  procedure  selected  will  depend  on  the  desired 
goals  of  the  investigation,  the  quality  of  the  field  data,  the  amount  of 
available  additional  geological,  geophysical,  and  hydrologic  data,  and  the 
time  and  computer  resources  available  for  the  interpretation  phase  of  the 
investigation. 

132.  Qualitative  interpretation  involves  preparation  of  data  profiles 
and  contours  and  visual  inspection  of  these  to  look  for  patterns  known  or 
thought  to  be  characteristic  of  seepage  flow  paths.  The  results  of  many  pre¬ 
vious  investigations  cited  in  Part  V  (as  well  as  the  Beaver  Lake  Dam  study 
described  in  Part  III)  and  of  quantitative  studies  using  the  techniques 
described  below  indicate  that  negative  SP  anomalies  often  are  seen  at  areas 
where  seepage  flow  is  entering  the  dam  and  above  seepage  paths  where  flow  is 
horizontal  or  descending;  and  that  positive  anomalies  often  are  seen  above 
areas  where  flow  is  ascending  toward  the  surface  or  where  surface  seepage  is 
occurring. 

133.  Such  qualitative  interpretation  has  proved  useful  in  many  cases 
where  the  SP  data  were  used  primarily  to  indicate  locations  for  more  intensive 
hydrologic  or  geophysical  investigations.  However,  the  use  of  geometric 
interpretation  techniques,  which  require  minimal  additional  effort,  can  help 
to  provide  information  about  flow  path  depth  and  configuration  as  well  as 
location. 

134.  Geometric  interpretation  involves  the  use  of  calculated  curves  and 
contours  generated  by  relatively  simple  SP  source  models  to  match  the  observed 
field  data.  The  available  models  include  polarized  points,  lines,  cylinders, 
spheres,  sheets,  and  other  geometric  forms.  Matching  of  field  data  to  the 
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curves  generated  by  these  sources  can  provide  useful  preliminary  information 
about  the  form,  depth,  and  orientation  of  inferred  seepage  paths. 

135.  Although  no  quantitative  information  about  seepage  flow  rates  is 
provided  by  these  techniques,  they  are  useful  not  only  for  the  source  param¬ 
eters  they  provide  but  also  for  helping  to  eliminate  SP  anomalies  caused  by 
sources  for  which  depth  or  configuration  is  inconsistent  with  known  geologic 
or  hydrologic  information.  Also,  the  preliminary  models  derived  from  these 
techniques  are  useful  as  input  to  the  quantitative  modeling  programs  dis¬ 
cussed  below. 

136.  The  listing  in  the  following  section  summarizes  a  number  of  source 
models  selected  from  the  geophysical  literature.  Because  most  of  the  algo¬ 
rithms  are  relatively  simple,  they  may  be  programmed  on  a  calculator  or  per¬ 
sonal  computer.  Appendix  A  includes  computer  program  SP1.BAS,  written  for  the 
IBM-PC  and  compatible  computers  for  the  calculation  of  anomalies  generated  by 
the  geometric  source  models  described  below,  as  well  as  a  user's  manual,  sam¬ 
ple  output,  and  the  program  listing.  SP1.BAS  can  be  used  directly  for  field 
SP  data  interpretation  using  geometric  source  models. 

137.  Quantitative  interpretation  of  SP  data  may  be  done  using  techni¬ 
ques  described  in  the  references  of  Nourbehecht,  Madden,  Fitterman,  Sill,  and 
others  cited  in  the  following  sections.  These  techniques  involve  the  use  of 
computer  programs  adapted  from  algorithms  originally  developed  for  calculation 
of  potentials  and  apparent  resistivities  for  two-  and  three-dimensional 
resistivity  distributions  in  the  earth.  Input  to  these  programs  include  the 
electrical  resistivity  structure  of  the  region  to  be  modeled,  values  of 
streaming  potential  coupling  coefficients  and  permeability  for  the  region,  and 
the  location  and  intensity  of  pressure  sources  and  sinks  representing  areas  of 
seepage  inflow  and  outflow. 

138.  These  computer  programs  are  complex  to  use  and  require  consider¬ 
able  memory  and  execution  time  compared  with  the  relatively  simple  geometric 
source  models  discussed  above.  In  many  cases  resistivity,  coupling  coeffi¬ 
cient,  or  permeability  values  may  not  be  available  and  must  be  estimated. 
Nevertheless,  these  quantitative  techniques  can  provide  a  powerful  tool  for 
interpretation  of  seepage-related  SP  data.  Unlike  the  simpler  techniques  dis¬ 
cussed  above,  they  can  (a)  account  for  complex  geologic,  electrical,  and 
hydrogeologic  structure;  (b)  distinguish  between  pressure  sources  and  sinks, 
and  (c)  provide  quantitative  estimates  of  seepage  flow  rates  and  velocities. 
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139,  Very  little  information  has  been  published  regarding  the  use  of 
these  computer  programs  (which  originally  were  developed  for  interpretation  of 
geothermal  SP  data)  for  analysis  of  seepage-related  SP  data.  The  subsequent 
section  on  quantitative  modeling  briefly  discusses  the  derivation  and  use  of 
one  &uch  program.  The  following  tabulation  Hots  references  for  each  source 
model. 

References  for  SP  Source  Models 


1.  Point  current  sources  (single  and  multiple) 

Alfano  1962 

Broughton,  Edge,  and  Laby  1931 

Corwin  1976 

Corwin  et  al.  1981 

DeMoully  and  Corwin  1980 

Heiland  1940 

Merkel  1971 

Morrison  et  al.  1978 

Morrison  et  al.  1979a 

Morrison  et  al.  1979b 

Paul  et  al.  1965 

Paul  and  Banerjee  1970 

Semenov  1974 

Stern  1945 

Telford  et  al.  1976 

Van  Nostrand  and  Cook  1976 

2.  Horizontal  lines  sources  (including  polarized  sheet 

sources  modeled  as  dipolar  line  pairs) 

Banerjee  1970 

Broughton,  Edge,  and  Laby  1931 

Laxman  et  al.  1986 

Meiser  1962 

Murty  et  al.  1985 

Paul  1965 

Rao  et  al.  1970 

Rao  et  al.  1983 

Roy  and  Chowdhury  1959 

Semenov  1974 

3.  Spherical  sources 

Bhattacharya  and  Roy  1981 
de  Witte  1948 
Heiland  1940 

lakubovskii  and  Liajov  1980 
Muoi  and  Quynh  1988 
Petrowsky  1928 
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Rao  et  al.  1970 
Semenov  1974 
Telford  et  al.  1976 
Yungul  1945 
Yungul  1950 

4.  Cylindrical  sources 

Bhattacharya  and  Roy  1981 
Murty  et  al.  1985 
Semenov  1974 

5.  Dipolar  sheet  source 

Fitterman  1979a 
Fitterman  and  Corwin  1982 
Fitterman  1984 

6.  Quantitative  modeling 

a.  Specific  for  SP  interpretation 

Fitterman  1976 
Fitterman  1978 
Fitterman  1979a 
Fitterman  1979b 
Fitterman  1979c 
Fitterman  1982a 
Fitterman  and  Corwin  1982 
Fitterman  1983a 
Fitterman  1984 
Harding  1981 
Hulse  1978 

Ishido  and  Mizutani  1981 
Nourbehecht  1963 
Nourbehecht  and  Madden  1970 
Sill  and  Johng  1979 
Sill  1981a 

Sill  and  Killpack  1982 

Sill  1982a 

Sill  1982b 

Sill  1982c 

Sill  1983a 

Sill  1983b 

b.  Fundamental  theory  for  quantitative  modeling 
Denbigh  1951 

Marshall  and  Madden  1959 
Mitchell  1976 
Onsager  1931 
Pourbaix  1949 
Prigogine  1955 
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Geometric  SP  Source  Models 


140.  This  section  summarized  SP  source  models  for  a  variety  of  source 
geometries.  Model  geometry  and  equations  for  calculating  SP  fields  are  shown 
in  Plates  A1  through  A5  of  Appendix  A.  For  metric  calculations  the  most  con¬ 
venient  units  are  milliamps  (mA)  for  current  I,  millivolts  (mV)  for  potential 
voltage  V,  ohm-meters  (ohm-m)  for  resistivity  p  ,  and  meters  (m)  for  length. 
Corresponding  English  units  are  mA,  mV,  ohm-feet  (ohm-ft),  and  feet.  The 
equations  used  for  the  algorithms  of  the  computer  program  SPl.BAS  included  in 
this  report  are  those  of  Plates  A1  through  A5.  A  summary  of  selected  refer¬ 
ences  for  each  geometric  source  type,  as  well  as  for  quantitative  modeling  of 
SP  data  ,  is  given  in  the  preceding  section. 

141.  It  should  be  noted  that  the  references  cited  above  for  geometric 
modeling  usually  include  not  only  derivations  and  equations  fro  calculating 
model  curves  but  also  interpretation  schemes  based  on  the  use  of  anomaly 
wavelengths  and  shapes,  characteristic  curves,  nomograms,  and  a  variety  of 
other  methods.  In  many  cases  use  of  these  interpretation  schemes  prior  to 
curve  matching  using  computer  program  SPl.BAS  in  Appendix  A  can  save  con¬ 
siderable  time  by  providing  reasonable  first  estimates  of  source  parameters; 
and,  in  some  cases,  these  easily  obtained  estimates  may  be  sufficient  for  the 
degree  of  interpretation  desired. 

Point  current  sources 

142.  The  geometry  and  modeling  equation  for  a  point  source  of  current 
in  a  uniform  half-space  is  shown  in  Plate  A1  (a  current  sink  is  defined  as  a 
negative  source).  A  point  source  or  sink  or  multiple  combinations  of  point 
sources  and  sinks  provide  a  powerful  and  flexible  geometry  for  modeling  of  SP 
anomalies.  Any  arbitrary  source  configuration,  with  any  arbitrary  charge 
distribution,  can  be  expressed  as  an  appropriate  spatial  distribution  of  point 
sources  and  sinks. 

143.  A  particularly  useful  application  of  a  single  point  source  or  sink 
model  is  to  provide  a  first  estimate  of  the  depth  to  the  source  of  a  circular 
or  nearly  circular  SP  anomaly.  Because  a  point  source  represents  the  minimum 
possible  source  size,  the  source  depth  of  the  anomaly  can  be  no  greater  than 
that  which  provides  a  reasonable  fit  to  a  point  source  model.  Thus,  fitting  a 
point  source  to  the  observed  data  can  quickly  indicate  the  maximum  source 
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depth.  The  size  of  the  source  region  then  must  increase  as  its  depth 
decreases  from  this  maximum  value. 

144.  The  hclf-wavelength  XH  (the  distance  from  the  origin  at  which  the 
anomaly  is  one-half  of  its  maximum  value)  of  an  anomaly  generated  by  a  point 
source  buried  at  depth  d  is  given  by 

XH  =  3Vd 

This  equation  is  helpful  for  quickly  estimating  the  depth  of  a  point  source. 

145.  As  numerous  analytical  equations  have  been  developed  for  calculat¬ 
ing  the  fields  generated  by  point  sources  in  inhomogeneous  media,  the  use  of 
single  or  multiple  point  sources  allowed  relatively  simple  calculation  of  SP 
fields  in  the  presence  of  geologic  structure  such  as  layers,  contacts,  faults, 
dikes,  etc.  More  complex  two-  or  three-dimensional  structure  may  be  modeled 
using  algorithms  developed  for  resistivity  interpretation. 

146.  Examples  of  the  use  of  multiple  point  sources  to  model  complex 
source  geometry  are  given  in  DeMoully  and  Corwin  (1980)  and  Corwin  et  al. 
(1981).  Morrison  et  al.  (1978)  present  a  computer  program  for  calculating  the 
SP  field  generated  by  an  arbitrary  array  of  point  sources  and  sinks  in  the 
presence  of  a  vertical  resistivity  contact.  The  computer  program  SP1.BAS 
included  in  Appendix  A  presently  calculates  point  source  fields  only  for  a 
uniform  half-space  but  could  be  adapted  relatively  easily  to  handle  more  com¬ 
plex  resistivity  structure. 

Horizontal  line  sources 

147.  The  geometry  and  modeling  equation  for  a  horizontal  line  source  of 
current  are  shown  in  Plate  A2.  The  source  is  parallel  to  the  y-axis,  has  a 
constant  current  I  per  unit  length,  and  is  located  in  a  uniform  half-space. 

148.  The  line  source  represents  the  simplest  geometry  for  modeling 
elongated  SP  anomalies.  As  for  a  point  source,  the  depth  to  a  line  source 
that  fits  the  field  data  represents  the  maximum  possible  source  depth  for  an 
elongated  anomaly. 

149.  More  complex  elongated  source  geometries  may  be  modeled  as  dis¬ 
tributions  of  multiple  line  sources  and  sinks.  In  the  literature,  many  models 
described  as  sheet  sources  actually  are  dipolar  line  pairs  located  along  the 
top  and  bottom  edges  of  the  "sheet".  Such  models  are  widely  used  for  inter¬ 
preting  SP  fields  generated  by  thin,  elongated  mineral  deposits.  For  this 
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study,  a  sheet  model  is  considered  as  one  having  uniform  charge  on  the  faces 
of  the  sheet  rather  than  charge  concentrated  along  the  upper  and  lower  edges. 

150.  The  analytical  equation  in  Plate  A2  is  valid  only  for  horizontal 
lines  having  constant  current  per  unit  length.  A  series  of  closely  spaced 
point  sources  may  be  used  to  approximate  lines  for  which  the  current  distribu¬ 
tion  is  not  constant  or  for  which  the  lines  are  not  horizontal  or  for  use  in 
areas  of  nonuniform  resistivity. 

Spherical  sources 

151.  The  geometry  and  modeling  equation  for  a  spherical  source  are 
shown  in  Plate  A3.  The  axis  of  polarization  of  the  sphere  is  inclined  at  an 
angle  a  to  the  vertical,  and  the  potential  along  the  surface  of  the  sphere 
decreases  cosinusoidally  from  its  "equator"  (where  the  charge  is  maximum) 
toward  the  axis  of  polarization.  The  sphere  is  located  in  a  half-space  of 
uniform  resistivity.  This  is  a  rather  restrictive  model  but  is  one  of  the  few 
spherical  models  that  can  be  handled  analytically  with  a  simple  closed  form 
solution . 

152.  This  model  has  proved  useful  in  interpreting  SP  data  for  mineral 
deposits  that  have  similar  dimensions  along  all  three  axes.  For  seepage  prob¬ 
lems,  this  type  of  dipolar  sphere  might  in  some  cases  represent  flow  through  a 
roughly  spherical  cavity.  The  sphere  model  also  could  be  useful  for  initial 
interpretation  of  approximately  circular  field  anomalies,  to  check  whether  the 
model  curve  shapes  are  more  characteristic  of  a  point  source  (indicative  of  a 
large  ratio  of  source  depth  to  source  size)  or  a  spherical  source  (indicative 
of  shallower  burial  depth)  . 

153.  To  approximate  a  charge  distribution  other  than  the  sinusoidal 
dipole  of  this  model,  a  closely  spaced  distribution  of  point  sources  and/or 
sinks  may  be  placed  on  a  spherical  surface.  As  noted  in  some  of  the  refer¬ 
ences,  at  burial  depths  that  are  large  relative  to  the  radius  of  the  sphere, 
the  SP  field  of  this  spherical  model  approaches  that  of  a  simple  dipole  con¬ 
sisting  of  a  point  current  source  and  sink.  For  example,  the  SP  field  gener¬ 
ated  by  a  sphere  having  an  inclination  angle  a  of  30  deg  and  a  depth/ 
diameter  ratio  of  5  deviates  by  no  more  than  1.5  percent  from  the  field  of  a 
point  dipole  having  the  same  inclination  angle. 

Cylindrical  sources 

154.  The  source  geometry  and  modeling  equation  for  a  horizontal 
cylinder  of  infinite  strike  extent  are  shown  in  Plate  A4 .  The  cylinder 
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carries  a  uniform  dipolar  charge  around  its  circumference,  and  the  angle  of 
polarization  8  is  measured  from  the  vertical  axis.  The  cylinder  is  located 
in  a  half-space  of  uniform  resistivity.  As  for  the  sphere  above,  this 
restrictive  model  is  one  that  has  an  analytical  solution  in  closed  form. 

155.  The  relation  of  the  cylindrical  model  to  the  line  source  is  ana¬ 
logous  to  that  of  the  sphere  to  the  point  source.  Cylindrical  models  having 
other  than  constant  dipolar  charge  may  be  approximated  by  a  series  of  line 
sources  placed  around  a  cylindrical  circumference ,  and  line  sources  also  can 
be  used  to  approximate  a  cylinder  of  finite  length.  At  burial  depths  that  are 
large  relative  to  the  radius  of  the  cylinder,  the  SP  field  of  the  cylinder 
approaches  that  of  a  line  dipole  having  the  same  inclination  as  the  angle  of 
polarization  of  the  cylinder  and  a  separation  equal  to  the  diameter  of  the 
cylinder. 

Vertical  dipolar  sheet  source 

156.  Plate  A5  shows  the  geometry  and  modeling  equation  for  a  vertical 
rectangular  sheet  source  having  a  constant  positive  charge  per  unit  area  on 
one  face  and  an  equal  and  opposite  charge  on  the  other  face.  The  resistivity 
of  the  earth  on  the  two  sides  of  the  sheet  (in  the  x-  or  y-direction)  may  be 
different . 

157.  This  model  is  particularly  useful  because  it  has  been  observed 
that  ground-water  flow  in  the  vicinity  of  vertical  discontinuities  of  resist¬ 
ivity  and/or  electrokinetic  coupling  coefficient  often  generates  dipolar 
charge  distributions  of  this  type,  and  anomalies  fitting  this  model  have  been 
observed  in  a  number  of  field  studies.  Although  most  of  these  anomalies  were 
related  to  the  movement  of  geothermal  fluids  in  the  vicinity  of  fault  or 
fracture  zones,  similar  anomalies  also  have  been  observed  above  vertical  or 
nearly  vertical  geologic  features  in  the  vicinity  of  flows  of  nonthermal 
ground  witer. 

158.  More  complex  source  distributions  can  be  modeled  using  techniques 
discussed  by  Fitterman  (1979a, b,c)  or  by  approximating  the  sheet  with  a  dis¬ 
tribution  of  point  or  line  sources  and  sinks.  Sheets  that  are  not  vertical 
also  can  be  modeled  using  either  techniques  described  by  Fitterman  (1984)  or 
by  using  approximations  with  point  or  line  sources. 
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Quantitative  Modeling  Techniques 

159.  The  quantitative  SP  modeling  techniques  discussed  previously  are 
based  on  concepts  of  irreversible  thermodynamics  and  coupled  flows  of  fluids, 
heat,  electrical  current,  and  chemical  diffusion  as  described  by  Onsager 
(1931),  Pourbaix  (1949),  Denbigh  (1951),  Prigogine  (1955),  and  others.  Appli¬ 
cation  of  these  concepts  to  flow  in  soils  is  discussed  by  Mitchell  (1976). 

160.  Specific  application  of  these  concepts  to  interpretation  of  SP 

data  was  first  studied  by  Nourbehecht  (1963),  followed  by  the  work  of 
Fitterman,  Sill,  and  other  investigators  listed  in  the  previous  section.  Of 
particular  interest  are  publications  by  Sill  (1983a)  and  Sill  and  Killpack 
(1982).  The  first  of  these  summarizes  previous  work,  presents  a  number  of 
useful  type  curves,  and  shows  a  field  example  for  which  quantitative  tech¬ 
niques  were  used  to  interpret  SP  data  for  a  geothermal  area  in  terms  of  heat 

and  fluid  flow.  The  1982  publication  describes  a  computer  program  (SPXCPL) 
for  quantitative  two-dimensional  modeling  of  SP  data  generated  by  the  flow  of 
fluid  and/or  heat  in  the  earth. 

161.  Efforts  presently  are  underway  by  M.  Wilt  of  the  Engineering 

Geoscience  group  at  the  University  of  California,  Berkeley,  to  adapt  and 

document  SPXCPL  for  use  on  personal  computers  and  to  use  the  program  to  model 

SP  anomalies  generated  by  dam  seepage.  Plate  22  shows  results  of  a  prelimi¬ 
nary  run  of  SPXCPL  to  determine  the  general  SP  pattern  associated  with  dam 
seepage  flow.  A  relatively  negative  anomaly  is  seen  over  the  seepage  inflow 
area  and  a  relatively  positive  anomaly  is  seen  over  the  seepage  outflow  area. 
As  discussed  previously,  this  general  pattern  agrees  with  that  often  seen  in 
field  data. 

162.  A  specific  effort  presently  is  being  made  to  use  SPXCPL  to  help 
interpret  the  SP  data  from  Beaver  Lake  Dam,  Arkansas  (Part  III).  Results  of 
this  study  are  expected  to  be  presented  at  a  symposium  to  be  held  at  Karlsruhe 
University,  West  Germany,  in  the  spring  of  1988. 
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PART  V:  SUMMARY  OF  FINDINGS  AND  RECOMMENDED  PROCEDURES 


Summary 

163.  Conclusions  of  the  electrode  and  Beaver  Lake  Dam  studies  are 
included  in  Parts  II  and  III  of  this  report.  The  paragraphs  below  briefly 
summarize  findings  for  the  entire  research  effort  and  present  recommended 
technical  procedures  for  the  use  of  self-potential  methods  for  dam  seepage 
investigations . 

164.  Nonpolarizing  copper-copper  sulfate  electrodes  using  a  gelled 
electrolyte  appear  to  be  technically  suitable  for  monitoring  of  self¬ 
potential  (SP)  signals  generated  by  water  flow  associated  with  dam  seepage. 
Noise  levels  generated  by  sources  such  as  rainfall  and  temperature  variations 
are  predictable  and  are  considerably  lower  in  amplitude  than  expected  signal 
levels.  Less  expensive  electrodes  of  CCS  appear  to  have  signal  levels  com¬ 
parable  to  copper-copper  sulfate  electrodes,  but  considerably  higher  noise 
levels.  Electrodes  of  commercial-grade  metallic  lead  installed  at  depths 
comparable  to  the  copper-copper  sulfate  or  CCS  electrodes  do  not  appear  to  be 
suitable  for  long-term  SP  monitoring. 

165.  The  most  important  results  of  the  field  investigation  at  Beaver 
Lake  Dam,  Arkansas,  included  furnishing  of  data  for  the  electrode  studies 
described  above;  providing  information  about  the  relation  between  SP,  hydro- 
logic,  geologic,  and  other  geophysical  data;  establishment  of  baseline  data  on 
natural  and  artificial  noise  sources  at  a  typical  dam  site;  and  indicating  the 
need  for  electrical  resistivity  data  as  a  component  of  any  long-term  SP  moni¬ 
toring  program. 

166.  Comparison  of  SP  data  with  the  extensive  information  obtained  from 
other  investigations  at  this  site  indicated  that  significant  SP  anomalies  were 
associated  with  the  seepage  flow,  with  negative  anomalies  seen  above  downward 
or  horizontal  flow  and  positive  anomalies  above  areas  of  upward  flow.  Varia¬ 
tions  of  the  areal  SP  anomaly  pattern  between  high  and  low  pool  levels 
occurred  at  fault  or  fracture  zones  inferred  as  seepage  flow  paths  by  hydro- 
logic  or  other  geophysical  data,  and  these  SP  variations  indicate  significant 
spatial  variations  in  the  flow  pattern  between  high  and  low  pool,  i.e., 
changes  in  flow  paths.  SP  magnitude  variations  correlated  to  flow  rate  varia¬ 
tions  are  complicated  and  masked  by  the  changes  in  subsurface  resistivity. 
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167.  The  computer  program  developed  as  part  of  this  study  provides  a 
rapid  technique  for  estimating  the  depth  and  configuration  of  seepage  flow 
paths  determined  from  SP  data.  A  preliminary  study  using  a  recently  developed 
quantitative  modeling  program  indicated  that  quantitative  SP  modeling  methods 
based  on  concepts  of  irreversible  thermodynamics  and  coupled  flows  can  be  very 
useful  for  more  detailed  interpretation  of  seepage-related  SP  data. 

168.  The  bibliography  and  data  base  included  as  part  of  this  report 
summarizes  much  of  the  published  information  available  on  the  use  of  SP  tech¬ 
niques  for  seepage  and  other  flow-related  investigations.  This  data  base 
should  provide  a  useful  starting  point  for  future  studies  or  applications. 

Conclusions  and  Recommendations 

169.  Based  on  the  findings  of  this  research  effort  and  of  the  refer¬ 
ences  cited  in  the  bibliography,  the  self-potential  method  can  provide  useful, 
and  sometimes  unique,  information  about  dam  seepage  flow.  Selection  of 
appropriate  data  acquisition  and  interpretation  techniques  is  important  for 
effective  use  of  the  SP  method  for  this  application.  A  recommended  field  and 
data  reduction  procedure  for  simple  SP  surveys  is  presented  in  Appendix  B,  and 
some  recommended  techniques  for  SP  seepage  detection  and  monitoring  networks 
are  discussed  below. 

170.  The  recommended  method  of  acquiring  SP  data  for  seepage  monitoring 

is  the  use  of  a  permanently  installed  array  of  buried  copper-copper  sulfate 

electrodes  with  gelled  electrolyte.  All  the  electrodes  should  be  hard-wired 
to  a  common  station  having  weather  protection  and  power  available  to  run  a 
small  computer  system.  The  installation  also  should  include  permanently 
installed  arrays  of  resistivity  monitoring  electrodes  with  cables  running  to  a 
transmitter-receiver  unit  at  the  common  station.  Measurement  of  resistivity 
along  with  SP  is  necessary  if  SP  magnitudes  are  to  be  correlated  to  flow 
rates. 

171.  As  described  in  the  body  of  the  report,  appropriate  filtering  of 

the  SP  readings  is  required  to  minimize  the  effects  of  high-frequency  noise 

such  as  natural  tellurics  and  that  due  to  60-Hz  sources.  Each  SP  reading 
should  be  averaged  over  a  period  of  1  or  2  min  to  account  for  the  effects  of 
longer-period  telluric  variations.  Rainfall,  temperature,  pool  level. 
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seepage  flow  rates,  and  other  important  data  should  be  recorded  along  with  the 
SP  and  resistivity  values. 

172.  The  computer  system  would  be  used  to  acquire,  reduce,  store,  and 
plot  the  SP  and  resistivity  data.  The  entire  procedure  could  be  done  auto¬ 
matically,  on  a  daily  or  weekly  schedule.  Once  such  a  system  is  installed, 
the  only  labor  needed  would  be  periodic  collection  and  interpretation  of  the 
data.  If  desired,  the  data  could  be  transmitted  via  telephone  lines  to  a 
location  remote  from  the  installation.  This  type  of  computer  system  is  well 
within  the  present  state  of  the  art,  and  could  be  assembled  using  available 
personal  computers,  data  acquisition  systems,  and  other  commercial  components. 

173.  There  are  two  major  technical  advantages  of  such  a  system.  First 
is  the  ease  of  acquiring  and  processing  data,  which  allows  frequent  measure¬ 
ments,  reduced  the  probability  of  missing  a  significant  seepage  event,  and 
allows  developing  events  to  be  intensively  monitored.  Second,  the  data  qual¬ 
ity  obtained  from  a  permanently  installed  array  will  be  considerably  better 
than  that  obtained  from  repeated  conventional  field  surveys.  Even  using  the 
most  careful  field  procedures  as  described  in  Appendix  B,  a  certain  amount  of 
irreproducible  error  is  accumulated  each  time  an  electrode  is  put  into  the 
soil.  This  error  is  minimized  by  one-time  installation  of  a  permanent  array. 

174.  The  initial  cost  of  the  system  described  above  could  be  reduced  by 
employing  manual  rather  than  computer  acquisition  of  the  data.  However,  the 
cost  of  the  computer  system  is  small  relative  to  that  of  the  electrodes, 
cables,  and  installation  labor,  and  the  initial  savings  would  be  less  than  the 
increased  labor  costs  for  acquisition  and  reduction  of  data  (especially  if 
reading  frequency  is  increased  to  follow  a  developing  event). 

175.  Initial  costs  can  be  eliminated  entirely  by  employing  repeated 
conventional  field  surveys  rather  than  a  permanently  installed  array.  This 
would  require  mobilizing  and  transporting  personnel  and  equipment  to  the  field 
site  each  time  a  measurement  is  made.  As  discussed  above,  the  disadvantages 
of  such  an  arrangement  include  reduced  data  quality,  the  possibility  of  miss¬ 
ing  seepage  events  due  to  increased  sampling  intervals,  and  the  difficulty  of 
monitoring  developing  events.  Also,  because  successful  SP  field  data  acquisi¬ 
tion  requires  considerable  training  and  experience,  it  could  be  difficult  to 
maintain  data  continuity  and  quality  over  a  period  of  several  years. 

176.  Data  reduction  and  interpretation  techniques  would  be  similar  how¬ 
ever  the  data  were  acquired.  The  SP  data  would  be  referred  to  a  common  base 
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station  and  plots  along  each  profile  line  would  be  prepared.  The  entire  data 
set  would  be  contoured,  and  profiles  and  contours  of  differences  between  suc¬ 
cessive  data  sets  also  would  be  prepared. 

177.  Using  the  computer  program  developed  as  part  of  this  study,  along 
with  all  available  supplemental  hydrologic  and  geophysical  data,  the  profiles 
and  contours  would  be  interpreted  in  terms  of  seepage  flow  paths  and  changes 
in  these  paths  with  time.  If  more  intensive  interpretation  in  terms  of  seep¬ 
age  flow  rates  and  path  locations  is  warranted,  a  quantitative  interpretation 
program  such  as  that  described  in  the  body  of  this  report  could  be  used. 


PART  VI:  BIBLIOGRAPHY  AND  DATA  BASE 


178.  This  section  presents  a  bibliography  and  data  base  of  references 
related  to  the  acquisition  and  interpretation  of  seepage-related  SP  data. 
Topics  covered  include  theoretical  and  laboratory  studies  of  streaming  poten¬ 
tial  phenomena,  field  studies  of  seepage  problems,  studies  of  electrode  per¬ 
formance,  and  modeling  and  interpretation  of  SP  data.  A  list  of  key  letters 
for  retrieval  of  specific  topics  is  given  below.  Each  of  the  listed  refer¬ 
ences  includes  one  or  more  of  these  key  letters  indicating  the  main  topic (s) 
of  the  reference. 


_ List  of  Key  Letters _ 

B  Extensive  bibliography 

E  Electrode  studies 

DS  Dam  seepage  investigation 

F  In-field  seepage  investigation 

G  Geothermal  investigation 

L  Laboratory  measurement  of  streaming  potentials 

M  Modeling  or  interpretation  of  SP  data 

T  Streaming  potential  theory 

TM  Telluric  current  measurement 
0  Other 
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Note : 


Approximately  15s  offset  traces  due  to  recording 
pen  offset. 


High-Frequency  Noise  Records 


( 


Sample  Daily  Records  for  Copper  Sulfate 
and  Lead  Electrodes 


PLATE  3 


MODERATE 


PLATE  4 


LEAD  ELECTRODES 


330  360 


o 

CO 

rO 

o 

rO 

rO 

O 

o 

fO 

o 

f'* 

CM 

O 

00 

LJ 

CM 

Q 

O 

Q£ 

210 

ion 

O 

y 

2 

UJ 

o  "5 

LJ 

CO 

t— 

CO 

< 

c 

u. 

— 

J 

o  t- 
LiO 

Z5 

(/) 

< 

a: 

LJ 

^  w 

CL 

O  >s 

Q. 

CM  O 

O 

T-  Q 

a 

o  O 

>—  CM 

I  I 


(S}|OAI||ILU)  dS 


oiDOioocnocno 

^fOrOoitN'-'-dd 

(ssipui)  ||D;iroy 


PLATE  5 


Monthly  Monitoring  Records 


o 

'  • 

LI  T  C  M  OP 

:  r> 

/ 

X 

/ 

\ 

i  &  ■ n 

i..  ■  y 

/ 

f.  . 

y  $ 

1  * 

t 

i  / 

i 

l EGEND 

\ 

\ 

Survey  station  « 

i  / 

Contour  Interval  -  IOmV 

/ 

SOLE  IN  FEET 


[Sailfl 


FAULT  ZONE 


L'*  “t  '  £ 

O.fK.xC*'  _,»■**"  ' 


m  :  '\7\Ai  \as\  pi  «f  \  \ 

•  ~>.rv*v:*  V '  *  *  •*  \* 


27+00 


i  (-7 


^ ^  1  *7<:  V7 ; 

•  7r —  ■,  0tv  -|  „  i' 7  ,  , 

•i:|*  -  7-;  .  f  ■  \  •  /  c." 


,  r  -£» 

-  • 

C*>  -• 

.  •  C  Ax 

k  _  !*" 


ferj  \jg :. 


\  V 


■29+00 


V»  «  X* 


\  V  \ 


•+4aWi' 


\\V\  ■  \  \ 

\  \  \  \  \ 

1  \  ' 

!  17  , 


3i+oo 


33+00 


35+00 


68+00 


66+00 


64+00 


Self  Potential  Contours,  February  1987 


PLATE  « 


30 


123d)  NCLLVA313 


(Kill  ■’] 


Record  of  Self  Potential  Between 
Stakes  C42  and  C62 

Beaver  Dam.  Arkansas 


8/86  2/87 


230 

150 

41 . 6 

773 

L  2700 

136 

i 

i 


Interpreted  Resistivity  Layering 

Reaver  Dam.  Arkansas 


PLATE  1 


n  i  i  i  i  ;  i  i  i  i  i  i  )  u  i  i  i  m  i  i  i  i  i 

STATION  ( FT  > 


Field  Data  for  Line  C,  August, 1986 

leaver  .'an  •Arkansas 


i  i  I  I  :  >1  l  i  i 

8809 


PLATE  14 


Self-Potential  vs.  Elevation,  Line  D 

Beaver  Dam,  Arkansas 


PLATE  16 


PLATE  19 


Unsmoothed  Data,  Line  C,  Februar 

leaver  Dam.  Arkansas 


" £  SMOOTHED 


4  *2  '  #_  I'!  5 


,  *  s  # 
’«sl « =»; 


TTTTTTTTTTTTT 

STATION  (FT) 


I  i i i  1 1  i  ii  1 1  ii 

8090 


CTAVxro  purtATurr. 
j  i  n  %j  rivv  Alt  ui/ 


rTTTTTTTTTTTTTTTTTTTTTTTTTTTl 


8003 

STATION  <  'P'T  ) 


Smoothed  Data,  Line  C,  February  1987 

Beaver  Dam,  Arkansas  _ 


PLATE  20 


hj</>  I  W5 


COPPER  SULFATE  Li  NS  MOOT  H  FT* 


6809  8SG0 

STATION  <FT> 


Unsmoothed  Data,  Line  B,  February  1987 

_ Beaver  Dam,  Arkansas _ 

PLATE  21 


APPENDIX  A: 


PROGRAM  SP1.BAS:  INTERPRETATION  OF  SELF-POTENTIAL  DATA 
USING  GEOMETRIC  SOURCE  MODELS 


Program  Description 


1.  Program  SP1.BAS  is  intended  to  assist  in  the  modeling  and  interpre¬ 
tation  of  self-potential  (SP)  data.  The  program  is  written  in  BASIC  for  use 
on  IBM-PC  and  compatible  personal  computers.  Interpretation  is  done  by  com¬ 
parison  between  field  data  and  calculated  curves  generated  by  relatively 
simple  geometric  models. 

2.  The  models  include  charged  points,  lines,  spheres,  cylinders,  and 
sheets.  Model  geometry,  algorithms,  and  references  are  shown  in  Plates  A1 
through  A5,  and  details  of  individual  models  are  discussed  in  Part  IV  of  this 
report.  The  general  interpretation  procedure  is  as  follows: 

a.  Entry  of  field  data  (x  and  y  coordinates  of  measurement  sta¬ 
tions,  elevations,  and  measured  SP  values). 

b.  Correction  of  field  data  by  gradient  and/or  constant  shifts 
(if  desired) . 

c.  Selection  of  model  type  and  parameters  and  calculation  of 
anomaly  curve  (s). 

d.  Comparison  between  field  data  profiles  and  calculated  curves 
using  both  data  listings  and  graphic  displays. 

e.  Adjustment  of  model  parameters  and  reconparison  until  a  sat¬ 
isfactory  fit  is  achieved. 

3.  As  discussed  in  Part  IV,  preliminary  interpretation  procedures  using 
nomograms,  characteristic  points,  half-wavelength  comparison,  and  other 
techniques  are  given  in  the  literature.  Use  of  these  procedures  can  be  help¬ 
ful  in  estimating  initial  model  parameters. 

4.  Neither  the  program  nor  the  graphics  are  highly  refined,  but  because 
the  algorithms  are  relatively  simple  the  program  runs  quickly,  and  the  graph¬ 
ics  are  adequate  for  visual  curve  matching.  Although  the  program  has  been 

t' sted  extensively,  it  may  contain  some  "bugs."  The  authors  would  appreciate 
receiving  user  comments  regarding  problems  and  general  operation  of  the 
program. 

5.  The  program  is  interactively  driven  and  largely  self-explanatory. 

It  is  assumed  that  the  user  is  familiar  with  the  operation  of  the  IBM-PC, 
general  field  data  entry  procedures,  and  the  procedures  required  to  run  BASIC 
programs  on  the  PC.  Operating  instructions,  a  sample  run,  and  the  program 
listing  are  given  in  the  following  sections. 
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Operating  Instructions 


6.  Before  entering  BASIC,  the  "Caps  Lock"  key  must  be  "on"  (i.e., 
printing  capital  letters)  for  the  program  to  run.  Also,  if  it  is  desired  to 
print  graphic  output  on  a  dot  matrix  printer,  the  computer  must  be  in  "GRAPH¬ 
ICS"  mode. 

7.  As  discussed  previously,  the  general  operating  sequence  consists  of 
field  data  entry,  data  correction  (if  desired),  model  entry  and  curve  calcula¬ 
tion,  and  graphic  display  (if  desired).  A  field  data  set  always  must  be 
entered  first  to  create  station  coordinates  for  model  calculations.  If  only 
model  curves  are  desired,  a  "null"  data  file  with  zero  SP  values  at  each  sta¬ 
tion  may  be  entered.  The  program  must  be  rerun  to  enter  a  new  field  data  set. 

8.  The  field  data  may  be  entered  from  the  keyboard  or  from  a  pre¬ 
viously  generated  file.  After  keyboard  data  entry,  provision  is  made  for 
writing  the  data  to  a  disk  file  if  desired.  Calculated  model  output  data  also 
may  be  written  to  disk  files. 

9.  Up  to  14  model  curves  may  be  generated  for  a  given  field  data  set, 
but  only  four  of  the  curves  may  be  plotted  on  a  single  graph.  The  field  data 
profiles  may  be  omitted  from  the  plots  if  desired. 

10.  Additional  operating  information  is  listed  below.  Further  operat¬ 
ing  notes  are  included  as  typed  comments  in  the  sample  run  in  the  following 
section.  The  model  profiles  shown  in  this  sample  run  were  deliberately  chosen 
to  be  a  poor  "fit"  to  the  field  data  so  that  the  curves  could  be  distinguished 
in  the  plotted  graphs. 

a.  Selection  of  x-  and  y-axis  geographic  orientation  is  arbi¬ 
trary.  However,  as  the  cylinder  and  line  models  must  be 
parallel  to  the  y-axis,  this  axis  generally  would  be  oriented 
in  the  assumed  direction  of  flow.  The  sheet  model  may  be 
parallel  to  either  the  x-  or  the  y-axis. 

b.  An  output  file  name  always  must  be  entered  before  a  model  type 
is  selected.  Failure  to  enter  a  file  name  will  result  in  a 
fatal  error. 

c.  Data  may  be  plotted  along  either  the  x-  or  y-axis  of  the  input 
coordinate  system. 

d.  The  data  reduction  option  is  included  to  allow  linear  gradient 
corrections  and/or  constant  shifts  to  be  applied  to  the  field 
data.  Removal  of  gradients  and  constant  shifts  allows  the 
field  data  to  be  leveled  and  centered  about  the  zero  millivolt 
axis  for  better  fits  with  calculated  model  profiles.  This 
option  also  can  be  used  to  apply  tie-in  and  electrode  drift 
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corrections  to  field  data.  Use  of  the  option  is  illustrated 
in  the  sample  run. 

e.  The  point  source  model  includes  an  option  for  automatic  source 
input.  Inputs  include  the  number  of  points,  initial  (x,y,z) 
and  current  values,  and  desired  increments  for  these  param¬ 
eters.  This  option  allows  rapid  input  of  multipoint  models. 

f.  Multiple  sources  may  be  entered  for  all  models  except  the 
dipolar  sheet. 

g.  Model  source  current  or  voltage  may  be  varied  to  adjust  cal¬ 
culated  model  voltages  to  match  field  anomaly  amplitudes.  If 
it  is  desired  to  maintain  constant  source  current  or  voltage, 
the  "scale  factor"  included  in  the  model  input  may  be  varied 
instead. 

h.  After  model  parameters  are  entered,  there  will  be  a  pause 
after  "Strike  any  key  to  continue"  is  displayed  while  profiles 
are  calculated.  Calculation  time  depends  on  the  number  of 
stations,  the  number  of  sources,  the  type  of  source  model,  and 
computer  speed.  For  a  286-based  computer  with  a  math  copro¬ 
cessor,  calculation  of  the  models  shown  in  the  following  sec¬ 
tion  takes  a  few  seconds.  For  a  model  having  50  stations  and 
20  source  points,  calculation  on  the  same  machine  requires 
about  20  sec,  and  if  the  points  are  replaced  with  line  sources 
the  same  calculation  requires  about  30  sec. 


Sample  Run 

11.  This  section  contains  a  sample  run  of  SP1.BAS.  The  run  illustrates 
the  following  operations: 

a.  Keyboard  entry  of  simulated  field  data  and  creation  of  disk 
data  file  TEST.DAT . 

b^.  Reentry  of  field  data  from  file  TEST.DAT. 

c.  Correction  of  a  portion  of  the  field  data  using  the  data 
reduction  option  and  creation  of  corrected  data  file 
TESTC.DAT. 

d.  Entry  of  data  file  TESTC.DAT  into  the  program. 

£.  Input  of  models  for  point,  line,  sphere,  cylinder,  and  dipolar 
sheet  sources. 

_f.  Creation  of  a  plot  showing  the  field  data  and  calculated  pro¬ 
files  for  the  point  and  line  sources. 

g.  Creation  of  a  second  plot  showing  the  field  data  and  the  cal¬ 
culated  profiles  for  the  sphere,  cylinder,  and  dipolar  sheet 
sources. 

12.  As  mentioned  previously,  fits  between  the  simulated  field  data  and 
the  model  profiles  were  designed  to  be  poor  so  that  the  various  curves  would 
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be  distinguishable  on  the  plots.  A  few  more  trials  with  any  of  the  models 
would  have  resulted  in  much  better  fits  and  lower  RMS  errors. 

13.  The  material  on  the  following  pages  was  printed  directly  from  the 
screen  display.  Typewritten  comments  are  intended  to  help  clarify  some  points 
and  procedures. 
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SAMPLE  RUN 


SP  MAIN  MENU 

Choose  one  of  the  -following  actions! 
FI  -  Input  data. 

F2  -  Model  data  -select 
F3  -  Graph  data 
F 4  -  Reduce  data 


hit  the  space  bar  to  e>:  i  t  the  program. 


11  1ST  2RUN  3LOAD"  4 SAVE "  5C0NT6 , "LPT  1  7TR0NBT RUFFI.EY  OSCKEEN 

/ 

This  line  indicates  bottom  of  screen 
except  for  graphic  displays 

IF  YOUR  DATA  AF(E  ALREADY  IN  A  DATA  FILE  THEY  MUST  BE  IN 
THE  FOLLOWING  ORDER  IN  THE  FILE: 

STATION  NUMBER  ELEVATION  X  Y  V-MEASURED 
(X  or  Y  and  V- MEASURED  aust  always  b«  entered;  ELEVATION  and/or  STATION  may  be 
oaitted) 

IF  THIS  IS  NOT  THE  FORMAT  OF  YOUR  DATA  THEN  ENTER  FROM 
THE  KEYBOARD  OR  EXIT  AND  EDIT  THE  FILE 

DO  YOU  WISH  TO  INPUT  DATA  FROM  THE  KEYBOARD  OR  FROM  A  DATAFILE?  (t  UR  D) :  I 

WILL  YOU  BE  ENTERING  ELEVATIONS"'  (Y  or  N)  :  Y  aboard  entry  example 

X  s?  :  Y 

Y's?  :  Y 

V-MEASURED?  :  Y 

ENTER  UNITS  (F  for  Ft  —  M  for  meters  -  K  for  kms) :  F 


II.  1ST  2RUN  3LOAD"  4SAVE"  5C0NT  6 ,  "LF  T 1  7TR0NBT  RDFFKEY  OSCKEEN 
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KEYBOARD  ENTRY  SCREEN 

Enter 

a  9099  in 

X  or  Y  ooIubo  to 

quit 

TAT ION 

ELEVATION  X 

Y 

V-MEASURED 

FT 

FT 

FT 

mV 

1 

0 

-5 

o 

-20 

2 

0 

-4 

0 

—50 

3 

0 

-3 

0 

-70 

4 

0 

o 

-90 

5 

0 

-l 

o 

-95 

6 

o 

0 

o 

-50 

7 

o 

1 

0 

o 

B 

0 

o 

0 

30 

V 

0 

if 

0 

50 

10, 

0 

4 

0 

50 

1  l‘ 

0 

5 

0 

20 

12 

0 

9999 

to  save 

this  data 

set  enter 

a  -filename?:  TEST 

.  DAT 

1LIS1  2RUN  3L0AD"  4SAVE"  5C0NT b , "LFT 1  7TR0NBTR0FFKE Y 


OSCREEN 


SP  MAIN  MENU 


Choose  one  of 
FI  - 
F2  - 

F  4  - 


the  Following  actions: 
Input  data 
Model  data 
Graph  data 
Reduce  data 


Space  bar  was  strucK  to  exit  program  for  data  file 
entry  example.  Program  then  was  re-entered 
using  F2  Key  or  typing  RUN 

hit  the  space  bar  to  exit  the  program. 

1L1ST  2RUN  3L0AD"  4 SAVE"  5C0NT6 , "LPT  1  7TR0NBTR0FFLE Y  OSCREEN 
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IF  YOUR  DATA  ARE  ALREADY  IN  A  DATA  FILE  THEY  MUST  BE  IN 
THE  FOLLOWING  ORDER  IN  THE  FILE I 

STATION  NUMBER  ELEVATION  X  Y  V-MEASURED 
(STATION,  X,  and  V-MEASURED  must  always  be  entered;  ELEVATION  and/or  Y  may  hr 
omitted ) 

IF  THIS  IS  NOT  THE  FORMAT  OF  YOUR  DATA  THEN  ENTER  FRUM 
THE  KEYBOARD  OR  EXIT  AND  EDIT  THE  FILE 

DO  YOU  WISH  TO  INPUT  DATA  FROM  THE  KEYBOARD  DR  FROM  A  DATAFILE?  (I;  OR  D) :  D 

WILL  YOU  BE  ENTERING  ELEVATIONS'’  <Y  or  N)  « 

X  '  5?  I 

Y 's?  i 

VtMEASURED?  ; 

ENTER  UNITS  <F  Tor  Ft  -  M  -for  meters  -  K  Tor 


Y 

Y 

Y 

Y 

kms  >  i 


data  file  entry  example' 


1L1ST  2RUN  3L0AD"  4SAVE "  5CONT 6 ,  “LF'T 1  7 TRONBTROFFKEY  OSCREEN 
ENTER  DATAFILE  NAME!  TEST . DA  T 

ARE  THERE  STATION  NUMBERS  IN  YOUR  DATAFILE?  (Y  or  N)  :  ,Y 


1L  1ST 


2RUN  3L0AD" 


4SAVE"  5C0NT6 , "LET  1  7TRONGTROFFKEY 


OSCREEN 
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printout  of  input  data 
DATA  FILE!  TEST . DAT 


STATION 

ELEVATION 

X 

Y 

V-MEASURED 

FT 

FT 

FT 

mV 

1 

0 

—5 

0 

-20 

2 

0 

-4 

o 

-50 

3 

O 

-3 

0 

-70 

4 

0 

_2 

o 

-90 

5 

0 

-1 

0 

-95 

6 

O 

O 

0 

-50 

7 

o 

1 

o 

0 

B 

o 

2 

0 

30 

9 

o 

3 

0 

50 

10 

(3 

4 

0 

50 

1 1 

0 

5 

0 

20 

Strike  any  key  to  continue.  .  . 


1LIST  2RUN  3L0AD"  4 SAVE"  5C0NT 6 , "LPT 1  7TR0NBT ROFFKEY  OSCREEN 

Choose  one  of  the  fol  lowing  actions: 

FI  -  input  data 
data  reduction  F2  _  Mod.,  data 

exan'Ple  F3  -  Graph  data 

select - -  _  Reduce  dat» 


hit  the  space  bar  to  e>:  1 1  trie  pr  ogram. 

Enter  -first  station  number,  correction  ^  1,5  +5  correction  on  station  1 

Enter  last  station  number,  correction  ?  6,1  +1  correction  on  Station  6 

Enter  constant  correction  ?  B  constant  correction  +8  on  stations  1-6 

Correct  another  set  of  stations  <Y/N>?  7  Y 

Enter  -first  station  number,  correction  7  7,0 

Enter  last  station  number,  correction  7  11,-10 

Enter  constant  correction  7  -3 

Correct  another  set  of  stations  (Y/N)7  7  N 

If  you  wish  to  save  this  data  set  enter  a  filename:  1ES1C.DA7 
1LIST  2RIIN  3L0AD"  4SAVE"  SCON T6 ,  "LPT  1  7TR0NBTR0EFKEY  OSCREEN 

Exit  program  and  re-run  to  enter  data  file  TESTC.DAT 
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SF  MAIN  MENU 


Choose  one  of 


F3  - 
F4  - 


the  following  actions: 
Input  data 
Model  data 
Graph  data 
Reduce  data 


hit  the  space  bar  to  exit  the  program. 


1LIST  2RUN  3L0AD"  4SAVE ”  5CCJN1 6 , "LPT  1  7TRON8TRUFFKEY  OECREEN 


IF  YOUR  DATA  ARE  ALREADY  IN  A  DATA  FILE  THEY  MUST  BE  IN 
THE  FOLLOWING  ORDER  IN  THE  FILE: 

STATION  NUMBER  ELEVATION  X  Y  V -MEASURED 
(X  or  Y  and  V-MEASURED  Bust  always  ba  entered;  ELEVATION  and/or  STATION  nay  be 
omitted) 


IF  THIS  IS  NOT  THE  FORMAT  OF  YOUR  UA1A  THEN  ENTER  FRUM 
THE  KEYBOARD  OR  EXIT  AND  EDIT  THE  FILE 

DO  YOU  WISH  TO  INPUI  DATA  FROM  THE  KEYBOARD  OR  FF<OM  A  DATAFILE?  (I.  OR  D)  :  D 


WILL  YOU  BE  ENTERING  ELEVATIONS"’  (Y  or  N)  : 

X  '  s'?  : 

Vs?  : 

V-MEA5UF<E!)^  : 

ENTER  UNITS  <F  for  Ft  -  M  for  meters  -  K  for 


Y 

Y 

Y 

Y 

\  ms )  : 


F 


1LIST  2 RUN  3L0AD"  4 SAVE "  5C0NT6 , MLFT 1  7TRQN8T RUFFKE Y  05 GREEN 

ENTER  DATAFILE  NAME:  TESTC.DAT  input  corrected  data  file  TESTC.DAT 
ARE  THERE  STATION  NUMBERS  IN  YOUR  DATAFILE?  (Y  or  N):Y 


All 


DATA  FILE 


TEST  C. DAT 


STATION 

1 


4 

5 
b 
7 
B 


ELEVATION  X 
FT  FT 

O  -5 

O  -4 

O  -3 

O  -2 

O  -I 

O  O 

0  1 

0  2 


V  V-MEASURED 

FT  mV 

0  -7 

0  -37 . BOOOOOO 1 1 97097 

O  -5B. 60000002 384 1 86 

0  -79. 40000009536743 

O  -85. 20000004769372 

O  -4  1 

O  —  3 

O  24.5 


9  O 

1 0  O 

1  1  O 


4 

5 


O 

O 


o 


4.2 

39.  5 


7 


Strike*  any  key  to  continue.  .  .  . 


1LIST  2RUN  3LUAD"  4 SAVE"  5CUNT6 , ”LP T 1  7TR0N8  TROFFl  EY 


OSCEEEN 


Begin  input  of  models 

SF  MAIN  MENU 


Choose 

select  - 


one  of 
F  1  - 
—  F2  - 
F3  - 
F  4  - 


the  following  actions: 
Input  data 
Model  data 
Graph  data 
Reduce  data 


lLisr 


hit  the  space  bar  to  e:<  1 1  the  program. 

2RUN  3L0AD"  4SAVE"  5CGN76,"LP71  77 R0NB7 RUT P LEV 


OSCRLLN 
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Point  source  entry 

Entnr  an  output  filenames  POINT. OUT  an  output  file  name  must  be 
select  source  model  types  entered  for  each  model 

ENTER  1  FOR  SPHERICAL  SOURCE 

ENTER  2  FOR  POINT  SOURCE - select  Notes  enter  numbers  1-5, 

ENIER  3  FOR  HORIZONTAL  LINE  SOURCE  not  keys  F1-F5 

ENTER  4  FUR  DIPOLAR  SHEET  SOURCE 
ENTER  5  FOR  HORIZONTAL  CYLINDER  SOURCE 

?  2 

Enter  resistivity,  no.  of  source  points  (<2Q0) ,  scale  factor 

?  100,2,1 

SOURCE  COORDINATE  ENTRY:  MANUAL  <1)  OR  AUTO  (2)^  2 
ENTER  INITIAL  X  SOURCE,  Y  SOURCE,  Z  SOURCE,  CURRENT 
?  0,0, 2, -15 

ENTER  INCREMENTS  FOR  X  SOURCE,  Y  SOURCE,  Z  SOURCE,  AND  CURRENT 
?  3,0,0,35 

SOURCE  SOURCE  Y  SOURCE  Z  SOURCE  CURRENT 

1  0  0  2  -15 

2  3  O  2  20 

Strike  any  key  to  continue 

AUTO  source  entry  used  here  for  illustrationi 
usually  used  only  for  entry  of  large  number  of  sources 


JL1ST 


2RUN  3L0AD"  4SAVE"  5C0NT6 , “LPT  1  7T RONBTROFFKE Y 


OSCREEN 


MODEL  NUMBER  = 

STATION  X 


Y  V-MEASURED  V- CALCULATED  7.  ERROR 


1 

-5.0 

0 . 0 

2 

-4.0 

o .  o 

- 

T 

-3.0 

0 . 0 

- 

4 

-2.0 

0 .  0 

- 

5 

-1.0 

o .  0 

- 

6 

C> .  0 

0 .  0 

- 

7 

1.0 

0 . 0 

8 

2 . 0 

0.  0 

9 

3 . 0 

0 . 0 

10 

4 . 0 

0.0 

1  1 

5.0 

c 

o 

-  7 . 0 

-5.7 

-18.  1 

37.8 

-9.7 

-74. 4 

58.6 

-15.9 

-72.9 

79.4 

-25. 3 

-68.  1 

85.2 

-35.6 

-58.  2 

4  1 . 0 

-31 . 1 

-24.2 

— 3 . 0 

5.8 

-292.5 

24.5 

57.9 

1  36 . 5 

42.0 

92.9 

121.3 

39.5 

89.0 

1  25.  2 

7.0 

68.  2 

874.4 

RMS  ERROR  =  289.0 

Strike  any  key  to  continue. 

1LIST  2RUN  3L0AD"  4  SAVE  ’’  5C0NT  6 ,  "LPT  1  7T  RONBT  RUEF IEY  OSCREEN 
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SP  MAIN  MENU 


Choose  one  of  the  following  actions: 
FI  -  Input  data 

select  -  F2  _  Modl?i  data 

F3  -  Graph  data 
F4  -  Reduce  data 


hit  the  space  bar  to  exit  the  program. 


1L1ST  2RUN  3LOAD"  4SAVE "  5C0NT6,"LPT1  77 RONBTRQFFkEY  OSCREEN 

Line  source  entry 

Enter  an  output  filename:  LINE. OUT 
SELECT  SOURCE  MODEL  TYPE: 

ENTER  1  FOR  SPHERICAL  SOURCE 
ENTER  2  FOR  POINT  SOURCE 

ENTER  3  FOR  HORIZONTAL  LINE  SOURCE - select 

ENIER  4  TOR  DIPOLAR  SHEET  SOURCE 
ENTER  5  FOR  HORIZONTAL  CYLINDER  SOURCE 

LINE  SOURCE (S)  PARALLEL  TO  Y-AXIS 

Enter  resistivity,  number  of  source  lines  (  <20  ),  scale  factor 

?  1 00 ,2,1 

Enter  source  line  center  point  X,  Y,  depth,  current,  line  half-length 
2  0,0, 2, -2, 10 
7  _< ,  O ,  ,7 , 1 0 

SOURCE  X  SOURCE  Y  SOURCE  Z  SOURCE  CURRENT  LENGTH 

1  0  0  2  -2  10 

2  7  0  2  2  10 

Strike  any  key  to  continue 


1LIST  2RUN  3L0AD"  4 SAVE"  5C0NT6 , "LPT  1  7T RONBTROF FEE Y  OSCREEN 
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MODEL  NUMBER 


2 

STATION  X  Y 


1 

-5.0 

0.  0 

2 

-4.0 

0.0 

-if 

-3.0 

0.0 

4 

-2.0 

0.0 

5 

-1.0 

0.0 

6 

0.0 

0.0 

7 

1.0 

0.0 

e 

2.  0 

0.0 

9 

3.0 

0.0 

10 

4.0 

o.  o 

l  i 

5.0 

0.0 

RMS  ERROR  =  191.4 

Strike  any  key  to  continue.  .  . 

1L1ST  2RUN  3L0AD"  4SAVE"  5C0NT6,“LP 


V-MEASlJRED  V- 

■CALCULATED 

7.  ERROR 

-7.0 

-22.  5 

221 . ; 

-37.8 

-26.9 

-28.9 

-58.6 

-3-2.2 

-45.  1 

-79.4 

-38.0 

-52.  1 

-85.  2 

-41.9 

—  50 . 8 

-41.0 

-36.2 

-11.8 

-3.0 

-14.5 

303.5 

24.5 

14.5 

-40.  8 

42.0 

36.2 

-13.9 

39.5 

41.9 

6.2 

7.0 

38 . 0 

443.6 

1  7TR0NBTR0FFI.E  Y 

OSCREEN 

SP  MAIN  MENU 


Choose  one  of 
FI  - 

select - F2  - 

FT.  - 
F4  - 


the  Following  actions: 
Input  data 
Mode)  data 
Graph  data 
Reduce  data 


hit  the  space  bar  to  exit  the  program. 


1LIST 


2RUN  TLOAD 


4SAVE "  5C0NT 6  ,  "LF'T  1  7TRON0T ROFF kE Y 


OSCREEN 
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Sphere  source  entry 
Enter  an  output  filename:  SPHERE. OUT 

SELECT  SOURCE  MODEL  TYPE: 

ENTER  1  FOR  SPHERICAL  SOURCE  -  select 

ENTER  2  FOR  POINT  SOURCE 
ENTER  3  FOR  HORIZONTAL  LINE  SOURCE 
ENTER  4  FOR  DIPOLAR  SHEET  SOURCE 
ENTER  5  FOR  HORIZONTAL  CYLINDER  SOURCE 

?  1 


Enter  number  of  spheres 

?  I 

ENTER  SPHERE  CENTER  POINT 
?  0,0,4,2,120,500 

SOURCE  X  Y  Z 

1  0  O  4 

Strike  any  key  to  continue 


<  <  200 ) 

<X,Y,Z>,  RADIUS,  INCLINATION,  AND  VOLTAGE 


RADIUS  I  NCI.  I  NAT  I  ON 
2  ^-120 


VOL  T  AGE 
500 


degrees 


1LIST  2RUN  3L0AD"  4 SAVE "  5C0NT6 , "LPT  1  7TRDNBTR0FF LEY 


OSCREEN 


MODEL  NUMBER  = 

STATION  X 


Y  V-MEASURED  V-CALCULATLD  7.  ERROR 


1 

-5.0  0.0 

-7.0 

-48.  2 

588. 9 

2 

-4.0  0.0 

-37.8 

—60.  4 

59.7 

3 

-3.0  0.0 

-58.6 

-73.6 

25 .  5 

4 

-2.0  0 . 0 

-79.  4 

-83.  5 

5.  I 

5 

- 1 . 0  0 .  o 

-B5.  2 

-B1 .8 

-  4 . 0 

6 

0 .  O  0 . 0 

-4  1 . 0 

-62.5 

52.  4 

7 

1 . 0  0 .  0 

-3.  i.i 

-32. 4 

978.5 

8 

2 . 0  0 .  O 

74.5 

-6.0 

-174.5 

9 

3 .  O  0 .  0 

42.  0 

9.6 

-77. 2 

J  0 

4 . 0  0 .  O 

39 . 5 

16.7 

-59.0 

1  1 

5.0  0.0 

7 .  O 

17.8 

153.6 

RMS  ERROR  =  351.6 

Strike  any  key  to  continue.  .  . 

1L1ST  2RUN  3L0AU"  4 SAVE"  5CONT6 , "LPT  1 

7T  RONBT ROFFt.EY 

OSCREEN 

A16 


rJi  m 


SP  MAIN  MENU 


Choose  one  of  the  following  actions! 
FI  -  Input  data 

select  - F2  -  Model  data 

F3  -  Graph  data 
F4  -  Reduce  data 


hit  the  space  bar  to  exit  the  program. 


1L1ST  2RUN  3LOAD"  4SAVE"  5C0NT6 , "LF'T t  7TR0NBTR0FELEY 

Cylinder  source  entry 
Enter  an  output  filename:  CYL.OUT 
SELECT  SOURCE  MODEL  TYRE: 

ENTER  1  FOR  SPHERICAL  SOURCE 
ENTER  2  FOR  POINT  SOURCE 
ENTER  3  FOR  HORIZONTAL  LINE  SOURCE 
ENTER  4  FOR  DIPOLAR  SHEET  SOURCE 

ENTER  5  FUR  HORIZONTAL  CYLINDER  SOURCE -  select 

7  t* 

Cylinders  of  infinite  strike  extent  parallel  to  y-asiis 

Enter  resi st i vi ty , number  of  cylinders  <<20>,  scale  factor 
1  00 ,1,1 

Enter  cylinder  axis  x-distance,  radius,  inclination,  depth, 


e  _  t  e  —  *>ri  *7  i 

SOURCE  X 

SOURCE  RADIUS 

1NCL I NAT  ION 

DEPTH  CURRENT 

1 

-1  2 

-3u 

2  1 

Strike  any  ley 

to  continue 

degrees 

OSCREEN 


current 


1LIST 


2RUN  3L0AD”  4SAVE 


5C0NT  6 , "LPT  1  7TRDNBT ROFF \E Y 


OSCREEN 


CO  * 


MODEL  NUMBER  =  4 

STATION  X 


V  V-MEASURED  V-CALCULATED 


I 

-5.0 

0.0 

-7.0 

-28.  3 

2 

-4.0 

0.0 

-37.  8 

-34.3 

3 

-3.0 

0.0 

-58.  6 

-40.8 

4 

-2.0 

o.  o 

-79.  4 

-41.3 

5 

-1.0 

0.  o 

-85.2 

-21 . 8 

6 

o 

o 

0.0 

-4  1 . 0 

6.5 

7 

1 . 0 

0.  0 

-3.  0 

1  9 . 0 

8 

2.0 

0.  o 

24.5 

20.9 

9 

3.0 

0.0 

42.0 

19.6 

10 

4.0 

0.0 

39.  5 

17.6 

1  1 

5.0 

c 

c 

7.0 

15.8 

MS  ERROR  =  24  7. 8 

trite  any  t ey  to  continue. 
1L1SI  2RUN  31.0AD"  4SAVE " 


5CONT 6  ,  "LPT  1  7TRCJN8T  R'EJFFI  EY 


OSCREEN 


SP  MAIN  MENU 


Choose  one  of 


select 


Ft  - 
F2  - 
F7  - 
F4  - 


the  <ol lowing  actions: 
Input  data 
Model  data 
Graph  data 
Reduce  data 


hit  the  space  bar  to  emt  the  program. 

1LIST  2r<UN  3LQAD "  4SAVE"  5C0NT 6  ,  "LF'T  1  7TR0NB1  ROT F I  E  Y 


(.'SCREEN 


7.  ERROR 


3o4.  I 
-9.2 
-30 .  4 
-47. 9 
-74.4 
-113,. 9 
-734.2 
-14.  fc 


1 


A18 


tj  U  L* 


Sheet  source  entry 

Enter  an  output  filename:  SHEET. OUT 
SELECT  SOURCE  MODEL  TYPE: 

ENTER  1  FOR  SPHERICAL  SOURCE 

ENTER  2  FOR  POINT  SOURCE 

ENTER  3  FOR  HORIZONTAL  LINE  SOURCE 

ENTER  4  FOR  DIPOLAR  SHEET  SOURCE  - Select 

ENTER  5  FOR  HOR I ZONTAL  CYLINDER  SOURCE 

?  4 

ENTER  SHEET  INPUT  PARAMETERS: 

ENTER  SHEET  TOP  DEPTH,  BOTTOM  DEPTH,  STRIKE  LENGTH  1,7,10 
ENTER  RH01  (Y  OR  X  NEG. ) ,  RH02  (Y  OR  X  PUS. ) ,  SCALE  FACTOR  10,5,250 
ENTER  SHEET  CENTER  COORDINATES  <X,Y>  1,0 

IS  SHEET  PARALLEL  TO  X  OR  Y  AXIS?  (ENTER  X  DR  Y >  Y 

AXIS  DTOP  DBOT  LENGTH  RH01  RH02  SCALE  XS  YS 

Y  1  7  10  10  5  250  1  O 

DO  YOU  WISH  TO  ENTER  A  NEW  MODEL  WITH  THE  SAME  STATIONS?  <1=Y,2=N0>2 


1LIST  2RUN  3L0AD "  4 SAVE "  5C0NT6 , "LPT  1  7TR0N8TR0FFKE Y  OSCREEN 


MODEL  NUMBER  =  5 

STATION  X  Y  V-MEASURED  V-CALCULATED  7.  ERROR 


l 

-5.  0 

0.0 

-7.0 

-42.7 

509.5 

-4.0 

0 .  0 

-77.8 

-50.  7 

77 . 11 

■J 

-  3 . 0 

<  ;> .  o 

-58.6 

-58.8 

0. 4 

4 

-2.0 

o .  i  ;> 

-79.4 

-67.2 

-15.4 

5 

-1.0 

o .  o 

-85.  2 

-  7  1 . 6 

-15.9 

6 

0.0 

o .  o 

-4  1 . 0 

-59 .  6 

45.5 

7 

1 . 0 

0 .  0 

-  7 . 0 

0 . 0 

-  1 00 .  1 1 

8 

2.0 

o .  o 

24.5 

29.0 

21.7 

9 

3.0 

o .  t  ;> 

43.0 

75.0 

-14.7 

10 

4.0 

o .  0 

39.5 

77.6 

-15.0 

1  1 

5.0 

0.  o 

7.0 

29.  4 

32o .  2 

RMS  ERROR  =  105.1 

Etrile  any  ley  to  continue.  .  . 

1LIST  2RUN  3L0AD"  4SAVE"  SC0NT6,"LPT1  7TR0NBT ROFFLE Y  OSCREEN 

End  input  of  models 


A19 


O  L* 


Begin  graphic  display  option 

SP  MAIN  MENU 

Choose  one  of  the  -following  actions: 

FI  -  Input  data 
F2  -  Model  data 

select  -  F3  -  Graph  data 

F4  -  Reduce  data 


hit  the  space  bar  to  exit  the  program. 


1L1ST  2RUN  3L0AD"  4SAVE "  5CONT 6 , "LPT  1  7TRDNBT ROFFKEY  OSCREEN 

PLOT  NEW  MODELS  ON  PREVIOUS  AXES  <Y/N>?  N 

ENTER  INITIAL  AND  FINAL  STATIONS  TO  BE  PLOTTED  1,11 

PI  OT  ALONG  X  OR  Y  AXIS  (ENTER  X  OR  Y)  X 

ENTER  XM I N , XMAX  FOR  GRAPH  X-AXIS:  -5,5 

ENTER  YM I N , YMAX  FOR  GRAPH  Y-AXIS:  -100,100 

ENTER  X-AXIS  GRID  LINE  INTERVAL:  1 

ENTER  Y-AXIS  GRID  LINE  INTERVAL:  20 

ENTER  X-AXIS  TICK  MARK  INTERVAL:  .5 

ENIER  Y-AXIS  TICK  MARK  INTERVAL:  10 

ENTER  A  GRAPH  TITLE:  POINT  AND  LINE  SOURCES 

ENTER  AN  X-AXIS  TITLE:  STATION 

ENTER  A  Y-AXIS  TITLE:  SP 

Choose  up  to  Pour  (4)  data  sets,  TOTAL.,  to  graph. 

GRAPH  THE  MEASURED  (FIELD)  DATA"?  (Y  OR  N)  :  ?  Y 

HON  MANY  CALCULATED  MODEL  CURVES  WOULD  YOU  LIIE  TO  GRAPH?  (0-4):  2 

FOR  CALCULATED  DATA  SET  1 

ENTER  MODEL  #:  1  ,  t 

plot  point  sources  (model  1),  line  sources 

cnrv  rfll  ri  li  n-rrrn  nnm  cct  o  (model  2)  and  field  data 


FOR  CALCULATED  DATA  SET  2 
ENIER  MODEL  #:  2 


A20 


FUR  CALCULATED  DATA  SET:  1 

WHICH  SYMBOL  WOULD  YOU  LIKE  : 

12  3  4 

Square  Dot  Triangle  Circle  (1,2,3,  or  4):  1 
DO  YOU  WISH  TU  CONNECT  THE  POINTS  (Y  OR  N)?t  Y 
WHAT  COLOR  LINE?  (1-3)  1 

FOR  CALCULATED  DATA  SET :  2 

WHICH  SYMBOL  WOULD  YOU  LIKE  I 
12  3  4 

Square  Dot  Triangle  Circle  (1,2,3,  or  4>:  3 

DO  YOU  WISH  TO  CONNECT  THE  POINTS  (Y  OR  N>7:  Y 
WHAT  COLOR  LINE?  (1-3)  2 

WHICH  SYMBOL  WOULD  YOU  LIKE  FOR  THE  MEASURED  DATA: 
12  3  4 

Square  Dot’  Triangle  Circle  (1,2,3,  or  4):  4 
DO  YOU  WISH  TO  CONNECT  THE  POINTS  (Y  OR  N)?i  Y 
WHAT  COLOR  LINE?  (1-3)  3 


POINT  AND  LINE  SOURCES 


StriKe  space  bar  to  clear  display 
and  return  to  main  menu 
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SP  MAIN  MENU 


Choose 

select  - 


one  of 
FI  - 
F2  - 
—  F3  - 
F4  - 


the  foil  owing 
Input  data 
Model  data 
Graph  data 
Reduce  data 


act i ons: 


hit  the  space  bar  to  exit  the  program. 


1LISI  2RUN  3L0AD"  4SAVE"  5C0NT 6 , "LF'T 1  7TR0NBTR0FFKEY  OSCREEN 


illustrates  use  of  previous  axes 

PLOT  NEW  MODELS  ON  PREVIOUS  AXES  <Y/N>?  Y  note  that  previous  title 

Choose  up  to  four  <4>  data  sets,  TOTAL,  to  graph.  a^so  is  used 
GRAPH  THE  MEASURED  (FIELD)  DATA”  <Y  OR  N)s  7  Y 

HOW  MANY  CALCULATED  MODEL  CURVES  WOULD  YOU  LIKE  TO  GRAPH?  (0-4) s  3 
FOR  CALCULATED  DATA  SET  1 
ENTER  MODEL  #:  3 

plot  sphere  (model  3),  cylinder  (model  4), 
F0F<  CALCULATED  DATA  SET  2  sheet  (model  5),  and  field  data 

ENTER  MODEL  #:  4 

FOR  CALCULATED  DATA  SET  3 
ENTER  MODEL  #:  5 


FOR  CALCULATED  DATA  SET:  1 

WHICH  SYMBOL  WOULD  YOU  LIKE  i 
12  3  4 

Square  Dot  Triangle  Circle  <1,2,3,  or  4>t  1 

DO  YOU  WISH  TO  CONNECT  THE  POINTS  (Y  OR  N)?s  Y 
WHAT  COLOR  LINE?  (1-3)  1 
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FUR  CALCULATED  DATA  SET  t  2 
WHICH  SYMBOL  WOULD  YOU  LIKE  j 
12  3  4 

Square  Dot  Triangle  Circle  <1,2,3,  or  4):  2 

DO  YOU  WISH  TO  CONNECT  THE  POINTS  <Y  OR  N)?l  Y 
WHAT  COLOR  LINE?  <l-3>  2 
FOR  CALCULATED  DATA  SET:  3 

WHICH  SYMBOL  WOULD  YOU  LIKE  : 

12  3  4 

Square  Dot  Triangle  Circle  <1,2,3,  or  4):  3 

DO  YOU  WISH  TO  CONNECT  THE  POINTS  <Y  OR  N>7:  Y 
WHAT  COLOR  LINE^  < 1 —3)  3 

WHICH  SYMBOL  WOULD  YOU  LIKE  FOR  THE  MEASURED  DATA: 
1  2-3  4 

Square  Dot  Triangle  Circle  <1,2,3,  or  4>:  4 

DO  YOU  WISH  TO  CONNECT  THE  POINTS  <Y  OR  N>7:  N 


note  incorrect  title  from  previous  plot 

POINT  AND  LINE  SOURCES 


100 

S 

P 


100 
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JU. 


Re-plot  entering  new  axes 

BP  MAIN  MENU 


Choose  one  of  the  following  actions: 
FI  -  Input  data 
F2  -  Model  data 

select  - F  3  -  Graph  data 

F4  -  Reduce  data 


hit  the  space  bar  to  e:<it  the  program. 

2RUN  3LCJAD" 


1LIST 


4SAVE 


SCONTfc,  "LF'Tl  7TRUNB1 ROFFKEY 


OSCREEN 


PLOT  NEW  MODELS  ON  PREVIOUS  AXES  (Y/N)?  N 

EN1ER  INITIAL  AND  FINAL  STATIONS  10  BE  PLOTTED  1,11 

PLOT  ALONG  X  OR  Y  AXIS  (ENTER  X  OR  Y )  X 

ENTER  XMI N , XMAX  FOR  GRAPH  X-AXIS:  -5,5 

ENTER  YM IN, YMAX  FOR  GRAFH  Y-AXIS:  -100,100 

EN1ER  X-AXIS  GRID  LINE  INTERVAL:  1 

ENTER  Y-AXIS  GRID  LINE  INTERVAL:  20 

ENTER  X-AXIS  TICK  MARK  INTERVAL:  .5 

ENTER  Y-AXIS  TICK  MARK  INTERVAL:  lO 

ENTER  A  GRAPH  TITLE:  MODELS  3,  4,  AND  5  commas  not  allowed  in  titles 
7Redo  from  start 

ENTER  A  GRAPH  TITLE:  MODELS  3-5 
ENTER  AN  X-AXIS  TITLE:  STATION  (FT) 

ENTER  A  Y-AXIS  TITLE:  mV 

Choose  up  to  four  (4)  data  sets,  TOTAL,  to  graph. 

GRAPH  THE  MEASURED  (FIELD)  DATA?  <Y  OR.  N):  ?  Y 

HOW  MANY  CALCULATED  MODEL  CURVES  WOULD  YOU  LIKE  TO  GRAPH?  (0-4):  3 

FOR  CALCULATED  DATA  SET  1 
ENTER  MODEL  #:  3 


plot  sphere  (model  3),  cylinder 
(model  4),  sheet  (model  5)  and 
field  data 


FOR  CALCULATED  DATA  SET  2 
ENTER  MODEL  4:  4 

FOR  CALCULATED  DATA  SET  3 
ENTER  MODEL  «l  5 

FOR  CALCULATED  DATA  SET:  1 

WHICH  SYMBOL  WOULD  \ OU  LIKE  : 

12  3  4 

Square  Dot  Trianqle  Circle  (1,2,3,  or  4):  1 

DO  YOU  WISH  TO  CONNECT  THE  POINTS  (Y  OR  N>?:  Y 
WHAT  COLOR  LINE"’  (1-3)  1 

FOR  CALCULATED  DATA  SET:  2 

WHICH  SYMBOL  WOULD  YOU  LIKE  : 

12  3  4 

Square  Dot  Trianqle  Circle  (1,2,3,  or  4):  2 

DO  YOU  WISH  TO  CONNECT  THE  POINTS  (Y  OR  N>?:  Y 
WHAT  COLOR  LINE-?  (1-3)  2 
FOR  CALCULATED  DATA  SET:  3 

WHICH  SYMBOL  WOULD  YOU  LIKE  I 
12  3  4 

Square  Dot  Triangle  Circle  (1,2,3,  or  4):  3 
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MODELS  3 


5 


Program  Listing 

14.  This  section  contains  the  complete  listing  for  SP1.BAS.  The 
program  is  written  in  BASIC  and  requires  about  17500  bytes  of  storage.  Data 
input  and  output  files  require  a  few  hundred  to  a  few  thousand  bytes,  and 
machine  memory -requirements  are  minimal. 
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JO  •  PROGRAM  SP1.BAS 
20  '  RFC/CM  JANUARY  19B8 
30  ' 

35  PI  -  3. 14159 

40  '  MAIN  MENU  FOR  SP  PROGRAM 

50  ' 

52  SCREEN  2: SCREEN  Oi CLOSE# 1 

60  ON  KEY  (1)  SOSUB  lOOO 

70  KEY  (1)  ON 

80  ON  KEY  (2)  GOSUB  2000 

90  KEY  (2)  ON 

lOO  ON  KEY  (3)  GOSUB  3000 

110  KEY  (3'  CN 

111  ON  KEY  (4)  GOSUB  7000 

112  KEY  <4 )  ON 

113  CLS 

120  LOCATE  3,35«PRINT  "SP  MAIN  MENU" 

130  LOCATE  5,25lPRINT  "Choose  one  of  the  following  actions:  " 


140 

LOCATE 

6 , 34  :  PR  I  NT 

"FI  - 

I  nput 

data" 

150 

LOCATE 

7  , 34 1 PRINT 

"F2  - 

Model 

data" 

160 

LOCATE 

8 , 34 1  PRINT 

"F3  - 

Graph 

data" 

162 

LOCATE 

9 , 34 : PRINT 

"F4  - 

Reduce  data" 

170 

LOCATE 

22, 10:PRINT 

"hi  t 

the  space  bar 

180  IF  INKEYS< >CHRS (32)  THEN  180 
200  END 
1000 
1010 
1020  • 

1030  '  DATA  ENTRY  SECTON 
1040  • 

1050  initialization 
1060  ■ 

1070  DEFDBL  V,X 

1080  DIM  ST AT ( 1 00) , ELEV (100) , X ( 100) , Y ( 1 00) , V ( 1 5 , 1 00 ) , I SYM (5) ,2 (loo) 

1081  DIM  LR(20> ,H(20>  ,CURR(20) 

1082  DIM  RAD  ( 200) , ALF (200) , E ( 200) ,  ALFR(200),  BETA  (20) 

1085  DIM  XG(IOO) .YG(IOO) 

1087  DIM  XS ( 200)  , YS (200 )  ,  2S  ( 200)  , CUR ( 200) 

1090  DIM  LCOLR ( 5 )  ,C ( 15)  ,RHO( 15)  ; XO ( I  5 )  , YO ( 1 5 >  ,20115)  ,ERRO (100)  , PERERR (100) 

1091  INCHECK-1 
1100  CLS 

1110  K=»  1  :  KK—  1 


1120  PRINT  "IF  YOUR  DATA  ARE  ALREADY  IN  A  DATA  FILE  THEY  MUST  BE  IN  " 

1130  PRINT  "THE  FOLLOWING  ORDER  IN  THE  FILE:  " : PR  I  NT 

1140  PRINT  "STATION  NUMBER  ELEVATION  X  Y  V-MEASURED " : PR  I  NT  "  (X  or  Y  and 
V-MEASURED  must  always  be  entered:  ELEVATION  and/or  STATION  may  be  omitted) 

1150  PRINT:  PRINT  "IF  THIS  IS  NOT  THE  FORMAT  OF  YOUR  DATA  THEN  ENTER  FROM  " : PR  I N 
T  "THE  KEYBOARD  OR  EXIT  AND  EDIT  THE  FILE" 

1160  PRINT  " 

1170  INPUT  "DO  YOU  WISH  TO  INPUT  DATA  FROM  THE  KEYBOARD  OR  FROM  A  DATAFILE'"  (>. 
OR  D):  " , IN# : PR  I  NT 

1180  • 

1190  variables  used  and  units  checks 
1200 

1210  INPUT  "WILL  YOU  BE  ENTERING  ELEVATIONS’?  (Y  or  N) : 

1220  INPUT  "  Xs?  : 

1230  INPUT  "  Y  s? 

1240  INPUT  "  V-MEASURED?  : 

1250  INPUT  "ENTER  UNITS  (F  for  Ft  -  M  for  meters  -  E  for 


"  ,EI  * 

",** 

"  ,Y* 

"  ,  V» 

kits):  ”,  UN  1  I  * 
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1260  IF  UNIT*=*"F“  THEN  LABL*=*"FT"  ELSE  IF  UNIT*="M"  THEN  LABL*=“M"  ELSE  LABL*="K 
M" 

1270  CLS 

1280  IF  INS«"D"  THEN  1680 
1290  ' 

1300  LOCATE  2,30:PRINT  "KEYBOARD  ENTRY  SCREEN- 

1310  LOCATE  3 , 25: PRINT  "Enter  a  9999  in  X  column  to  quit" 

1320  ' 

1330  LOCATE  4,1 

1331  PRINT  TAB (121  “STATION"! 

1332  PRINT  TAB <241  "ELEVATION" ; 

1333  PRINT  T AB ( 36 )  "X“» 

1334  PRINT  TAB ( 48 >  “Y“| 

1335  PRINT  TAB (58)  "V-MEASURED" 

1337  PR I NT -TAB (24)  LASL*; 

1338  PRINT  TAB (36)  LABL#: 

1339  PRINT  T AB ( 48 )  LABL* ( 

1340  PRINT  TAB (58)  "mV" 

1350  ' 

1360  LOCATE  (K+61,1 
1370  PRINT  TAB (12)  KK; 

1380 

1390  IF  EL*<>"Y"  THEN  1420 

1400  LOCATE  (K+6) , 24: INPUT  “ " , ELEV (KK ) 

1410  ’ 

1420  IF  X*0"Y"  THEN  1460 

1430  LOCATE  (K+6) , 36 : I NPUT  "" ,X(KK> 

1440  IF  X (KK) “9999  THEN  1560 
1450  ' 

1460  IF  Y*<>"Y"  THEN  1500 

1470  LOCATE  (K+6) , 48: INPUT  , V ( Kk ) 

14eO  IF  X*="N"  AND  Y(KK>*=9999  THEN  1550 
1490  • 

1500  IF  V*<>"Y"  THEN  1530 

1510  LOCATE  (K+6) ,58: INPUT  "",V(OtKK> 

1520  ' 

1530  IF  K= 15  THEN  K=O.FOR  1=6  TO  22:L0CATE  1,1 sPRINT  " 

"  :  NE  X  T  I 

1531  K=k+liKK«KK+l 
1540  GOTO  1360 
1550  ' 

1560  NUMPTS=»KK-1 

1561  ’ 

1562  datafile  save 
1570 

1580  INPUT  "If  you  wish  to  save  this  data  set  enter  a  filename:  ",DATFIL* 

1590  IF  DATF IL#=" "  THEN  1980 
1600  OPEN  DATF I L*  FOR  OUTPUT  AS  #1 
1610  FOR  L-l  TO  NUMPTS 

1620  IF  X*="Y"  AND  Y*="Y"  AND  EL*»"Y"  AND  V*=”Y"  THEN  DT A*=STR* (ELEV ( L ) ) + "  " + 

STR*(X  (L>  )+■"  "+STR*  (Y(L)  )  +"  "+STR*  ( V  (O  ,L)  ) 

1630  IF  Y*«'“N"  AND  EL*«"Y"  THEN  DTA*=BTR*  (ELEV  (L  )  )  +  "  "  +  STR* ( X  (L >  > ♦  "  "+STR*(V 

<0,L)  ) 

1640  IF  EL* = "N"  AND  Y*="N"  THEN  DT A*=STR* ( X (L ) ) ♦ "  " +STR* (V (O ,L ) 1 
1645  IF  EL*«"N"  AND  X*="N"  THEN  DTA*=STR* (Y (L) ) +"  " + STR* ( V (0 , L ) ) 

1650  PRINT41  ,L; DTA* 

1660  NEXT  L 
1670  GOTO  1980 
1680  ‘ 

1690  '  datafile  read 
1700  ' 

1710  INPUT  "ENTER  DATAFILE  NAME:  " , DATF IL* 
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1  720 
1  7  _'.0 
1  7  40 
1  7  SO 
1751 
1  752 
1 753 
1  754 
1755 

1757 

1758 

1759 
1  760 

17  BO 
1790 
1  BOO 
1810 
1 320 
IB  30 
1040 
1 B50 
I860 

18  70 
1  330 
1890 
1  900 

19ol 
1  9(.i.? 
19o7 
19o4 
1  9i.)5 
1  9  1  0 
1  97o 
1971 
1922 
1937 

19  2  4 
1970 
1940 
1  950 
1960 
19  71.) 
1900 
1981 
1990 
7oOO 

200  i 
2020 
7070 
7040 
704  1 
2013 
704  4 
2016 
7'.)  17 
7040 

204  9 
705 1 
2052 

205 


OPEN  DATFILi  F C1K  INPUT  AS  #1 

INPUT  "APE  THERE  STATION  NUMBERS  IN  YCJUR  DAIAFILE^  (Y  or  N):",S1 
CLS-.LOCA1E  2 , 20  :  PFt  I  NT  "DATA  FILE:  " .DATEILJ 
LUCA  IE  4,5 

PRINT  T  AR (17)  "SIATION": 

PRINT  I AR ( 7  4 )  "LI  EVA  1  ION"  ; 

PRINT  TAP (76)  "X"; 

PRINT  TAP (40)  "Y"; 

PRINT  TAB (SB)  " V- MEASURED" 

PRINT  TAB  (24)  LABL.T; 

PRINT  TAB  (76)  LABI,  t ; 

PRINI  TAB (48)  LABL1 ; 

PRINT  TAB  <  5B )  "mV" 


FOR  1=1  TO  lOoO 
IF  EOF  < 1 )  THEN  1920 

IP  St~  "-Y  "  THEN  INPUTH1  ,STAT  (  I  ):  LOCATE  <  I  I +6  >  ,  1  1  :  PRINT  SlAl(I): 

IF  EL1  =  "Y"  THEN  I NPUT # 1  , ELEV ( I )  : LOCAT E  ( I  I  +  6 )  , 23: PR I  NT  ELEV ( I ) 

IF  X  t-=  "  Y  "  THEN  INPUT # 1 , X ( I  1 : LOCATE  <  1 1 +  6 )  ,  35:  PR  I N  I  X(I); 

IP  Y  T  =  " Y "  THEN  INPUT #1 , Y ( 1 ): L0LA1 E  ( I  1 +6 )  , 4 7 : F R I  NT  Y(I> 

IP  VT="Y"  THEN  I  NPUT  #  1  ,  V  <  O ,  I  )  :  LOC  A  I F.  < 1 1 + 6 > , 57 : r R I N I  V(n,I) 

IP  11  HUD  15-0  AND  11  14  THIN  19(-7  ELSE  19o5 

F'RIIII  "Sir  1 1  p  tiny  I  t;y  to  continue.  .  .  " 

IP  INI  EY  *=■""  THEN  1907 
I.ULAIE  6,1:11=0 
11=11+1 
NE  X  I  I 


NIJMF  15=1-1 

Bl  AN!  =NUMPTS  MUU  1 5:  F  I LL=  1 5- BL ANI  +1 
FUR  1=1  TO  FILL 
T PINT " 

NEXT  I 


LOCATE  22.1 

PRINT  "Strilf?  riny  ley  to  continue.  .  .  . 

IF  I  N1  E  Y  f  "  THEN  1960 


CLOSE# 1 
SOI  0  52 
RETURN 

'in)*****###***# 

'  MUDEL  INPUT  SECTION 


CLS 

IP  INCHED  =0  mr;N  2042  1.1  5E  3"47 

FOR  1  =  1  TO  200: L0CA1 E  22,1:PRINT  "NO  DA  1  A  EXISTS  YE T ■  • " : N1  X T  1:0070 


MOI)F.L  =  MUI)Fl  +  1 

INF-’IJI  "Enter  mn  output  ti  lename:  ",UUTF1H 

OPEN  OUTFILt  FUR  OUT PUT  AB  #1 
PRINI  "SELECT  SOURCE  MUDEL  TYPE:" 

PRINT  "  ENTER  1  FOR  SPHERICAL  SOURCE" 

F'RINT  "  ENIER  2  FUR  POINT  SOURCE" 
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2054 

FRINT  " 

EN  1  ER 

3  FOR 

HORIZONTAL  LINE  SOURCE1 

2055 

PRINT  " 

ENTER 

4  FOR 

DIPOLAR 

SHEET  SOURCE" 

2057 

FRINT  " 

ENTER 

5  FOR 

HORIZONTAL  CYLINDER  SOURCE" 

2061 

INPUT  STYE 

2062 

IE  ST iP  = 

1 

THEN 

GOSIIB 

2346 

2063 

IE  STYE  = 

o 

THEN 

GOSUR 

2400 

2064 

IE  STYE’  = 

3 

THEN 

GOSUFi 

2600 

2065 

IE  ST  YP  = 

4 

THEN 

GOSIIB 

2  BOO 

2066 

IE  STYE  = 

5 

THEN 

GOSUB 

2500 

21  00 

Ffi  I  N  MM  , 

" 

MODEL  NUMBER: 

**  ,  MUDf  L 

2110 

FR I N 1 #  1  ,  " 

" 

2 1  20 

FRINMM  , 

"STATION 

X 

V -MEASURED 

O-CALCULAT  ED 

V.  EKROfV  M 

2135 

FOR  J  =  1 

TO  NONETS 

2136 

IF  V<0,J) 

O  THEN  DERR  =  .1 

ELSE  DERR  =  0(0, 

,  J  ) 

2140  ERREKJ  )  =  <V  (MODEL. ,  J  >  -V<«>.  J  )  »  /  DERR  : FERERR (J>  =ERR0  I  1  >  »  t"0 

2171  F'RINntl  .USING  "««#««  ####**.#  #«»«««.«  FMMMMHI.W  «««««« .  H" ;  J  ,  X  ( J 

>  ,  V  ( O .  J  )  ,  V  (  MODEL  ,  J  )  ,  FERERR  ( J  ) 

2 1  f JO  NFXT  .1 

2 1 05  SUMS  -  0 

2170  FOR  1=1  TO  NLIMF  IS 

22<  '0  SUMS  =  SUMS  +  FERERR  IF  )  2 

2220  ME XT  1 

22T.ii  RMS-  (SUMS/MUMF1S)  .5 

22  TO  PRINT  HI  ,  "KMS  ERRUR=  “  ;  RMS 
2250  ’ 

2260  ELS: LOCATE  3,1 

2261  PRINT  "  MODEL  NUMBER  =  ".MODEL 

2263  FRINT 

22/0  FRINT  "  STATION  X  V  V-MEASURED  O-CALCUIAIF 

D  7.  ERROR  " 

22RO  LOCATE  10,1 
2270  FOR  L= 1  TO  NUMT 1 S 

73' >0  FRINT, USING  "#«##«  ««««##.»  «#«#*«.«  #«««##. 8  HtUHtHtt.H  IMHMH* 

tt.  H  "  iL  ,  X  (L>  ,  Y  <L  )  ,7(0,1)  ,  V  (MODEL  ,L  >  .FERERR  (L  ) 

2301  II  1  MUD  13  =0  THEN  23-02  ELSE  2"  I" 

2302  FRINT  "St.rile  .inv  I  pv  to  continue.  .  .  " 

73-' >3  H  T N1  EY  J  =  "  "  THEN  2303 

2  Tut  LOCATE  10,1 
2  T 1 '  >  NEXT  L 

2311  FI  AN!  =MUMF’TS  MOD  13  :  F  1  LL=  1  3 - RL  ANI  +  1 
2T1  2  FOR  1-1  10  F  Il.l. 

7313  FRINT" 

231 T  ME  XI  I 

7  T  1  5  FRINT  "RMS  ERROR  =  ”  :  :  f  Tv  I  NT  UR  I  NIT "  H  8#  8  88  8  8  .  8  "  ;  RMS 

2316  FRINT  "Stri  I  e  ,my  I  !?y  to  continue.  .  . 

2317  IF  INI  EYJ  =  ""  T 1  It  N  231/ 

2  7'.( ) 

2340  LI  USE  It  1 

2341  GO lO  52 

2342  RETURN 

2346  Model  input  For  sinole  point  source 

2  350  ’  Begin  sphere  model  input 

2352  PRINT  "  Enter  number  oF  spheres  <  2001” 

235)4  INPUT  NSniJE 

2356  FRINT  "ENTER  STRIKE  CENIER  F'UINI  (X,Y,Z>,  RADIUS,  INCLINATION,  AND  UUL 1  AGE " 
2350  FOR  1  =  I  TO  NSOUR 

2360  INEU!  XS(I),  YS ( I ) , ZS ( I ) ,RAD < I ) , ALF ( I ) ,E ( I ) 

2361  ALFR(I)  =  ALF ( I ) /57, 2950 

2362  NEXT  I 
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£364  PRINT  "  SOURCF  X  Y  Z  RADIUS  INCLINATION  VOLTAGE" 

2 '•■66  FOR  I  =  1  TO  NSOUR 

27-60  PR  INI  “  "  ;Is"  ";XStI>j"  “;YS(I);"  ":ZS<I);"  " :  T.nD  (  I  )  ; 

"  5  Al  F  <  I  )  ;  "  ";E(I) 

2770  NEXT  I 

2772  PRINT  "Strile  any  ley  to  continue": 

2774  IF  IMLEYf =" "  1  HEN  2774 
2778  PRINT 

2780  FOR  J  =  1  TO  NUMPTS 
23R2  SV  =  <j 

2784  FOR  I  =  1  TO  NSOUR 

2386  NUMER  =  ( X  <  J  )  — XS  ( 1 )  )  *S  I N  ( Al.FR  (  I  )  )  +  ZS  (  I  )  « COS  ( AIT  R  (  I  )  > 

2388  DEN  =  <  < X < J ) -XS < II  )  '2  +  ( Y ( J ) -YS < I ) ) ' 2  +  <ZS(I>>'2  )  1.5 

2390  VFART  =  E ( I ) * (RAD ( I ) " 2) *NUMER/DEN 
2792  SV  =  SV  +  VFART 

2794  NEXT  1 

2795  V(NODEL..,J)  =  SV 

2796  NEXT  J 

2797  RETURN 

2400  Multiple  point  sources  model 
2404  F=I 
2406  C4  =  i;> 

2403  Pegin  model  parameter  input 

2410  ERINI  "  Enter  resistivity,  no.  of  source  points  (  2o0),  stale  factor 
2112  INPUT  RT  It)  ( MODEL  )  ,  NSOUR.  T 

2414  INPUT  "SOURCE  COORDINATE  EN1RY:  MANUAL  (1)  liR  All  I O  <2>":C7 
24  16  IF  C  7  =  1  THEN  2430 
24  18  IF  C 7  -  '•  THEN  24 '-H 

2419  Manual  source  coordinate  entry 

2420  FR INI  "ENIER  SOURCE  X,  Y,  Z,  AND  CURRENT" 

2422  FOR  I  =  1  TO  NSOUR 

7424  INPUT  XSI1),  YS<I),  ZS(I),  CUR ( I ) 

2  4  26  Nr:  XI  I 

7478  ERINI  "  SOURCE  X  SOURCE  Y  SOURCE  2  SOURCE  CURRENI" 

2470  FOR  I  ir  v  TO  NSOUR 

24’ ?  PRINT  "  "5  I:"  " ; XS ( I ) ; "  ";>5<I>;”  ":7S(I):" 

: CUR < I ) 

2474  NEXT  I 

2475  PRINT  "Strile  ari>  ley  to  continue"; 

2476  IF  INIEtT=”"  THEN  2436 

2477  GOTO  2470 

2478  Auto  source  coordinate  entry 

2440  PRINT  "ENTER  INITIAL  X  SOURCE,  Y  SOURCE ,  Z  SOURCE,  CURREtll " 

7442  INPUT  XS1 ,  YS 1 ,  7S1,  CUR1 

2444  PRINT  "ENTER  INCREMENTS  FOR  X  SOURCE,  V  SOURCE,  Z  SOURCE,  AND  CURRENI 
2446  INPUT  DXS,  DYS ,  DZS,  DCUR 
7448  FUR  I  -  1  TO  NSOUR 


245') 

XS  <  I  ) 

=  XS1  r 

( I  - n  »nxs 

2452 

YS  I  1  ) 

=  YS1  + 

(1-1)  «I)YS 

2454 

ZS(  I  ) 

=  ZS1  + 

< I-l ) »DZS 

7456  CIJR(I)  =  CUI\  1  +  (  1  -  1  )  »DCUR 

2458  NEXT  I 

2460  PRIM  I  "  SOI  11(1  E  X  SOURCE  Y  SUURCE  Z  SOURCE  CURRENI" 

247,7  FOR  I  =  1  TO  NSOUR 

247.4  PRINT  "  "?Ii"  "  ;  X  S  (  I  )  s  "  "  :YS(I>;"  ";7SII>:  " 

: cur  < i ) 

2466  NEXT  I 

2467  P'R  I  NT  "Strile  any  ley  to  c  on  1 1  n  ue  "  ; 

2468  IE  INI  E  Yf  •-  "  "  THEN  2468 
2470  FOR  J  -  1  TO  NUMPTS 
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f)  M  M  U 


2472  S V  =  O 

2474  FOR  I  =  1  TO  NSDUR 

24  76  R  =  SOR< <X (J)-XS(l) >  2  + (Y < J > -YS < I > >  2  -»  (  Z  ( J  ) -25  <  I  >  > '  2) 

24  78  VF  ART  =  (F *CUR < II *RHO ( MUDEL M / < 2«P I *R > 

2 4 BO  SV  *■  SV  +  VF  ART 

24B2  NEXT  I 

2404  V (MODEL, J)  =  SV 

2485  NEXT  J 

2494  RETURN 

2500  Cylinder  model  input 
2502  F  -  1 

2504  F'RINT  "Cylinders  of  infinite  strife  extent  parallel  to  y-a.:  i  s  "  :  F  F<  IN  1 
250fc  PRINT  "Enter  r es i it l v i t y , number  of  cylinders  (  20 ) ,  scale  factor 
2508  I  NT  II I  RHO(MODEI.  )  ,NCYLS,F 

25 1 0  PRINT  "Enter  cylinder  a::i5  x-distance,  radius,  inclination,  depth,  current" 
2512  FOR  1  =  1  TO  NCYI.S 

2514  INM.lt  XSU  )  ,RAD  (  1  )  ,L<ET0<  I)  ,H<  I  )  ,CURR<  I) 

25  Id  NEXT  I 

2518  PRINT  "  SOURCE  X  SOURCE  RADIUS  INCLINATION  DEPTH  L'URRf  111" 

2519  FUR  I  =  1  TO  NL'YLS 

2520  PRINT  "  "sis"  ” ;XS(1)5"  " ; RAD ( 1 )  : "  " : HE  1  A ( 1  >  ;  "  ";ll 

(f);"  ” ; CURR ( I ) 

2522  NEXT  I 

2524  F'RINT  "Slrile  any  ley  to  continue"; 

2576  IF  1  NT  EY  t  =  " "  THEN  2526 
2528  FUR  J  =  1  1U  NUMF I S 
2570  SV  =  0 

2572  rUR  1  =  1  111  NLYLS 

2574  RETAR  =  RE T A ( 1) /57. 2958 

2576  FI  =  2«CURR  (  I)  *RHU  (MODrt  )  •RAD  <  1  ) /F  I 

2578  NUN!  R  ~  (  X  ( J  )  -  XS  (  I  )  ULUS  (  EF  T  ART  +  H  <  I  )  »S  I M  ( 1<F.  I  HR  > 

2540  DEMON  =  (X(J)-XS(l))  2  *■  Hill  2 

2547  VPART  =  F *F  1  * NUI1E R/ DENOM 
254  4  F.V  -  SV  +  VPART 

2546  NEXT  I 

2548  V ( MUDEL , J )  =  SV 
2550  NEXT  J 

7552  RETURN 

26o0  Line  source  model 
2("  i 2  F  -  1 

2ti>7  TRIM  "l  INC  S'  HIM  t  (5)  THEM. ML  TU  Y  -  A  X  15":  Tl  INI 

26<'5  PRINI  "fnt.er  resistivity,  number  of  source  lines  1  7o  ),  scale  facior  " 
261"  INI  111  RHU<  MUD!  L  >  ,  Nl  INI'S,  ! 

2615  F'RINT  "Enter  sour  te  line  center  point  X,  Y,  depth,  current,  line  ha  1  f -  1 mp  t 
h" 

2620  F  UR  1  e  1  TCI  III.  INF  5 

2675  1 1 11  T!  1  XS  (  I  )  ,  YS  (  I  )  ,H  <  1  1  ,  CURR  (  I  )  ,LR  l  l  ) 

267"  HE  XI  1 

26~5  PRINT  •'  SOURCE  X  SOURCE  <  SOURCE  7  SUURT  E  CURRFN1  I  LN5  1H  • 

264"  FUR  1  =  1  TU  NL I  MRS 

2645  PRINT  "  1  •  "  ";XS(I):"  "  ;  X  S  (  1  )  ;  ”  "  ;  H  M  I  ;  "  ";C"I. 

R  (I  >  5  "  " 5  LR ( I ) 

265't  NEXT  I 

2655  F  k  1  NT  "St  rile  any  ley  to  continue"; 

766"  IF  INI  EYT  =  " "  THI  N  766" 

7665  ri)R  0  =  1  TO  NUMF' T  S 
670  SV  =  O 

675  FDR  1  =  1  TO  Nl  INFS 

68"  FI  =  CURR  (  1  )  «RHF)  (MUDEL)  /  (2»P  I  ) 

685  YF'L  =  Y  ( J  ) -YS  (  I  )  +LR  <  I  > 


A32 


M  n  tj  *j  -*j  fj 


2690  YML  -  Y<J) -YS( I ) -LR( I > 

2695  NtltiFR  =  YFL  +  SPR  (YFl.  2  +  (X(J)-XS(I))  2  +  M  <  I  )  2) 

27o0  DENTJM  -  YML  +  SL'R  ( YML  2  +  (X(J)--XS(I))  2  +  H(l)  2) 

2705  Or  AM  =  ri  «  l  UO  (NUMER/DENOM)  «  F 

27  tO  SO  =  SV  +■  VP  ART 
2715  NEXT  I 

2720  V< MODEL, J  >  =SV 
2  7  25  NE X  T  J 
2750  RE  T  I.IRN 
2800  PRINT 

2802  ’  Fitterman  Short  Modnl 

2805  PRINT  "ENTER  SHEET  INFUT  PARAMETERS:" 

2808  FRINT 

2010  INFUT  "ENTER  SUET  T  TOP  DEPTH,  POT  TOM  DEPTH,  STRIIE  LEMUIH  "  ,  A ,  T' ,  l 

2815  INFUT  "ENTER  RHU1  (Y  OR  X  NEG.l,  RH02  <Y  UR  X  PUS. > ,  SCALE  T AC T UR  "  ,R1,R2, 

NF 

2820  INFUT  "ENTER  SHEET  CENTER  COORDINATES  <X,Y>  ",XS,YS 

2822  INFUT  "IS  SHEET  FARALLEL  TO  X  UR  Y  AXIS-1  TENTER  X  OR  Y)",FAXi 

2825  FLINT 

2824  FRINT  "  AXIS  DT  OF  DBOT  LENGTH  RH01  Rl  102  SCALE  XK 

YS“ 

2825  PRINT  "  "sTAXT:"  ";Aj"  ";B;"  ";Li"  "jRls"  ":R2:" 

" : NF ;  "  " :  X  S ;  "  "  i  YS 

20  '2  NIUT  -  NUMFTS 

2877  IF  FAXi="X"  THEN  2875  ELSE  2920 

28  75  F OR  I  -  1  TO  NTLIl 
2840  YD  =  Y  < 1 ) - YS 
2845  XU  -  X ( I ) -X5 

050  IF  YD  =  O  THEN  2855  ELSE  2865 
855  V  O-IODEL  ,  IT  *0 
860  GOTO  2915 

865  IF  VD  0  THEN  FE  =  NF  .  ( F  I  «  (  1  <  FV2  /R 1  >  > 

970  IF  YD  o  THEN  FE  -  NF / (F  1 » t 1 +R l / R2 ) ) 

8  75  XF  =  X  (  I  )  -XS-U.  72 
880  XM  =  X ( t ) - XS-L/2 

2085  FA  =  AIN  (  CD* XF) / T YD»5UR l XT  2  +  YD  2  * L  7>  )  ) 

2090  FD  -  AIN  (  ( A*  Xl  >  /  ( YD*SOR  (  XT  2  *  Y I)  2  +A  2>  )  ) 

7895  re  -  ATN  (  <D*XM) / <YD*SOR(XM  2  *  YP  2  U'  7  >  )  I 

29oo  FD  -  ATN  (  (A*XM)  /  (  v  D  *  SDR  <  X  M  2  *  YD  2  +A  2)  )  ) 

'79o5  V  (MODEL,  I  )  =  FE  »  (FA  —  F  E<  -FL  *  F  l> ) 

2915  NEXT  I 

2920  FUR  I  ~  1  TO  NIUT 
2972  YD  =  Y  < I ) -YS 
7974  XD  -  X ( I ) -XS 

2926  IF  XD  -  O  THEN  2928  ELSE  2972 
7928  V ( MU DEI  ,  I  )  =0 
2970  GOTO  2957 

79  72  IF  XD  0  THEN  FE  =  NT  /  (  T  I  *  (  1  *  R7/ R  1  )  ) 

2974  IF  XD  O  THEN  FE  =  NF /(Fl»(l+Rl/R2)> 

79  76  Y F-  =  Y  (  I  >  -  Y S *  L  /  7 

29  73  YM  -  Y  <  I  > -YS--L/2 

294o  FA  =  ATN  (  <  R#YT)  /  (  XD«5PR  (  YF'  7*  XD  2*1'  71)' 

29  4  7  F  Ei  =  A  T  N  (  ( A*  YF )  /  ( X  D»SOR  <  Yf '  2tXl>  2*A  7)1) 

2944  FL'  =  ATN(  (R»YM)  /  <XD<SUR<YM  2+XD  7*  Ft  2)1) 

2946  FD  n  A  TN  (  ( A*  YM  )  /  (  XD*  SDR  ( YM  2+XL)  2+A  2))) 

794(3  V  (MODEL  ,  I  1  =  FE  *  (F  A-f  B- F  C  +  F  D ) 

2752  NEXT  I 

2905  INFUT  "DU  YUU  WISH  TO  ENTER  A  NEW  MODEL  WITH  HIE  SAME  E !  A  T  I UNV.  '  (1  =  Y,7‘N(0" 
,C.F2 

2990  IF  CF2  =  1  THEN  ;  (TOO  ELSE  2995 


A3  3 


2995  RETURN 

300 1  '  PLOTTING  ROUTINES 

3015  CI  S:  I.EV  OFF 

3016  INPUT  "PLOT  NEW  MODELS  ON  PREVIOUS  AXES  <Y/N>-  "  ,  REF  L  T  f 

3017  IF  REFLIt  -  "Y"  GOTO  3070  ELSE  GOTO  3020 
3021.1  ’  CM  9/ 28/ 86 

3022  INPUT  "ENTER  INITIAL  AND  FINAL  STATIONS  TO  RE  PLOTTED  ",STPl,SfP2 

3024  NT  L  T  -■  SPP2  -  STPl  +  1  * 

3026  INPUT  "FLOT  ALONG  X  OR  Y  AXIS  (ENTER  X  OR  Y)  " , T  AX I 

3028  IF  FAX*  =  "X"  THEN  3030  ELSE  3036 

3030  FOR  I  =  1  TO  NFLT 

3032  XG ( I )  =  XISTPl  +1-1) 

3034  NEXT  I 

3036  FOR  I~1  TO  NFLT 
30*8  YG(I)  =  Y(STP1  +1-1) 

3040  NEXT  I  • 


3047 

INPUT 

"FITTER 

XMIN.XMAX  FOR  GRAPH  X-AXIS: 

" , XM1N, 

XMnX 

304  4 

INPUT 

"EN1ER 

YMIN, YMAX  FOR  GRAPH  Y-AXIS: 

" , YMIN, 

YMAX 

304  6 

INPUT 

"ENIER 

X-AXIS  GRID  LINE 

INTERVAL : 

" , XGIN1 

3043 

INPUT 

"ENTER 

Y-AXIS  GRID  LINE 

INTERVAL: 

" . YGINT 

I  NPl  i  r 

"ENIER 

X-AXIS  T  I  Cl  MART 

INTERVAL: 

",  XT  INI 

303  3 

INPUT 

"ENIER 

Y-AXIS  TICI  MARI 

INTERVAL: 

"  , YT I  NT 

3«  >S4 

I  Ml  7  1 1 

"EN I ER 

A  GRAP'H  T11LE: 

" , T  ITLE  1 

3  oVj6 

INPUT 

"EN PER 

AN  X-AXIS  TITLE: 

",  X  T  I  1 1. 1 

INF  LI  I 

"EN TER 

A  Y-AXIS  TITLE: 

".YTITI.I 

Zu  7<.» 

PRINT 

"Choose 

i*(i  to  Pour  (4) 

data  eets. 

IOTA!  ,  to 

cjr  aph  .  " 

31 080 

INPUT 

"GRAF  H 

THE  MEASURED  (TIILD)  DA  1  A-1 

(Y  OR  N>  : 

";f  i 

3,1  y?i\ 

I  NP'U  I 

"HOW  MANY  CALCULATED  MODEL  CURVES 

WOULD  YOU 

LIKE  TO  GRAIH 

, MODLN 

3100  IF  MODLN  =  0  THEN  3210 

7110  FUR  1  +  1  TO  MODLN :  PR  I N  T  "F0(<  CALCULATED  DATA  SET  "j  1:  INPUT  "ENIEE  MODEL  « 
" ,MODL < I > : PR  INI  "  " : NE X  T  1 

7120  FOR  ]=■*]  TO  MODI  N 

7170  PRINT  'TOR  CALCULATED  DATA  SET:  ";1 

7140  PRINT  "WHICH  SYMBOL  WOULD  YOU  LItE  : 

7150  PRINT  "  1  2  7  4" 

7  1  A*1  I  NT  UT  "  Squar  e  Dot  1  r  i  ana  1  P  C  i  r  c  1  p  (1,2,7,  or  4  )  :  "  ,  1  B  XII  (  I  ) 

7170  INPUT  "DO  YOU  WISH  TO  CONNECT  IMF  FU1N1S  (Y  UR  N>":  ",CUNN*(l> 

7 1  flu  If  CUNIIt  (  I  )  -  "n  ”  OR  CUNNt  (  I  )  -  "N"  T HP  N  7210 

719"  INPUT  "WHA1  COLOR  LINE*1  (1-7)  "  ,LLULP<  (  I  ) 

72*  >o  NEXT  I 

*310  If  FI  "V"  THEN  7*00 

722*7  PRINT  "WHICH  SYMDUL  WOULD  YOU  Lilt  FOR  1  HE  MEASURED  DAIA: 

777*7  I  R1 M1  "  1  2  7  4" 

774*7  INPUT  "  Squpr  p  Dot  Triangle  Pirrlo  (1,2,7,  c.r  4):  ".ISvMOD 
7250  I NP  LI  I  "DU  YOU  WISH  TO  CONNECT  TOE  I  (UNIS  (  r  OK  N>  "  ,  CUNIIT  ( o ) 

726*7  IP  CUNIIt  O.i) -"o"  OR  CONNJ  ( *  i )  -  "N"  TUT  N  -.•(To 
7.770  INPUT  "WHaT  COLUR  LINE"*  (1-3)  "  ,  LCUl.  R  ( *  * ) 

7771  ’ 

72  7  2  ’ 

7780  ' 

7790  ’  calculate  scaling  Factors 
74*7*7 

7410  XSLl  -716/ (XMAX-XMIN)  :  '  716-  total  length  of  ■<  an  n 

7430  YSCL=  1 04  /  ( YMAX  -  YMIN)  :  '  lo4=  total  length  of  y  a::  l  s 

7470  SCREEN  1 : COLOR  0,1:CLS 
3440  ’ 

345*.i  draw  ti  cl  s 

3460  ’ 


A34 


’•'170 
34B0 
3500 
35  1 0 
752*> 
3570 
3540 
3550 
3570 
35B0 
3590 
3600 
36 1 0 
3620 
3630 
3640 
365*7 
3660 
3  6  BO 
3670 
7 1 1)0 
3710 
7  770 
7770 
3740 
3760 
37  70 
3780 
7  79*  > 
3800 
7010 

3  coo 
70  70 
3840 
7850 
3860 
30  7*7 
700*7 
7090 
79*  ii 7 
39  1  * ) 
39'.  *  > 
393*7 
794*  I 
395*7 
7960 
79  7*1 
790*7 
799'* 
4*7*70 
4*7  1  *7 

4  *' 3*7 
4  07*7 
4*74*7 
4*75*7 
4060 
4070 
408*7 
4090 
4  1 OO 


Y  INC  - 5*7 

IF  Y1NC  155  THEN  3530s' 

LINE  ( 6*  * , Y I NO ) - (55 , Y I NC ) :  ' 
YINC=Y1NC+YSCL* YT INI 
GOTO  3480 

X  IWC=6*7 

IF  X  INC  277  THEM  36*7*7: 

LINE  ( X  INC, 154) - ( X  INC, 159) : 

X 1NC-X 1NC  +  XBCL.*XT INI 
GOTO  3550 

'  draw  grid  1  i  nes 

IF  Y6INI-0  AND  X15INT=0  THEN  3790 
IF  YG 1NT =0  THEN  3720 
YG INC-50 

IF  YGINC  155  THEN  7710 

L INE  (60, YGINC) - (276, YGINC) ,  l  :  ' 
YG I NC  =  YG I NC  +  YSCL* YG I NT 
GOTO  366*7 


154-bottom  o f  graph 
hor  i  r  ontal 


276-r  i  alit  side  of  graph 
vertices* 


hor i r  on t  a 1 


XGINC-60 

IF  XG IN  T  —  0  THEN  3790 
IF  XG  I  NC  277  THIN  379*) 

LINE  ( XG I NC , 1 54 ) - ( XG I NC , 50  >  , 1 :  vertical 

XG1NC-XG1NC-*  XSCC*XGINT 
GOTO  3740 

draw  a"P5 

CINE  (60, 154)  -  (276, 154)  ,3:  '  ::  axis 

LINE'  (60,  154)-(6*7,50>  ,3:  '  y  axis 

label  endpoints 

TENET ( 1 > =S1RT (YMAX >: TEHIT <2) =B1RT (YMIN) : '  sot  up  temp,  arrays 
T  F.MFT  (7)  -SIR!  (  XMIN)  :  TF.MI  T  (4  )  =  STRT  <  XMAX  > 

XFL1S  (  1  )  -7:  YFUS  (  1  >  -  6:  XF  05,  (2)  -2*7:  Yf  05  (2)  -6 
XF'UO  (  3)  =22:  YPOG  (  7)  =8:  XF'OS  (4  )  =22:  YF'OS  (4  )  =  75 

FOR  1=  1  10  4:  r i ght  justify  4  labels 

JFOB-O 

FOR  Jr  1  111  HIM  IE  nr  1  (  I  )  )  -1 
AT  =11  IDT  (  IT.  MR  I  <  I  )  ,  LEN  <  TEITf  I  (  I  )  )  -,1F  05 , 1  > 

L  OCA  I K  XROS < I ) , YFUS < I > -JF OS: PRINT  AT 
JFOS  =  JF  OS+ l 
NE  X  T  J 
NE  X  I  I 


grapb  title 

xr OS=INI  (27-  (l  EIJ(T 1TLET) 72)  )  :  •  center  title 

LOCATE  7, XF US: FLINT  T1TLET 

a::es  tit]  es 

XFUS=1NT  (23- (LEN(XTITLT)  72)  >  :  ’  center  ::  axis  label 

LOCATE  23,  XF’OS:  PRINT  XTITLT 
FOR  I  =  1  TO  15: 


write  y  axis  vertically 


4110  Ai=MIDr (YTITL* , I , 1 ) 

4120  LOCATE  1+8,3 
4  1  To  PR INI  A* 

4140  NEXT  I 
4150  LOCATE  23,1 
4  1  60  ’ 

4170  '  post  data  points 

4 1  BO  ' 

4161  IF  F  4  = '  Y "  THEN  J  =  0  ELSE  J  =  1 
4182  FOR  JJ=J  TO  MODLN 

4190  FOR  1=1  TO  NFLT 

4191  VG(  I  )  =V  <  MODI-  (  JJ  >  ,  I+STF1-1  ) 

4200  XPT= (XG ( I > -XMIN) +XSCL+60 

4210  YRT=200-( ( YG ( I ) -YMIN) *YSCL+46> 

4220  ' 

4240  IF  ISYMtJJX  l  THEN  4280 

4250  LINE  1 XFT-1 , YF T-l > - (XPT+1 , YPT  +  1 )  ,2,B:  ’  rectangle 

42f.n  GOTO  4381 

42  70  ' 

42B0  IP  ISYM(JJ)  2  THEN  4320 

4290  PSET  (XPT, YPT) ,2: '  dot 

4300  GOTO  4381 

4  3 1  O  ’ 

4720  IP  ISYN(J.l)-  3  THEN  43B0 

4330  LINE  (XPT-7,  YT"r+7>  -  (  XF  T  ,  YFT-7)  :  '  triangle* 

4340  LINE  - (XPT + 7, YPT +3) 

4330  L I ME  -<XPT-3,YPT+3> 

4360  0010  4381 

43  71.1 

4700  CIRCLE  ( XPT , YF T ) , 3 , , , , 1 : ’  circle 

4781  NEXT  I 
4390  NEXT  JJ 
4  400  ' 

4410  '  connect  the  data  points 

4420  ' 

4421  IP  Ft  =  "Y"  T HF N  J  =  0  ELSE  J  =  1 

4422  FOR  JJ=J  TO  MODI. IT 

4423  IF  CONIT  t  ( J J ) = "N "  OR  CONNI  ( J  J  )  =  "n "  THEN  4520 

4424  FOR  1  =  1  TO  NfT.T  -1 

44  75  XG  <  I)  -0  (MODI  (  J  .1 )  ,I+STF'l  —  1)  : YG  < 1  +  1 >  =V  (MLJDL  (JJ  )  ,  I+S1P1  ) 
4  460  XPT  =  ( XG  U  > -XMIN)  *XSCL  +60 

4  4  7o  YPT =200  <  < YG ( I > - YM IN) « Y SOL  +  46) 

4  480  XFT 2= ( XG ( I + 1  >  -XMIN) *XSCI  * 6<  > 

44  90  YPT  2=20o-  (  (  Y ti  (  I  +  1  )  -YMIN)  «  Y SCI .  +  46) 

4500  LIME  (XF T , YP 1) - ( XP T 2 , YPT 2 ) ,LCOLR(JJ) 

4510  NEXT  I 

4520  NEXT  JJ 

4521  ' 

453o  ’  legend 

4540  ' 

4550  ' 

456o  '  pause 

4570  ' 

4500  IT  IN!  E YT  = " "  THEN  4580 
4585  FEY  ON 
4590  ' 

4600  GOTO  52 
4610  RETURN 
5000  ' 

5005  KEY  ON 


5010  GOTO  52 
5020  RETURN 

7000  ■*•*******•«***•*******•**•*•*«•««*««•**•««*********«*** 

7001  CLS 

7005  'Data  reduction  optioni  Gradient  and  constant  correction 

7010  INPUT  "Enter  -first  station  number,  correction  "jSF.Cl 

7020  INPUT  “Enter  last  station  number,  correction  "jSL,C2 

7030  INPUT  "Enter  constant  correction  “;C3 

7040  NS  »  SL-SF+1 

7050  NCOR  =  SL-SF 

7060  DCOR  =  (C2-C 1 > /NCOR 

7065  COUNT  =  0 

7070  FOR  I  =  SF  TO  SL 

7075  COUNT  =  COUNT  ♦  1 

7080  V(0,I>  =  V(0,I)  +C1  +DCOR« (COUNT- 1 >  +C3 

7081  'Round  off  calculated  values  using  integer  division 

7082  V  <  0 ,  I  >  =  V ( O , I > \ 1 
7090  NEXT  I 

7100  INPUT  "Correct  another  set  of  stations  (Y/N)"?  "(MSTAT# 
7110  IF  MSTAT*  =  "Y”  THEN  7000  ELSE  1562 
7120  RETURN 


A37 


(0,0,0) 


x 


where  Vp  =  generated  voltage 

=  source  or  sink  current  intensity 
p  =  resistivity  of  medium 

r.  s  A /( x  -  x  +  (y  -  y  *  (z  -  z 
1  p  s^'  KJp  7s^J  v  p  s^ 


References 

TeLford  et  al. ,  1976 
Corwin  et  ai. ,  L98L 

Point  Current  Sources 


PLATE  Al 


Line  source  is  parallel  to  y-axis,  centered  at  (xs,  ys) 


(0,0,0) 


where  Vp  =  generated  voltage 

I  =  current  per  unit  Length 
p  =  resistivity  of  medium 


References: 

Rao  et  al..  1970 

Broughton  Edge  and  Laby,  1931 

Semenov,  1974 

Horizontal  Line  Source 


PLATE  A2 


(0,0,0) 


x 


V  (x  v  .  0)  >  E  r2  xo  sin  •  -  h  c°s  ° 
P  P  P  0  0 - 


(X, 


2  .2.3/2 

yp  -  h  ) 


where  Vp  =  generated  voLtage 

E  =  maximum  potential  at  surface  of  sphere  (at  9=0) 
r0  =  radius  of  sphere 
h  >>  r0 

resistivity  of  sphere  =  0 


References 

Petrowsky,  1925 
Semenov,  1974 

_ _  Spherical  Source 

PLATE  A3 


Cylinder  is  of  infinite  length  and  parallel  to  y-axis. 
Polarization  is  uniform. 


where  Vp  =  generated  voitage 

I  =  current  per  unit  of  circumference 
p  -  resistivity  of  medium 
r0  =:  radius  of  cylinder 
h  >>  r0 

resistivity  of  cylinder  -  0 

References 

Murty  and  Haricbaran,  1985 
Bhattacharya  and  Roy,  1981 

Semenov,  197 h  Horizontal  Cylindrical  Source 


PLATE  AA 


1  2 


P  =  pi  for  y  >  0 

y 


tan-1  F^  +  tan'1 


%> 


where  =  generated  voltage  in  region  1  (y  >  0)  or  2  (y  <  0) 

2 

S  =  source  strength  b  >  a  >  0 

f  _  _ Lx  ♦  1/2)  b _ 

L  y  [ (x  4.  l/2)2  +  y2  +  b2]I/2 

F _ (x  +  t/21  a _ 

2  y  [(x  f  J./2)2  4-  y2  4-  a2 ] 1/2 

F  _  _ (x  -  1/21  b _ 

3  y  t(x  -  t/2)2  4-  y2  4-  b2]1/2 

F _ (x  -  i/2)  a _ 

4  y  [(x  -  9. / 2 ) 2  4-  y2  t-  a2]l/2 
Notes:  1.  f-’L,  F2,  F 3 ,  F4  are  expressed  in  radians 

2.  V  is  positive  for  x  and  y  positive  and  s  positive 

3.  Interchange  x  and  y  for  plane  parallel  to  y-axis 

References:  Fitterman,  1979 

Fitterman  and  Corwin,  1982  Vertical  Dipolar  Sheet  Source 
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APPENDIX  B 


SELF-POTENTIAL  FIELD  PROCEDURE  AND  DATA  REDUCTION 


1.  The  self-potential  (SP)  survey  procedures  described  below  have 
been  found  to  give  data  of  acceptable  reproducibility  for  seepage  studies, 
mineral  and  geothermal  exploration,  and  other  applications  of  the  SP 
method.  More  elaborate  procedures  involving  multiple  base  or  measurement 
stations,  repeated  reoccupation  of  stations,  monitoring  of  electrode 
temperatures,  etc.  may  result  in  better  data  quality,  but  the  increased  time 
and  expense  involved  must  be  compared  with  anticipated  improvements  in  data 
quality  to  determine  whether  such  procedures  will  be  cost-effective  for  a 
particular  application. 

2.  The  procedures  given  are  necessarily  general  in  nature,  and  may 
have  to  be  modified  somewhat  for  special  conditions  such  as  frozen  or  very 
rocky  soil,  high  telluric  or  stray  current  noise  levels,  etc.  It  is 
recommended  that  a  telluric  monitor,  consisting  of  a  battery-operated  strip 
chart  recorder  continuously  measuring  the  potential  between  a  pair  of 
stationary  electrodes,  be  installed  in  the  survey  area  to  detect  SP  time 
variations  that  could  erroneously  be  recorded  as  spatial  anomalies. 

3.  The  recommended  "fixed-base"  survey  configuration  has  been  found 
to  give  much  more  reproducible  results  than  the  "leapfrog"  or  "gradient" 
configuration,  in  which  a  length  of  wire  eqiJ^.  to  the  station  separation  is 
moved  along  the  survey  line  and  the  total  SP  reading  at  a  given  station  is 
calculated  by  successive  addition  of  individual  dipole  readings.  The 
cumulative  error  inherent  in  this  process  has  been  found  to  generate 
significant  spurious  "anomalies",  and  it  is  recommended  that  this  type  of 
configuration  be  used  only  in  situations  where  it  is  impossible  to  deploy 
the  longer  wire  lengths  needed  for  the  fixed-base  technique. 
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Survey  Procedure 

(using  copper-copper  sulfate  electrodes  and  digital  multimeter) 

1.  Record  initial  information  on  data  sheet  (see  sample  following 
text). 

2.  Reel  and  meter  check 

a.  Check  wire  insulation  for  breaks. 

b.  Measure  and  record  reel  wire  resistance  (infinite  reading 
indicates  open  wire). 

c.  Measure  resistance  between  either  end  of  wire  and  frame  of 
reel  for  short  circuits  between  wire  and  frame.  Reel  must  not 
contact  ground  during  readings  if  short  is  present. 

3.  Initial  electrode  polarization  check 

Measure  and  record  voltage  between  all  three  electrode  pairs  in 
copper  sulfate  bath  as  shown  on  sample  data  sheet.  These  values 
will  be  used  later  for  drift  corrections.  Measurements  should  be 
made  on  meter  range  having  0.1  or  1.0  mV  resolution. 

4.  Base  electrode  installation 

a.  Select  base  station  site  where  conditions  are  unlikely  to 
change  during  the  survey  period.  If  tying  into  previous 
station,  re-occupy  previous  hole. 

b.  Tie  off  end  of  reel  wire  and  attach  to  base  electrode  terminal. 
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c.  Oig  hole  into  moist  soil  below  surface  layer.  Do  not  add 
water  to  hole.  Clean  loose  soil  and  rocks  from  bottom  of  hole. 

d.  Install  base  electrode  by  pushing  down  and  rotating  to  ensure 
good  soil  contact.  Pile  soil  around  electrode  and  tamp  down 
firmly.  Be  sure  that  no  soil  contacts  electrode  terminal. 

e.  Place  sun  shade  over  electrode  and  weight  down  with  soil  or 
rocks. 

f.  Be  sure  to  note  base  station  location  on  data  sheet. 


5 .  Measurement  procedure 

a.  Minimum  equipment  needed  by  survey  operator  includes 
multimeter,  data  sheet,  digging  tool,  flagging  tape,  and 
"portable  reference  electrode"  carried  in  a  bath  of  copper 
sulfate  solution. 

b.  Transport  reel  and  equipment  above  to  first  measuring  station. 

c.  Select  station  location  in  area  of  uniform,  undisturbed  soil, 
if  possible.  Dig  hole  to  moist  soil  below  surface  layer  and 
clean  loose  soil  and  rocks  from  bottom  of  hole. 

d.  Insert  measuring  electrode  in  hole  with  firm  downward  pressure 
and  rotate  to  ensure  good  soil  contact.  Hold  electrode  firmly 
in  place  or  pack  soil  around  it  to  hold  it  upright. 

e.  Connect  positive  lead  of  meter  to  electrode  and  negative  lead 
to  connector  on  reel.  Shade  electrode  from  sun. 

f.  Make  SP  reading: 

1)  Read  potential  on  meter  range  having  0.1  or  1.0  mV 
resolution. 

2)  Read  value  for  at  least  10  to  20  seconds  to  check  for 
drift  or  telluric  fluctuations.  If  reading  fluctuates, 
read  long  enough  to  obtain  a  reasonable  average  value. 
Record  final  value  in  mV,  with  +  or  -  sign  to  indicate 
polarity. 

3)  Measure  contact  resistance  between  electrodes  using 
appropriate  kOhm  range  of  multimeter.  Resistance  should 
be  measured  for  only  a  few  seconds  or  less  to  avoid  undue 
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polarization  of  electrodes.  If  resistance  differs 
considerably  from  previous  values,  try  re-seating 
electrode  or  relocating  station.  Note  that  large  SP 
values  may  produce  erroneous  resistance  readings  on  the 
meter.  In  such  cases,  measure  resistance  with  meter  leads 
reversed  and  record  average  value.  Return  meter  to 
voltage  range  after  reading  resistance. 

4)  Note  soil  conditions,  vegetation,  cultural  features,  etc. 
in  "Comments"  column  of  data  sheet. 

5)  Remove  electrode  and  clean  loose  soil  from  ceramic  tip. 
Carry  electrode  in  shaded  location. 

6)  If  there  is  any  chance  that  the  station  will  be  a  future 
tie-in  point  or  will  be  re-occupied  for  any  other  reason, 
put  flagging  tape  in  hole  to  mark  it  and  re-fill  it  to 
keep  the  soil  from  drying. 

g.  Move  to  next  station  and  repeat  procedure  above.  About  once 
every  half  hour  to  hour,  or  at  the  end  of  a  survey  line, 
measure  the  potential  between  the  measuring  electrode  and  the 
portable  reference  electrode  in  the  copper  sulfate  bath  and 
record  as  shown  on  the  sample  data  sheet.  These  readings  will 
be  used  for  later  drift  corrections,  as  discussed  below. 

h.  After  completing  survey  line,  reel  wire  in  and  return  to  base 
station.  Remove  and  clean  base  electrode  and  re-measure 
polarization  between  all  three  electrode  pairs  in  copper 
sulfate  bath  as  shown  on  sample  data  sheet. 


Data  Reduction  Procedures 

4.  It  is  assumed  that  the  polarization  between  the  base  and  the 
portable  reference  electrodes  remains  constant,  and  that  the  measuring 
electrode  drifts  relative  to  these  two.  Therefore  two  electrode  corrections 
are  superimposed  on  the  measured  values:  an  initial  polarization  between  the 
base  and  measuring  electrodes  (which  is  assumed  to  remain  constant 
throughout  the  survey)  and  a  drift  of  the  measuring  electrode,  which  is 
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determined  by  periodically  measuring  its  potential  relative  to  the  portable 
reference.  If  in  a  given  case  it  is  apparent  that  the  portable  reference 
has  drifted  with  respect  to  the  base,  then  drift  corrections  should  be  made 
using  only  the  initial  and  final  readings  between  the  base  and  measuring 
electrodes  instead  of  the  procedure  below. 

5.  On  the  sample  data  sheet,  the  initial  polarization  between  the 
base  (#1)  and  measuring  (#3)  electrodes  is  +5  mV.  This  is  corrected  by 
subtracting  3  mV  from  all  the  measured  values. 

6.  At  the  first  reading  between  the  portable  reference  and  measuring 
electrodes  at  0933  the  #2(-)  #3(+)  reading  (#2  electrode  to  the  negative 
voltmeter  lead;  #3  electrode  to  the  positive  lead)  has  gone  from  an  initial 
reading  of  +3  mV  to  a  new  value  of  +10  mV,  i.e.,  the  measuring  electrode  is 
7  mV  more  positive  than  it  was  at  the  beginning  of  the  survey.  Therefore, 
an  additional  7mV  is  subtracted  from  the  initial  -5  mV  correction  for  a 
total  electrode  correction  of  -12  mV.  Electrode  corrections  between  drift 
measurements  are  interpolated  linearly. 

7.  At  1005,  the  #2(-)  #3(+)  reading  is  +8  mV,  so  the  measuring 
electrode  now  is  5  mV  positive  with  respect  to  the  initial  reading  of  +3 
mV.  Therefore,  5  mV  is  subtracted  from  the  original  -5  mV  polarization 
correction,  giving  a  total  electrode  correction  of  -10  mV.  Again,  electrode 
correction  values  are  interpolated  back  to  the  previous  drift  reading. 
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8.  If  the  base  station  was  at  a  point  where  SP  had  been  measured 
previously,  the  SP  value  at  that  point  must  be  added  to  all  the  measured  and 
drift -corrected  values  as  a  tie-in  correction  (e.g.,  the  +12  mV  value  on  the 
sample  data  sheet).  The  final  corrected  SP  value  is  the  sum  of  the 
measured,  electrode  drift  correction,  and  tie-in  values. 


SAMPLE 


DATE 


2  Of. til  l°lt° 
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PAGE 
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Comprehensive  Seepage  Assessment:  Beaver  Dam,  Arkansas 


J.L.  Llopi* 

U  S.  Army  Engineer!  Waterways  Experiment  Station,  Vicksburg, 
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C.M.  Deaver 
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Mississippi 

S.C.  Hartung 
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SYNOPSIS:  A  general  philosophy  of  the  role  of  engineering  geology  end  engineering  geophysics  in 
seepage  assessaent  is  presantad.  Practical  application  of  this  philosophy  is  illustrated  by  e  case 
history.  A  large  dike  continues  to  havs  anomalous  seepage  in  spite  of  pre-construction  and  post¬ 
construction  grouting.  The  dike  is  founded  over  a  graben  of  cavernous  limestone  with  about  a  200-ft. 
vertical  offset  along  the  bounding  fault  zones,  which  art  horizontally  separated  by  about  1000  ft. 
Objectives  of  the  seepage  assessment  program  were  to  define  the  geological  and  hydrological  condi¬ 
tions  beneath  the  dike  in  sufficient  detail  to  allow  rational  remedial  planning. 

Integration  of  results  of  a  geophysical  investigation  with  the  overall  assessment  program  is 
emphasized:  preliminary  interpretation  of  the  geophysical  results  Is  used  to  site  new  piezometers; 
detailed  analysis  of  the  geophysical  results  Is  used  to  site  exploratory  borings;  feedback  from  ex¬ 
ploratory  borings  and  new  piezometers  is  used  to  refine  geophysical  interpretation. 


INTRODUCTION 

Background 

Earth  dams  and  dikes  are  expected  to  seep,  and 
thair  designs  include  drainage  systems  to  col¬ 
lect  and  discharge  seepage  water  into  the 
downstream  channel.  Sometimes,  however,  seepage 
occurs  in  an  unplanned  manner,  exceeding  the 
capacity  of  the  drainage  system  or  along  a  path 
not  considered  in  the  seepage  design.  Excessive 
unplanned  seepage  may  be  just  unsightly  (though 
possibly  disconcerting  to  the  public) ,  or  it  may 
threaten  the  integrity  of  the  embankment.  In 
these  cases  it  may  be  necessary  to  conduct  a 
seepage  assessment  program  to  detect  and  map 
seepage  paths  in  order  to  more  rationally  plan 
remedial  measures. 

Dike  1,  at  Beaver  Dam  has  been  experiencing  a 
general  increase  In  seepage  rates  since  initial 
reservoir  filling  in  1966.  Recently  however,  the 
proliferation  of  seepage  exits  along  the  toe  of 
Dike  1  has  prompted  the  Little  Rock  District, 
V.S.  Army  Corps  of  Engineers  (SWL)  to  undertake 
a  comprehensive  seepage  assessment  program.  This 
program  consisted  of  examining  the  project  his¬ 
tory,  mapping  and  topographic  survsying,  surface 
geophysical  testing,  extending  the  piezometer 
network  (including  drilling,  sampling  and 
testing),  exploratory  drilling,  seepage  flow 
meesureoents,  planning  for  and  installing  an 
automated  piezometer  and  flow  measurement  data 
acquisition  systsm,  and  remsdial  measurs 
analysis.  In  support  of  this  effort  the  U.S. 

Army  Engineer  Waterways  Experiment  Station  (WES) 
was  raquestad  to  perform  a  detailed  geophysical 
Investigation  of  the  dike  and  its  foundation. 

PMTP9M 

The  purpose  of  this  paper  is  to  present  the 
general  philosophy  of  a  seepage  assessment 
prograa  conducted  at  Beaver  Dam,  Arkansas. 
Described  are  the  various  phases  of  the  prograa 
and  how  thay  are  Integrated  to  allow  for  a  more 
rational  approach  to  remedial  planning. 


Sits  Location  and  Description 

Beaver  Dam  la  located  on  the  Whits  River  at 
rivar  mils  609.0  in  Carroll  County,  Arkansas, 
spproxiaataly  6  alias  northwest  of  Eureka 
Springs,  Arkansas.  Beaver  Dam  is  a  straight, 
gravity-type,  concrete  structure  flanked  to  the 
north  by  an  earth  embankment  and  thraa  saddle 
dikes.  The  location  of  Dike  1  relative  to  the 
concrete  dam  and  main  embankment  is  shown  in 
Figure  1.  The  reservoir  (Beaver  Lake)  is  used 
for  flood  control,  power  generation,  and  water 
supply.  Construction  of  the  dam  was  started  in 
November  1960  and  ended  in  June  1966.  Dike  1  is 
approximately  1,000  ft  in  length  and  30  ft  high. 
The  top  of  the  conservation  pool  is  elevation 
1,120  ft  whila  the  top  of  the  dike  is  elevation 
1,162  ft.  Dike  1  is  founded  on  severely 
weathered  1 imestone  and  ia  experiencing  seepage 
from  various  axlts. 


GEOLOGY 

Ssneral  Qtolpqy 

Beaver  Dam  and  reservoir  area  are  located  in  an 
area  known  as'the  Ozark  uplift,  a  region  con¬ 
sisting  of  flat-lying  sedimentary  rocks  composed 
chiefly  of  limestone  and  dolomitic  limestone. 

The  strata  are  nearly  horizontal  over  the 
greater  part  of  the  area  but  ara  locally 
deformed  by  simple  dislocations  along  southwest- 
northeast  trending  normal  faults  and  shallow 
basins  that  in  places  of  ara  of  considerable 
magnitude. 

Bhvsloqraphy 

The  upland  area  around  the  dam  is  a  part  of  tha 
Springfiald  Plataau,  tha  surface  of  which  la 
developed  at  approximata  elevation  1500  ft  tha 
cherty  limastona  of  the  Boone  Formation.  In  tha 
dam  and  raservoir  arsa,  tha  Whits  Rivsr  has  cut 
a  channsl  approximately  600  ft  in  depth.  This 
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Figure  1.  Location  of  Dike  1  relative  to  concrete  dam  and  eain  embankment 


inciaion  into  the  Plateau  surface  has  resulted 
in  a  deeply  and  intricately  dissected  type  of 
topography.  The  entrenched  river  follows  a 
meandering  course  across  the  area. 

Stratigraphy 

Five  formations  are  exposed  at  the  dam  site. 

They  are  (moving  upsection)  the  Powell  Forma¬ 
tion,  the  Cotter  and  the  Jefferson  City  Forma¬ 
tions  of  the  Jefferson  City  Group  which  is  of 
Ordovician  age,  the  Chattanooga  Formation  of 
Devonian  age,  and  the  Boone  Formation  of  Missis- 
sippian  age.  The  Chattanooga  and  the  Boone  For¬ 
mations  are  generally  above  reservoir  level  ex¬ 
cept  in  the  vicinity  of  the  left  abutment  of  the 
dam  and  Dike  1  where  the  units  are  downfaulted. 
In  the  vicinity  of  the  dam  site,  the  Boone  For¬ 
mation  caps  the  higher  ridges  and  forms  the 
sides  of  the  valley  down  to  approximate  eleva¬ 
tion  1200  ft.  Beneath  this  lies  the  Chattanooga 
Shale  member  (Chattanooga  Formation) ,  which  in 
turn  is  underlain  by  its  Sylamore  Sandstone  mem- 
ber.  Beneath  these  and  forming  the  valley  walls 
below  elevation  iibo  ft  and  underlying  the 
greater  part  of  the  valley  bottom  are  limestones 
and  dolomitic  limestones  of  the  Jefferson  City 
Group  (Design  Memorandum  No.  5,  1959). 

Structural  oaoloov 

The  general  structural  geology  of  the  region  is 
that  of  flat  lying  rocks  which  are  locally 
deformed  by  simple  dislocations  along  southwest- 
northeast  trending  normal  faults  that  extend  for 
considerable  distances,  and  by  monoclines,  low 
domes,  and  shallow  basins.  The  Beaver  Dam  site 
lies  near  the  northeast  end  of  a  very  gentle, 
shallow,  elongate,  northeast-southwest  trending 
structural  basin  known  as  the  Price  Mountain 
syncline.  This  basin  is  often  faulted  in  areas 
where  the  downfolding  is  most  pronounced.  In  the 
9r**ter  part  of  the  lover  end  of  the  reservoir. 


Ordovician  strata  underlie  the  valley  floor  and 
extend  up  the  sides  of  the  valley  to  about 
elevation  X.1S0  ft.  Overlying  these  and  almost 
everywhere  above  pool  level  are  formations  of 
Devonian  and  MisBissippian  age.  In  localized 
areas,  these  units  have  been  downfaulted  to  form 
a  part  of  the  foundation  under  the  most 
topographically  desirable  dam  sites  in  the  val¬ 
ley.  This  is  the  case  at  Dike  1. 

Dike  1  Foundation  Materials 

Figure  2  shows  the  foundation  materials  underly¬ 
ing  Dike  1.  Dike  1  is  founded  on  a  downfaulted 
block  of  the  Boone  formation.  This  downfaulted 
block  (graben)  extends  approximately  between 
station  63+00  at  the  northern  end  to  approximate 
station  75+00  at  the  southern  end,  a  total  dis¬ 
tance  of  approximately  1,200  ft.  The  graben  is 
bounded  by  steeply  dipping  normal  faults  on 
either  side  trending  roughly  in  a  northeast- 
southwest  direction.  The  vertical  displacement 
of  these  faults  is  approximately  200  ft.  Cores 
of  the  rock  adjacent  to  the  northern  fault  zone 
show  evidence  of  fracturing;  however,  the  frac¬ 
tures  appear  to  be  filled  or  cemented  and  sound. 
Boring  information  from  the  southern  fault  zone 
area  indicates  the  presence  of  many  clay-filled 
cavities.  The  southern  fault  gouge  does  not  ap¬ 
pear  to  have  the  sue  degree  of  soundness  as  the 
northern  fault  zone.  The  Boone  Chert  which  makes 
up  the  foundation  of  Dike  1  can  be  divided  into 
two  distinct  sub-units.  The  upper  sub-unit  of 
the  Boone  Chert  (estimated  thickness,  100  ft)  is 
composed  of  calcium  carbonate  and  chert  which 
upon  weathering  has  resulted  in  the  removal  of 
calcium  carbonate  and  left  a  spongy,  vuggy, 
residual  material  that  is  predominately  chert. 
The  lover  sub-unit  of  the  Boone  Chert  (estimated 
thickness,  60  ft)  Is  also  composed  of  calcium 
carbonate  and  silica;  however,  this  sub-unit  is 
characterized  as  being  slightly  weathered  to  un¬ 
weathered  and  contains  more  crystalline  calcium 
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Figure  2.  Foundation  materials  underlying  Dike  1 


carbonate.  The  lower  sub-unit  is  noderataly  to 
closely  jointed  and  this  jointing  has  allowed 
the  passage  of  water  which  has  led  to  the  dis¬ 
solution  of  calcium  carbonate  and  in  turn  has 
resulted  in  open  channels  and  cavities. 

Underlying  the  Boone  Chert  Unit  is  the  St.  Joe 
Limestone,  described  as  non-cherty,  gray  to 
green-gray,  crystalline,  very  fossiliferous,  and 
containing  numerous  thin  shale  seams  and  part¬ 
ings.  Underlying  the  St.  Joe  Limestone  is  the 
Chattanooga  Shale  described  as  black,  firm,  and 
fissile.  The  shale  is  considered  to  be  an  effec¬ 
tive  barrier  to  any  downward  movement  of  ground 
water . 


SEEPAGE  HISTORY  OF  DIKE  1 

fre-Congtruction  grout  Curtain 

The  foundation  materials  of  Dike  1  ware  recog¬ 
nized  as  being  susceptible  to  seepage  during  the 
early  phases  of  the  site  selection.  In  June  1959 
it  was  decided  that  an  economical  solution  to 
prevent  a  potential  seepage  problem  was  to  in¬ 
stall  a  grout  curtain.  The  grout  curtain  con¬ 
sisted  of  two  lines  of  holes  spaced  5  ft  apart 
with  10-ft  hole  spacings  which  extended  to  a 
depth  of  5  ft  below  the  top  of  sound  or  un¬ 
weathered  rock  at  all  locations  except  between 
stations  72+70  and  74+70,  where  the  grout  cur¬ 
tain  was  extended  deeper  (16  to  65  ft)  into 
sound  rock  (Figure  2).  A  total  of  284  holes 
(26,200  linear  ft)  were  drilled  in  this  grouting 
program.  The  grout  (31,000  cu.  ft)  was  placed  by 
gravity  flow  (Reconnaissance  Report,  Beaver  Das, 
1984) . 

During  initial  filling  of  the  reservoir  (April 
1966)  seepage  was  detected  in  a  small  valley 
downstream  of  Oiks  1.  The  reservoir  pool  eleva¬ 
tion  at  this  point  was  1,102+  ft  and  the  seep 
was  flowing  at  a  rate  of  150-200  gpm  By  June 
1966  the  reservoir  elevation  waa  1114  ft  and 
eight  additional  seeps  wars  detected  with  a  com¬ 
bined  flow  rata  of  approximately  400  gpm  By  the 
time  remedial  grouting  operations  were  initiated 
in  1968,  the  combined  flow  rata  of  these  seeps 


had  risen  to  approximately  800  gpm  Conclusions 
from  studies  conducted  at  Dike  1,  including  flow 
measurements  and  dye  tracing,  indicated  seepage 
was  coming  from  the  lake  through  two  possible 
passages,  either  beneath  the  grout  curtain 
through  open  cavities  in  the  foundation  rock,  or 
along  the  top  of  rock  or  both.  Seepage  was  oc¬ 
curring  along  the  entire  length  of  Dike  1  with 
the  most  concentrated  flow  occurring  in  the 
vicinity  of  station  71+00  near  the  southern  por¬ 
tion  of  the  dike  (Reconnaissance  Report,  1984)  . 

Several  possible  explanations  why  the  pre¬ 
construction  grout  curtain  did  not  perform 
satisfactorily  are  as  follows: 

a.  Grout  holes  were  not  drilled  deep  enough 
to  sound  rock  to  intercept  open  joints. 

b.  Since  grout  was  placed  by  gravity  flow, 
it  is  possible  many  small  cavities  and 
joints  were  not  filled. 

c.  Grout  was  too  thick  to  enter  some  of 
the  cavities  and  joints. 

d.  Since  drilling  was  performed  with 
tricone  roller  bits  using  compressed  air 
to  remove  cuttings,  the  cuttings  could 
have  plugged  some  of  the  cavities 
preventing  them  from  being  grouted. 

a.  Many  cavities  and  joints  could  have  been 
missed  altogether  because  of  the  grout 
hole  spacing. 

Early  Seepage  Flow  Studies 

Flow  measurements,  exploratory  drilling,  pres¬ 
sure  tests,  and  dya  and  temperature  tests  were 
conducted  from  the  time  of  leakage  (1966)  until 
1968  to  determine  the  extent  and  routes  of 
seapaga  through  and  beneath  the  dike  and  to  for¬ 
mulate  possible  remedial  measures.  These 
measurements  ware  accomplished  by  installing  two 
weira,  a  Parshell  flume,  and  twenty-seven 
piezometers.  The  data  suggested  that  the  leakage 
was  issuing  both  through  cavities  below  the 
original  grout  curtain  and  along  the  top  of 
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crystallina/weathered  rock  interface.  It  was 
concluded  that  seepage  occurred  along  the  entire 
length  of  Dike  1  and  to  the  fault  zone  beneath 
the  sain  embankment  at  station  71+00,  with  the 
greatest  seepage  occurring  along  the  shortest 
flow  path  in  the  vicinity  of  station  71+00. 

Remedial  Grout  Curtain 

During  the  period  July  1968  to  December  1971  an 
extensive  grouting  program  was  conducted  in  an 
effort  to  abate  the  seepage  occurring  at  the 
dike.  The  program  consisted  of  10,040  linear  ft 
being  drilled  in  226  holes.  Also,  18,900  cubic 
ft  of  grout  solids  were  pressure  injected  into 
these  holes  with  the  heaviest  grout  takes  occur¬ 
ring  in  an  area  between  stations  70+50  and  72+00 
(Figure  2) .  Problems  encountered  during  the 
grouting  operations  were  collapsing  boring  walls 
(cave-in),  insufficient  seating  of  casing,  and 
incapability  of  grout  pump  to  grout  some  large 
cavities  to  refusal. 

As  a  result  of  the  remedial  grouting  program, 
seepage  was  reduced  to  approximately  450  to  500 
gpm  for  mid-pool  elevations  (1120-1110  ft),  a 
decrease  in  flow  of  10  to  IS  percent.  During 
the  period  1971  through  1984  piezometers  were 
manually  read  approximately  twice  a  year  by  SWL 
personnel  while  the  Parshall  flume  was  read  on  a 
monthly  basis  by  project  personnel.  During  a  pe¬ 
riodic  inspection  in  1980,  a  new  seepage  area 
was  located  on  the  downstream  right  abutment  of 
the  dike.  This  prompted  SWL  personnel  to  under¬ 
take  an  effort  to  locate,  inspect,  and  describe 
all  known  seepage  exits. 


Dan  Safety  Assurance  Program 

When  the  U.S.  Army  Engineer  Southwest  Division's 
(SWD)  Division-wide  Master  Plan  for  the  Cam 
Safety  Assurance  Program  was  submitted  in  1981, 
Beaver  Dam  was  listed  as  requiring  studies  for  a 
Reconnaissance  Report  under  designated 
priorities  of  spillway  adequacy  and  major 
seepage.  The  Reconnaissance  Report  (May  1984) 
concluded  that  seepage  at  Dike  1  would  increase 
to  near  pre-grouting  flows  (800+  gpm.)  during  a 
Spillway  Design  Flood  (Probable  Maximum  Flood, 
pool  elevation  1,139.9  ft)  and  continue  flowing 
at  this  rata  even  after  the  flood  receded  due  to 
expansion  of  existing  cavities.  This  conclusion 
was  proven  to  be  valid  on  23  December  1984  when 
a  Pool  of  Record  (el  1,130.4  ft)  occurred. 

During  the  emergency  flood  procedure  inspection 
on  that  date  the  project  superintendent  observed 
a  new  seepage  exit  500  feet  downstream  from  Dike 

1  with  a  flow  rate  of  approximately  25  gpm.  The 
alarming  factor  at  the  newly  discovered  exit 
however,  was  the  large  amounts  of  detrital 
material  (sediment),  ranging  from  clay-  to 
gravel-size  being  discharged  in  the  flow,  i.e. 
muddy  water.  Another  new  seep  was  discovered  on 

2  January  1985  near  the  left  dike/abutment  con¬ 
tact  at  approximate  elevation  1,106  ft  (Figure 
3) .  Water  from  this  new  seep  was  described  as 
jetting  vertically  with  a  flow  rate  of  ap¬ 
proximately  7  gpm.  at  pool  elevation  1,125.1  ft 
(Feature  Design  Memorandum,  1987).  The  1984 
Reconnaissance  Report  recommended  that  a  seepage 
investigation  be  undertaken  to  determine  the 
location  and  extent  of  seepage  and  develop 
remedial  measures  to  control  seepage  at  Dike  1. 


Figure  3.  Location  of  seepage  exits  and  proposed  piezometer  and  exploratory  borings. 
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A  combined  SWD/SWL/WES  meeting  was  held  at 
Beaver  Dam  on  14-16  January  1985  to  discuss  the 
proposed  seepage  investigations,  which  were  in¬ 
itiated  in  February  1985. 


SEEPAGE  INVESTIGATIONS 
General 

The  new  muddy  seepage  exit  below  Dike  1  dis¬ 
covered  during  a  pool-of-record  (1,130.4  ft)  in 
December  1984  not  only  substantiated  the  need 
for  seepage  investigations,  but  also  added  an 
element  of  urgency  and  a  necessity  to  expedite 
the  investigations,  and  recommendations  of 
measures  to  control  seepage.  In  1985  the 
monitoring  /inspection  of  instrumentation  and 
seepage  was  revised  to  more  frequent  scheduling, 
especially  for  pool  levels  above  elevation  1,128 
ft.  The  action  having  the  greatest  impact  on 
project  operations,  due  to  severe  seepage,  is 
the  request  and  approval  for  a  deviation  (loss) 
of  the  authorized  flood  storage  pool  in  Beaver 
Lake  from  elevation  1,130  to  1,128  ft  until  the 
seepage  the  seepage  problem  is  resolved.  The 
maj  or  «l«nents  of  the  seepage  investigation  are 
described  below. 

Geophysical  Investigations 

In  March  1985  WES  personnel  conducted  a 
geophysical  investigation  at  Dike  1.  Several 
geophysical  methods  were  used  for  this  study  in¬ 
cluding  self-potential  (SP) ,  electrical  resis¬ 
tivity,  electromagnetic  induction  (EM) ,  seismic 
refraction,  magnetic  profiling,  and  borehole 
water  conductivity/temperature  measurements.  The 
objectives  of  the  geophysical  investigation  were 
to  (a)  detect,  map,  and  monitor  seepage  through 
the  foundation  of  Dike  1,  (b)  delineate  geologic 
structure  beneath  and  immediately  adjacent  to 
Dike  1,  and  (c)  provide  input  to  the  planning  of 
remedial  measures. 

The  geophysical  methods  necessary  for  a  seepage 
analysis  are  not  difficult  to  use.  However,  a 
geophysical  survey  program  must  be  planned 
based,  to  the  maximum  extent  possible  on 
knowledge  of  the  (1)  surface  geometry  of  the  dam 
and  associated  features,  (2)  design  and  con¬ 
struction  details  of  the  structure,  and  (3)  the 
geology  of  the  foundation  and  abutments  (Butler, 
1985) . 

The  primary,  long-term  geophysical  method  was 
the  self -potential  (SP)  method,  which  was 
monitored  throughout  the  duration  of  the  inves¬ 
tigation.  SP  data  were  obtained  during  various 
pool  levels  to  determine  relationships  of 
seepage  flows  and  pool  levels.  The  SP  arrays 
were  installed  by  SWL  personnel  in  February 

1985.  Initial  SP  readings  vere  made  by  WES  in 
March  1985.  Subsequent  readings,  during  various 
pool  levels,  were  taken  by  SWL  personnel  and 
forwarded  to  WES  for  interpretation.  Detailed 
results  of  the  geophysical  studies  pertaining  to 
the  seepage  investigation  are  presented  in 
Supplement  No.  1  of  the  Reconnaissance  Report, 

1986. 

Additional  geophysical  studies  were  conducted  at 
Beaver  Dam  in  conjunction  with  the  Corps  of 
Engineers'  Repair,  Evaluation,  Maintenance,  and 
Rehabilitation  Program  (REMR) ,  These  REMR  spon¬ 


sored  geophysical  tests  included  high-resolution 
seismic  reflection,  ground-penetrating  radar, 
microgravimetry,  and  additional  SP  investiga¬ 
tions.  The  ground-penetrating  radar  survey  vas 
conducted  in  September  1985  with  an  additional 
survey  being  conducted  in  February  1986.  In 
August  1986  a  high-resolution  seismic  reflection 
and  a  "low  pool  level*  SP  survey  was  conducted. 
Detailed  results  of  these  geophysical  surveys 
are  presented  in  the  Feature  Design  Memorandum 

1987. 

The  geophysical  investigation  was  successful  in 
delineating  the  fault  zones  bounding  Dike  1, 
which  are  believed  to  act  as  channels  for  lake 
water  to  exit  downstream,  as  well  as  identifying 
other  faults  which  were  not  previously  known  to 
exist.  The  tests  also  identified  fractured  and 
saturated  zones  as  well  as  determining  the  ver¬ 
tical  extent  of  the  weathered  Boone  Chert.  The 
tests  also  indicated  that  seepage  is  occurring 
along  the  entire  length  of  the  dike.  The 
geophysical  tests  suggest  that  both  axial  and 
transverse  seepage  flows  are  occurring  along  the 
south  fault  zone,  but  that  the  north  fault  zone 
is  relatively  tight  (impermeable)  to  those 
flows.  Based  on  results  of  the  geophysical  tests 
an  integrated  seepage  map  was  produced  showing 
that  seepage  flows  are  moving  primarily  in  an 
east-southeasterly  direction  with  the  greatest 
flows  occurring  between  stations  69*00  and 
73+00,  and  along  the  south  fault  zone  (Figure 
4)  • 

Exploratory  Borinas 

Twenty-five  exploratory  borings  were  drilled 
along  the  upstream  crest  of  Dike  1  and  its  abut¬ 
ments  during  the  period  April  1986  to  August 
1987.  The  primary  purpose  of  these  borings  was 
to  delineate  the  limits  and  geologic  charac¬ 
teristics  of  the  downthrown  faulted  block  of  the 
Boone  Formation  beneath  Dike  1  and  the  North  and 
South  fracture  zones  that  bound  the  Dike. 
Originally,  the  boring  locations  were  selected 
based  on  areas  that  had  experienced  high  grout 
takes  during  the  previous  grouting  program. 
However,  locations  for  the  borings  were  later 
changed  to  take  advantage  of  information  ob¬ 
tained  from  geophysical  testing.  Based  on 
results  of  the  SP,  resistivity,  and  other 
geophysical  testing  and  also  considering  pre¬ 
vious  grout  takes,  fault  locations,  and 
piezometer  data,  WES  submitted  a  list  of 
proposed  locations  for  exploratory  borings  to 
SWL  for  approval.  Figure  3  shows  the  WES  sug¬ 
gested  exploratory  boring  locations. 

Extensive  investigations  were  conducted  on  each 
of  the  borings,  typically  included  soil  sam¬ 
pling,  diamond  core  drilling,  detailed  descrip¬ 
tive  logging  of  rock  core,  dye  testing  at  zones 
of  drill  fluid  loss,  pressure  testing  of  rock, 
downhole  geophysical  logging,  inspection  with 
downhole  video  equipment,  and  laboratory  testing 
of  rock  core  samples. 

The  investigations  conducted  in  the  exploratory 
borings  determined  that  the  northern  fault  zone 
has  a  vertical  offset  of  230  ft  while  the 
southern  fault  zone's  vertical  displacement 
measures  approximately  146  ft.  The  unsound  na¬ 
ture  of  the  fault  zones  was  evidenced  during 
drilling  by  noting  the  complete  loss  of  drill 
fluid  and  large  core  losses.  This  condition  was 
substantiated  by  SWO  laboratory  personnel  while 
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Figure  4 .  Integrated  seepage  nap 


performing  "down-looking"  and  "side-wall 
looking"  observations  with  a  down-hole  video 
camera.  Numerous  open  cavities,  channels, 
joints,  and  intensely  fractured  zones  were  en¬ 
countered  in  the  the  fault  zones  as  well  as  in 
the  upper  cherty  Boone  Formation.  Subsurface 
flows  through  channels  in  rock  were  apparent  in 
several  borings  where  normally  suspended  fines 
could  be  seen  moving  rapidly. 

Piezometers 

There  were  26  open  (well  point)  piezometers  at 
Dike  1  prior  to  the  seepage  investigation.  A 
review  and  analysis  of  locations  and  depths  of 
the  existing  piezometers  was  made  to  determine 
key  areas  (and  depths)  where  piezometric  data 
was  inadequate  for  analyzing  the  overall 
groundwater  (seepage)  flow  network  beneath  Dike 
1.  The  new  piezometers  were  located  in  a  direc¬ 
tional  alignment  pattern  (grid)  with  existing 
piezometers  to  facilitate  preparation  of  cross 
sections  through  the  piezometers,  both  parallel 
and  perpendicular  to  the  dike.  Piezometers  were 
also  located  to  give  broader  coverage 
(north/south)  of  the  fault  zones  and  downstream 
seepage  areas  (east/west).  The  new  piezometers 
were  dual-tipped  and  were  designed  such  that  the 
lover  tip  was  placed  at  or  near  elevation  1040 
ft,  which  is  within  ths  zone  between  knovn 
seepage  exits  and  Parshall  flumes  located  fur¬ 
ther  downstream,  and  the  upper  tip  placed  within 
the  upper  weathered  Boone  Formation,  at  an 
average  depth  of  11  feet  below  top-of-rock. 


Four  of  the  new  piezometer  sites  were  relocated 
in  May  1985  based  on  results  of  SP  and  resis¬ 
tivity  geophysical  tests  conducted  by  WES. 

Figure  3  shows  the  WES  proposed  piezometer  loca¬ 
tions.  Thirty  piezometers  were  installed  at  Dike 
1  between  the  period  May  and  September  1985, 
giving  a  total  of  56  piezometers  at  the  struc¬ 
ture  . 

The  piezometer  borings  drilled  in  1985  at  Dike  1 
were  sampled  and  tested  to  determine  subsurface 
conditions  prior  to  installing  the  piezometers. 

A  common  difficulty  was  heavy  loss  of  drilling 
(circulation)  fluid,  with  most  borings  having  a 
total  circulation  loss  at  some  point  during 
drilling.  A  downhole  camera  lowered  into  several 
of  the  piezometers  in  August  1985  indicated  rock 
characteristics  and  features  which  contribute  to 
subsurface  seepage  such  as  open  cavities,  chan¬ 
nels,  intensive  fracturing,  and  weathering. 

Seepage  Flow  Measurements 

Prior  to  the  seepage  investigation  there  was 
only  one  Parshall  flume  used  for  measuring 
seepage  flow  rates  downstream  of  Dike  1.  The 
frequency  of  the  flow  measurements  were  taken 
based  on  pool  level.  Measurements  were  made  by 
reading  the  water  level  on  a  scaled  gauge  on  the 
interior  wall  of  the  flume  and  converting  the 
readings  to  gpm.  In  October  1985,  a  flow  re¬ 
corder  was  installed  on  the  Parshall  flume  al¬ 
lowing  the  seepage  to  be  monitored  continuously. 
In  November  1985,  a  second  Parshall  flume  and 
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recorder  were  installed  approximately  170  ft 
downstream  of  the  first  flume.  The  necessity  for 
a  flume  at  this  second  location  came  from  the 
appearance  of  the  new  muddy  seepage  flow  which 
bypassed  the  first  flume. 

Topographic  Surveys  and  Mapping 

Field  control  for  the  seepage  investigations, 
boring  locations,  and  geophysical  surveys  was 
established  by  installing  a  200-foot  survey 
point  grid.  Also,  the  topographic  map  of  Dike  1 
was  updated  by  a  new  planetable  survey,  with  a 
final  plan  on  2-foot  contour  lines. 

Seepage -Siudv  Findings  and  Recommendations 

In  April  1986  SUL  reported  the  findings  of  the 
seepage  study  in  Supplement  No.  1  of  the  Recon¬ 
naissance  Report.  The  report  concluded  that  the 
foundation  beneath  Dike  1  was  in  an  advanced 
stage  of  deterioration,  and  that  seepage  could 
be  generally  described  as  pervasive.  Also,  the 
risk  factor  and  potential  existed  during  a  high 
(flood)  pool  condition  for  one  of  the  numerous 
seepage  flows  to  seek  a  new  and  larger  exit 
path,  by  removal  of  detrital  material  (cavity 
clays,  etc.),  and  "blowout"  through  the  overbur¬ 
den  in  the  downstream  area.  Finally  the  report 
concluded  that  the  element  of  time  was  both  a 
critical  and  a  debatable  factor  on  a  seepage 
problem  such  as  this. 

The  report  also  investigated  various  alterna¬ 
tives  for  controlling  seepage  beneath  the  dike. 
The  seepage  control  alternatives  considered  at 
Dike  1  ware  construction  of  an  additional  grout 
curtain,  a  cutoff  wall,  a  downstream  berm, 
placement  of  an  upstream  blanket,  or  do  nothing 
and  continue  monitoring  the  seepage.  The  recom¬ 
mended  remedial  technique  was  a  concrete  diaph¬ 
ragm  cutoff  wall  installed  upstream  of  the  cen¬ 
terline  of  the  dike.  The  report  concluded  that 
this  was  the  most  feasible  method  to  adequately 
provide  a  positive  cutoff  of  the  seepage.  The 
other  methods  were  considered  to  be  only  tem¬ 
porary  measures  to  control  seepage  and  in¬ 
adequate  for  providing  a  positive  cutoff,  which 
was  deemed  necessary  from  seepage  investiga¬ 
tions  . 

Automated-Piezometer  System 

During  the  period  November  1986  through  April 
1987  an  automated  monitoring  network  was  in¬ 
stalled  at  Beaver  Dam  to  read  all  (88)  open  well 
point  piezometers  at  Dike  1  and  the  main  embank¬ 
ment  and  Dike  3.  The  system  transmits  the 
piezometer  information  via  telephone  modem  to 
the  District  Office  (250  miles).  Readings  are 
routinely  taken  every  «  hours  and  can  be  read 
with  a  higher  frequency  if  needed.  Since  the  in¬ 
stallation  of  the  automated  system  a  high  degree 
of  interconnection  between  piezometers  has  been 
detected.  This  was  evidenced  during  periods  of 
drilling  or  performing  down-hole  tests  when 
piezometers  were  being  monitored  at  intervals  as 
short  as  1  hour.  The  automated  system  should  aid 
in  constructing  more  accurate  piezometric 
profiles  of  the  site  since  short  term 
piezometric  head  versus  pool  level  can  be  deter¬ 
mined. 


PLAN  FOR  REHABILITATION 

In  September  1987  a  Feature  Design  Memorandum 
was  prepared  by  SWL.  The  report  described  the 
recommended  design  for  a  concrete  diaphragm 
cutoff  wall. The  plan  for  the  wall  consists  of 
constructing  the  wall  through  the  embankment  and 
permeable  zone  of  the  foundation  rock.  The  wall 
will  be  a  minimum  of  1,600  feet  long,  2  feet 
wide,  and  vary  in  depth  from  130  to  205  feet. 

The  estimated  cost  of  constructing  the  cutoff 
wall  is  $16,000,000. 

A  rock-mill  type  excavation  system  will  be  used 
to  excavate  the  cutoff  wall  trench,  using  ben¬ 
tonite  slurry  to  stabilize  the  trench  during 
both  excavation  and  concrete  placement.  The 
rock-mill  was  determined  to  be  the  most  effi¬ 
cient  and  cost  effective  method  to  construct  the 
proposed  wall  due  to  the  amount  and  characteris¬ 
tics  of  rock  that  will  be  encountered.  More 
detailed  information  on  this  excavation  method 
is  given  by  Hess,  1985. 

Also  included  in  the  Feature  Design  Memorandum 
was  a  recommendation  by  WES  to  install  and  main¬ 
tain  an  automated  geophysical  monitoring  network 
to  monitor  seepage  before,  during,  and  after  im¬ 
plementation  of  a  remedial  measure  (such  as  a 
concrete  cutoff  wall).  The  result  of  the 
monitoring  network  analysis  will  be  an  assess¬ 
ment  of  the  effectiveness  of  the  remedial 
measure.  The  computer  controlled  network  is  en¬ 
visioned  as  consisting  of  a  permanently  in¬ 
stalled  SP  array  and  borehole  resistivity  probes 
with  the  capability  of  scanning  the  network  at 
any  desired  time  interval. 


SUMMARY  AND  CONCLUSIONS 

The  results  of  the  seepage  investigation  indi¬ 
cated  that  the  foundation  beneath  Dike  1  was  in 
an  advanced  stage  of  deterioration,  and  that 
seepage  can  be  generally  described  as  pervasive. 
Also,  the  risk  factor  and  potential  exists 
during  a  high  (flood)  pool  condition  for  one  of 
the  numerous  seepage  flows  to  "blowout"  through 
the  overburden  in  the  downstream  area.  Finally 
the  report  concluded  that  the  element  of  time  is 
both  a  critical  and  debatable  factor  on  a 
seepage  problem  such  as  the  one  above. 

The  investigation  also  recommended  that  a  con¬ 
crete  diaphragm  wall  be  installed  upstream  of 
Dike  1  as  a  mean  of  controlling  seepage. 

By  conducting  a  comprehensive  seepage  program 
such  as  the  one  performed  at  Beaver  Dam  it  has 
been  demonstrated  that  integration  of  results 
from  various  phases  of  the  investigation  has  led 
to  a  more  rational  approach  to  remedial  seepage 
planning.  In  a  program  of  this  magnitude  it  is 
very  important  to  consider  the  geophysical  sur¬ 
veys  as  an  integral  part  of  the  seepage 
analysis.  It  is  also  important  for  the  project 
engineer  and  the  geophysicist  to  communicate 
with  each  other  and  share  their  knowledge  of  the 
project  in  order  to  make  more  meaningful  inter¬ 
pretations  of  test  data  and  to  more  efficiently 
plan  any  future  testing. 
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